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Figure 19-12 

Stages of resonance, (a) The 
first compression meets the 
fork and gives it a tiny and 
momentary push; (b) the 
fork bends and then 
(c) returns to its initial posi- 
tion just at the time a rare- 
faction arrives and 
id) overshoots in the oppo- 
site direction. Just when it 
returns to its initial position 
(e) the next compression 
arrives to repeat the cycle. 
Now it bends farther 
because it is moving. 


If the forks are not adjusted for matched frequencies, the timing of 
pushes is off, and resonance will not occur. When you tune your radio 
set, you are similarly adjusting the natural frequency of the electronics 
in the set to match one of the many surrounding signals. The set then 
resonates to one station at a time instead of playing all stations at once. 

Resonance is not restricted to wave motion. It occurs whenever suc- 
cessive impulses are applied to a vibrating object in rhythm with its nat- 
ural frequency. In 1831, cavalry troops marching across a footbridge near 
Manchester, England, inadvertently caused the bridge to collapse when 
they marched in rhythm with the bridge’s natural frequency. Since then, 
it has become customary to order troops to “break step” when crossing 
bridges. A more recent bridge disaster was caused by wind-generated 
resonance (Figure 19-13). 

The effects of resonance are all about us. Resonance underscores not 
only the sound of music, but the color of autumn leaves, the height of 
ocean tides, the operation of lasers, and a vast multitude of phenomena 
that add to the beauty of the world about us. 



Interference 


Sound waves, like any waves, can be made to exhibit interference. Recall 
that wave interference was discussed in the last chapter. A comparison 
of interference for transverse waves and longitudinal waves is shown in 
Figure 19-14. In either case, when the crests of one wave overlap the 
crests of another wave, increased amplitude results. Or when the crest 
of one wave overlaps the trough of another wave, decreased amplitude 
results. In the case of sound, the crest of a wave corresponds to a 
compression, and the trough of a wave corresponds to a rarefaction. 
Interference occurs for both transverse and longitudinal waves. 

An interesting case of sound interference is illustrated in Figure 
19-15. If you are an equal distance from two sound speakers that emit 
identical tones of constant frequency, the sound is louder because the 
effects of the two speakers add. The compressions and rarefactions ol 
the tones arrive in step, or in phase. However, if you move to the side 
so that the paths from the speakers to you differ by a half wavelength, 
then the rarefactions from one speaker will be filled in by the compres- 


2 



The superposition of two identical transverse waves in phase produces a wave of increased amplitude. 


Figure 19-14 

Wave interference for trans- 
verse and longitudinal 
waves. 




The superposition of two identical longitudinal waves in phase produces a wave of increased intensity. 



Two identical trans- erse waves that are out of phase destroy each other when they are superimposed. 


d 





Two identical longitudinal waves that are out of phase destroy each other when they are superimposed. 


sions from the other speaker. This is destructive interference. It is jusi 
as if the crest of one water wave exactly filled in the trough of another 
water wave. If the region is devoid of any reflecting surfaces, little or 
no sound will be heard! 


Predicting Flow-Induced 
Vibrations in a Convoluted 
Hose 

A composite model was constructed from two 
less accurate models. 

Marshall Space Flight Center, Alabama 

A composite mathematical model predicts approximately 
the frequencies and modes of vibrations induced by flows of 
various fluids in a convoluted hose (see figure), This model was 
constructed from two older approximate models that were 
complementary to each other. One of the older models did not 
adequately account for interactions between the hose struc- 
tures and the fluids and therefore gave accurate predictions 
only in cases in which the fluids had little effect or in which the 
fluids acted as simple additional mass loads that were parts of 
the hose structures. The other older model accounted for 
fluid/structure interactions but was applicable only to hoses 
made of steel. 

Like the second-mentioned older model, the composite 
model is based partly on a spring-and-lumped-mass repre- 
sentation of dynamics that involve the springiness and mass of 
a convolution of the hose and the density of the fluid in the 
-hose. The Laplace-transformed equations of coupled motion 
of the lumped masses and springs lead to the following fourth- 
order equation for the complex frequencies of the vibrational 
modes: .''4 

In this equation, s is the Laplace-transform complex fre- 
quency, K is the spring stiffness of a convolution, rf is the 
mass density of the fluid, A is a cross-sectional area of a con-_ 


volution, C is called the “fluid capacitance” (analogous to elec-' 
trical capacitance) and is a measure of the energy and 
momentum associated with the change in the amount of fluid 
in a convolution during flexing of the hose, and L is called “fluid 
inertance” (analogous to both inertia and electrical inductance) 
and is a function of the dimensions of a convolution. The solu- 
tion of this equation consists of two sets of complex roots 
associated with two frequencies. The resonant frequencies 
predicted by this model are generally lower than those 
observed in experiments. 

This work was done by Stuart A. Harvey of Rockwell 
International Corp. for Marshall Space Flight Center. For 
further information, write in 101 on the TSP Request Card. 
MFS-29969 


Vibrations 
In Both 
Stretching 
and 
Bending 
modes can 
occur when 
the flow of 
fluid in the 
hose inter- 
acts with the 
convolu- 
tions. 
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HARMONIC VIBRATIONS. 


Experimental Science, Vol. 1 , George M. Hopkins, 25 th ed, NY: Munn & Co, 1906 
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in Fig. 173 consists of a pair of tubes made to slide 
telescopically one within the other, the inner one being 
graduated to indicate the different lengths required for dif- 
ferent pitches, and consequently for different speeds. As 
the fan revolves, the air drawn in through the holes at the 
center of the casing is thrown ontward by centrifugal force, 
thus maintaining a pressure of air at the periphery of the 
fan. The holes in the rim of the fan allow the air to escape 
in regular puffs, the frequency of which depends upon the 
velocity of the fan. These puffs produce sounds varying in 
pitch and intensity with the speed of the fan, and the reso- 
nating tube re-enforces the particular note to which it is 
tuned, so that when a speed is reached corresponding with 
the adjustment of the tube, the fact is known by the 
superior strength of that particular note. Any change 
of speed may be detected by the lessening of the intensity 
■of the sound and the change of pitch. 

The siren is shown in Fig. 175 in connection with me- 
chanism for driving it by hand. It is provided with a rev- 
olution counter and with a trumpet-shaped resonator. It is 
designed to be used in the same manner as the siren of 
Cagniard Latour, and, like that instrument, it yields sounds 
under water. 
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Experimental Investigation of a Cylindrical Resonator 

James H. T. Wu,* P. Perry Ostrowski, 4 Rafik A. Neemeh^ and Patrick H. W. Lee§ 

McGill University, Montreal, Canada 


Similar to a conventional resonance tube, a cylindrical resonator is investigated experimentally for producing 
high pressure and temperature gases. A quartz pressure transducer was used to determine pressure at the center 
of the resonator disk and a schlieren system to observe internal and external shock waves. A peak pressure 
nearly 3 times the jet stagnation pressure is observed and the internal shock waves are shown to be somewhat 
nonsymmetric and eccentric. 


Introduction 

S T! 'Y of fluid oscillations associated with the operation of 
a . ;ple resonance tube was initiated by Hartmann. 1 This 
device consists of a closed-end tube mounted with its open end 
facing an underexpanded fluid jet. Spark schlieren 2 and shadow 
photography 3,4 has shown that the resonant cycle is character- 
ized by alternate periods of inflow and outflow from the tube 
which correspond to the appearance of an internal travelling 
shock wave and expansion waves. Thompson 4,5 has introduced 
a wave diagram analysis based on appropriate boundary con- 
ditions at the tube mouth. A more recent investigation 6 has 
improved Thompson’s theory and used streak schlieren photo- 
graphs to demonstrate that the traveling shock wave is internally 
generated. 

It is well-known that the structure of the exciting jet strongly 
affects the performance of a simple resonance tube. 3 At a pressure 
ratio less than 4, underexpanded choked jets have a repetitive 
cellular appearance and resonance occurs only if the tube mouth 
is placed in the downstream half of each jet cell. This half cell 
is termed “the intervals of instability". Smith and Powell 3 have 
shown that higher harmonic resonant frequencies are obtained 
with the cavity mouth located near the midpoint of a jet cell and 
that the fundamental frequency is obtained at the downstream 
end. 

Current interest 7,8 in aerodynamic resonators is focused on 
their use as a high-temperature device or igniter by taking 
advantage of the repeated irreversible shock heating of the in- 
digenous fluid which remains permanently trapped inside a 
straight or tapered tube. It is noted that end wall temperatures 
of the order of 400°C can be achieved with gaseous nitrogen 7 
and ! 200 ’C with helium. 8 

The present investigation demonstrates a working device based 
on a combination of principles given by Hartmann 1 and 
GuderleyC Amplification of the internal shock wave can be 
accomplished by employing a cylindrical resonator cavity ex- 
cited by a circumferential fluid jet. The resonant cycle is then 
characterized by a cylindrical implosion-explosion process which 
has a significantly greater potential for generation of high tem- 
perature and pressure than the conventional plane shock wave. 
Since all significant shock amplification occurs in a com- 
paratively small region surrounding the implosion center, the 
cylindrical symmetry of the imploding wave is crucial to the 
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successful operation of this device. Accordingly, the aim of this 
investigation is to examine the operating characteristics of such 
a cylindrical configuration. To this end, spark schlieren and 
shadow photography is employed to observe the nature of the 
external flowfield (exciting jets) as well as the internal flowfield 
(traveling shock wave). The resonant frequency and pressure 
amplitude are measured with a pressure transducer located at the 
geometric center of one resonator disk. Acoustic theory is used 
for a rough estimate of the resonant frequency. 


Acoustic Theory 

Despite the fact that strong traveling shocks may appear inside 
the resonator cavity, acoustic theory may be employed to give 
an estimate of the resonant frequency. Essentially, the internal 
waves are taken to be acoustic in nature so that the wave motion 
is then described by the linear wave equation. In cylindrical co- 
ordinates this procedure yields the well-known solution for the 
perturbation velocity potential, 10 which is written in terms of the 
familiar Bessel functions. Application oT the boundary con- 
ditions (the particle velocity at the center must be zero and the 
perturbation pressure is zero at the resonator mouth i.e., r = r 0 ) 
yields the condition Jo(2nfrJc) = 0 where / is the frequency 
and c the local sound speed. The first root of the Bessel 
function then gives the fundamental frequency of acoustic 
vibration: / = c/2.61 5r 0 . 


Experiments 

A cross-sectional view of the cylindrical resonator is presented 
in Fig. 1. Two circular disks 2.65 in. in diameter with a lip 
chamfer of 45° are separated by a gap to form the resonator 
cavity. These disks are supported by a hollow (for schlieren 
viewing) threaded rod which is held in place by a threaded 
tightening ring so that the disk gap is adjustable. The resonator 
cavity is surrounded by a cylindrical plenum chamber 7.5 in. in 
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Fig. 1 Schematic of the cylindrical resonator. 




Fig. 2 A typical spark shadowgraph showing jet cells (pressure 

ratio = 3.94). 

diameter fitted with five compressed air inlets connected to a 
common supply. A narrow slit 0.030 in. wide with a \ in. 
rounded entrance in the inner wall of the plenum serves as a 
nozzle for the exciting jet which is directed radially towards the 
resonator cavity. The device ss constructed entirely of aluminum 
except for the central portion of the resonator disks which are 
fitted with Plexiglas windows for schlieren viewing. 

The most important physical parameters which affect the 
acoustic or thermal output of the resonator are the nozzle width 
and diameter, jet stagnation pressure, resonator disk gap and 
diameter, and the resonator material. The present study is largely 
a qualitative one and an extensive parametric program to maxi- 
mize the resonator output has not been undertaken. Instead, a 
nozzle diameter of 2.95 in. and nozzle width of 0.030 in. were 
arbitrarily chosen and resonator disks of 2.75 and 2.65 in. were 
tested separately. The latter were found to yield a consistently 
higher resonant pressure amplitude for all jet pressures. 
Similarly, a resonator disk gap of 0.068 in. was also found to 
be superior to any other value in this regard. These values were 
then held fixed for the balance of the tests despite the fact that 
they were not determined from an exhaustive parametric study. 

With the resonator geometry thus fixed, the test procedure 
was largely concerned with adjusting the jet stagnation pressure. 
As the pressure is increased, each jet cel! is lengthened and the 
effective location of the resonator mouth is then altered. For the 
present investigation, three series of tests were performed. First, 
flash shadowgraphy was employed to examine the jet cell 
dimensions with the resonator disks removed. Second, a 5 in. 
quartz pressure transducer (peb Piezotronics Inc., Model 101 A) 
was mounted in the center of one resonator disk so that the 
pressure history could be recorded on an oscilloscope (Tektronix 
Model 555). Finally, the resonator disks were fitted with Plexi- 
glas windows and the resonant cycle was examined via spark 
schlieren photography. The w'indows gave a 1.5-in.-diam field 
of view at the geometric center of the resonator cavity and the 
optical axis was normal to the resonator disks. 

Results and Discussion 

The structure of the exciting jet alone (resonator disks re- 
moved) is shown in Fig. 2 for a stagnation pressure ratio of 
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Fig. 3 Jet cell radius vs jet pressure ratio (nozzle radius = 1.475 in. 
and nozzle width = 0.03 in.). 

3.94. In this view, three jet cells are visible with the first exhibiting 
the highest degree of circular symmetry. Since cel! symmetry is 
crucial to the operation of this device, it appears to be best to 
design the resonator to operate on the first or outermost cell of 
the exciting jet. 

Figure 3 shows the jet cell measurements taken from several 
shadowgraphs similar to that of Fig. 2 for jet pressure ratios 
between 3.0 and 4.5. It can be seen that cell length increases 
for higher pressure ratios and the second cell is longer than 
the first. Now, with the resonator disk in place, resonance 
is observed to commence at a pressure ratio of 3.76. When this 
operating point is superimposed on the figure, it surprisingly falls 
into the first half of the second jet cell which is theoretically 
outside the “intervals of instability." Yet resonance persists. The 
possible reason for this is that the jet structure is altered when the 
resonator disks are in place, as the disks are chamfered with 
some bluntness. This type of behavior has been observed for 
simple resonance tubes as well . 3 

The pressure history at the center of one resonator disk is 
presented in Fig. 4 for jet pressure ratios between 3.69 and 5.19. 
Resonance begins very abruptly at a pressure ratio of 3.76 and 
thereafter the pressure amplitude (T,,,.,, P„ u „) falls rapid ly and 
becomes only about one fourth its initial level at a pressure ratio 
of 5.19, as in Fig. 5. At the onset of resonance, the pressure 
trace is seen to exhibit periodically a very steep rise which 
implies the existence of a strong internal imploding-exploding 
shock wave. The peak pressure of approximately 163 psia is 
nearly three times the jet stagnation pressure of 55 psia. In 
addition, the subsequent expansion phase of the resonant cycle 
is seen to take place in three distinct stages because of wave 



Fig. 4 Pressure measurements at the disk center. 
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Fig. 5 Pressure amplitude vs jet pressur ratio (3.76 is a minimum value 

for resona r •). 

interaction within the confined •„ verging region. As the jet 
pressure ratio is increased, the waveform gradually becomes dis- 
torted until eventually the motion appears to be entirely acoustic. 
It is implied from Fig. 3 that this behavior is due to lengthening 
of the jet cell structure so that eventually the resonator mouth 
is effectively located in the upstream half of the cell and 
resonance cannot be maintained in this stable interval. 

According to the oscilloscope traces, the resonant frequency 
is 3.33 kcps, which is about 15% lower than the acoustic 
frequency of 3.93 kcps with c= 1135 fps and r„ — 1.325 in. 
Attempting to account for shock heating effects by raising c 
above its ambient value makes the discrepancy even larger. Thus 
it appears that the shock-wave speeds with respect to the disk 
are subsonic over a considerable segment of the resonant cycle. 
Not surprisingly, the frequency measured from Fig. 4f is only 
4% lower than the theoretical acoustic value (see Fig. 5). 

Schlieren photographs, taken randomly but arranged in time 
sequence to illustrate pertinent features of the resonant cycle for 
a jet pressure ratio of 3.76, are presented in Fig. 6. These clearly 
show a nearly circular internal shock wave. Although whether 
this wave is imploding or exploding has not been determined 
experimentally, Figs. 6b and 6c are assumed to show the im- 
ploding wave as the enclosed flowfield appears quiescent, while 
Figs. 6d and 6e are thought to show the exploding shock for the 
opposite reason. Figures 6a and 6d show in- and outflow of 
the resonant cycle as the external shock is pushed far upstream 
against the impulse of the exciting jet, respectively. 

Figure 6 also yields some qualitative information regarding the 
resonator performance. The internal shock, although nearly cir- 
cular, is not perfectly so and is characterized by the appearance 
ol several Mach reflections 11 around its periphery. The shock 
wave does not appear to collapse to the geometric center of the 
disks, therefore it seems likely that the resonator performance 
(i.e.. pressure amplitude or temperature) can be increased by re- 
fining the design. It is noted that the appearance of a burn close 
to the center of the Plexiglas windows (small dark area in the 
photographs) indicates that the internal shock creates high tem- 
perature and collapses into a fairly small region. 

Conclusions 

A cylindrical resonator operating on an imploding-exploding 
cycle has been demonstrated experimentally. To achieve resonant 
conditions, the resonator disks must be precisely aligned relative 
to the excitingjet nozzle and the disk gap must be adjusted until 
the pressure amplitude is maximized. This type of “tuning” is 
required because a systematic experimental program to define the 
effect of various geometric parameters such as nozzle width and 
radius, disk gap and radius, etc., has not been completed. 

The higher acoustic potential of the present device over con- 
ventional resonance tubes is demonstrated as a pressure rise of 



Fig. 6 Schlieren photographs showing shock w'aves (jet pressure 
ratio = 3.76 and total pressure = 55 psia). 

nearly 3 times the jet stagnation pressure has been measured 
at the center of the resonator cavity. Acoustic theory yields a 
fair estimate of the resonant frequency as the measured value is 
about 15% lower. 

Shadow and schlieren photographs show that 1) without the 
resonator disks, the jet from the circular nozzle has clearly de- 
fined jet cells of nearly circular shape, and 2) with the disks in 
place, the imploding and exploding shock waves in the central 
region of the disks are somewhat eccentric and characterized by 
the appearance of Mach reflections along the main shock profile. 
Although the cylindrical resonator is ultimately intended for high 
temperature generation, temperature measurements have been 
delayed until the shock wave eccentricity is improved. 
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Aerothermodynamic Characteristics 
of a Resonance Tube 
Driven by a Subsonic Jet 

C. E. G. Przirembel* 

Rutgers University, New Brunswick, N.J. 
and 

L. S. Fletcher! 

University of Virginia, Charlottesville, Va. 

Introduction 

A RESONANCE tube, in its simplest form, is a cylindrical 
or rectangular tube closed at one end and facing a 
coaxial fluid jet at the open end. If appropriate geometric and 
fluid conditions are satisfied by the jet/tube combination, the 
fluid within the tube is excited to violent harmonic motion. If 
the tube is sufficiently long, the intense periodic oscillations 
give rise to strong thermal effects which result in resonance 
tube endwall temperatures significantly higher than the 
corresponding jet stagnation temperatures. For subsonic and 
low Mach number supersonic jets, the fundamental resonance 
frequency is independent of the jet Mach number Af , and can 
be approximated by /= C/4L, where C is the speed of sound 
of the jet and L is the length of the resonance tube. 

Previous experimental and analytical work, reviewed by 
Przirembel and Fletcher,' was concerned primarily with 
resonance tubes driven by supersonic and underexpanded, 
choked jets. For subsonic jets, only very limited experimental 
data have been reported. Savory 3 has observed fluid 
oscillations for subsonic jets, when the resonance tube is 
supported axially by a cylindrical rod passing through the 
nozzle. The fluid dynamics of this flowfield have been ana- 
lyzed by Broucher, Maresca, and Bournay. 3 Sprenger 4 was 
able to obtain oscillating flows with jet Mach numbers as low 
as 0.52 by placing a thin thread across the jet axis. 
Vrebalovich 5 observed oscillations in a cavity placed in a 
subsonic flowfield by locating a ring trip upstream of the tube 
mouth. More recently Przirembel, Fletcher, and Wolf 6 have 
observed resonance in two-dimensional resonance tubes 
driven by a subsonic, two-dimensional jet, without the in- 
troduction of any tripping device. 

Experimental Program 

A simple schematic of the axisymmetric nozzle/resonance 
tube model and pertinent geometrical parameters are shown 
in Fig. 1. The air jet was generated by a simple, converging 
nozzle with a contraction ratio of 3 to 1. The nozzle exit 
diameter was 0.375 in. (0.953 cm). 

The resonance tube had an inside diameter d of 0.375 in. 
(0.953 cm). The length L to diameter d ratio was 30 for all 
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tests. The resonance tube endwall was designed and con- 
structed to accommodate iron-constantan thermocouples for 
temperature measurements, and an orifice for pressure 
measurements. 

Results and Discussion 

The variable test parameters were the ratio of the 
separation distance between the nozzle exit and the resonance 
tube inlet to the nozzle exit diameter S/D, the ratio of the jei 
stagnation pressure to the ambient pressure P 0 /P a , the ratio 
of the diameter of the trip to the nozzle exit diameter b/D, 
and the separation distance of the trip wire from the nozzle 
exit x. The nominal jet stagnation temperature was the at- 
mospheric temperature. 

Of principal interest in this experimental investigation was 
the determination of the temperature and thermal effects in 
the vicinity -of the resonance tube endwall. Figures 2 and 3 
summarize some typical results for jet flows with and without 
a trip. From Fig. 2, it can be seen that there is no measurable 
increase in the endw'all temperature above the jet stagnation 
temperature for a simple subsonic jet (A/ JC , =0.73). However, 
as a thin wire is placed across the center of the jet exit plane, 
there is a significant temperature increase. As the jet Mach 
number is increased in Fig. 3, there is some limited resonance 
tube heating w'ithout the trip. The shape of curves for flow in 
the presence of the trip are quite similar. The peak endwall 
temperature increases as the jet Mach number increases. 
These results are consistent with those reported by Sprenger. 4 
Broucher et al. 3 have argued that the presence of the trip 
results in a stagnation pressure deficiency in the jet, which 
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Fig. 3 Endwall temperature variation with S/D for Pq/P^ ~ 1-90- 


allows the resonance tube outflow to penetrate the jet. From 
wake calculations for the trip immersed in the core of the jet, 
it can be shown that the maximum velocity defect is on the 
order of 20% for the maximum endwall temperatures. At the 
inlet plane for the resonance tube, the trip wake width is on 
the order of one-half the resonance tube diameter. By 
comparing the frequency of oscillation (286 Hz) and the 
shedding frequency (73,400 Hz) associated with the trip, there 
is no indication that the resonance condition is driven by the 
wake shedding process. The endwall pressure measurements 
showed a significant correlation with the endwall temperature 
curves, i.e., pressure increases occurred at the same S/D 
values as temperature increases. 

The effect of the relative location of the trip wire x is shown 
in Fig. 4. As can be seen, there is a significant decrease in the 
maximum endwall temperature, as well as a shift in S/D as 
the wire is moved away from the jet exit plane. The jet cen- 
terline velocity defect increases very slightly at the resonance 
tube inlet. However, the mass flow of the jet entering the 
resonance tube decreases. 

In order to determine the influence of trip diameter on the 
maximum endwall temperature, a series of experiments were 
performed which resulted in curves for different trip 
diameters similar to those shown in Figs. 2-4. For these ex- 
periments, the trips were placed at the exit plane of the nozzle. 
A summary of these results is presented in Fig. 5, where the 
maximum endwall temperatures are plotted for various trip 
diameters. The primary observation is that there is no distinct 
peak for a particular combination of wire diameter and jet 
velocity. Vortex shedding frequencies based on a Strouhal 
number of 0.21 do not provide any additional information on 
the mechanism of resonance since the estimated peaks appear 
independent of jet velocity. It should be noted that the 
calculated shedding frequencies are on the order of 200 times 
the frequency of oscillation in the resonance tube. From these 
observations it can be concluded that the oscillations in the 
resonance tube are not directly driven by the vortices shed 
from the trip. 
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Fig. 5 Maximum endwall temperature variation with trip diameter. 


As indicated earlier, Broucher et al. 3 have suggested that 
the best method to obtain oscillations in a resonance tube 
driven by a subsonic jet is to decrease the kinetic energy of the 
jet in the neighborhood of the axis in order to act upon the 
radial pressure distribution at the resonance tube entrance. 
Without detailed jet velocity and static pressure 
measurements near the resonance tube, it is not clear at this 
point if the present experimental results support their 
analysis. The wake defect of the present cylinder is still 
significant for the separation distances at which maximum 
endwall temperatures were measured. However, it is also 
possible that the trip wake vortices cause disturbances in the 
outer region of the jet and thereby initiate instabilities which 
ultimately lead to oscillations in the resonance tube. 
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peratures well in excess of the exciting jet stagnation tem- 
perature could be achieved with a uniform tube. Although it 
was suspected that this temperature rise resulted from 
nonlinear wave motion, this was not verified until spark 
shadowgraphs 2 clearly showed the existence of internal 
traveling shock waves. Streak schlieren photography 3 has 
shown that these are generated during the inflow phase of 
each cycle and that when higher harmonics are driven, several 
traveling shocks are present simultaneously. 

Once it was recognized that thermal effects in the simple 
Hartmann-Sprenger tube resulted mainly from repeated 
shock heating of the fluid trapped inside, it became evident 
that the thermal performance could be enhanced either by 
strengthening the internal shock or increasing the frequency 
of the oscillations. This may be accomplished by a variety of 




Fig. 1 Spark schlieren photographs of 
the shock motion within the tapered 
cavity, P oj /P a = 1.92, /i/r/ = 0.85. 



Fig. 2 Spark schlieren photographs of 
the shock motion within the logarithmic 
spiral cavity, P 0J /P a = 1.92, h/ii= 0.85. 


Thermal Performance of a 
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Introduction 

T HE simple Hartmann-Sprenger tube consists of a closed- 
end tube that is excited by an underexpanded gas jet 
facing its open end to produce violent, cyclic fluid oscillations 
internally and externally. Although this phenomenon has 
been known for some time, the heat generation capability of 
this device has been exploited only relatively recently . 1 
Sprenger experimented to demonstrate that end-wall tem- 


inimizing the tube length, employing a 
light driving gas, or contouring the tube profiles to accelerate 
the traveling shock. Thus far, it appears that no attempt has 
been made to formally optimize the internal geometry in order 
to maximize the thermal performance; therefore, tube shapes 
hate been selected rather arbitrarily. Tapered tubes have 
commonly been employed to yield end-wall temperatures as 
high as 725°C using air 4 and 540°C with hydrogen. 3 A 
cylindrical resonator that operates on a shock implosion- 
explosion cycle has been devised by Wu et al. 6 in order to take 
advantage of the significant shock amplification that may be 
achieved through symmetrical implosion. 

The present Note investigates the thermal performance of a 
Hartmann-Sprenger tube with a two-dimensional logarithmic 
spiral interna! contour. This particular profile was selected on 
the basis of prior shock tube studies, 7 which have shown that 
the shock motion in such area contractions yields cylindrical 
implosion and is, therefore, attractive for the present ap- 
plication. For experimental convenience, only two-dim- 
ensional shapes are tested, and performance is judged on the 
basis of spark schlieren photographs, end-wall temperatures, 
and piezoelectric pressure measurements. 
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Experiments 

All tests were conducted with air as the working fluid. Flow 
visualization was accomplished with tubes of rectangular 
cross section (5.1x25.4 mm at the entrance) which were 
excited by an underexpanded jet from a converging nozzle 
with the same exit dimensions. Cavities with rectangular, 10- 
deg tapered, and logarithmic spiral wall contours with 



respective lengths of 15.2, 13.1, and 13.3 cm were tested. A 
single spark schlieren photograph was obtained for each test. 
Piezoelectric pressures were measured at different locations. 

Temperature measurements were conducted with a second 
set of tubes, 17.8-cm long and 10.2-mm square at the en- 
trance, whose walls were constructed with an asbestos 
compound. The tapered and spiral tubes have instrumentation 
pons 2.5-mm wide by 12.7-mm long for which chromelalumel 
thermocouples are located at the port midpoints. Each cavicy 
is excited by an underexpanded air jet of diameter d = 10.2 
mm. 


Results and Discussion 

Spark schlieren photographs of the shock propagation in 
the tapered and logarithmic spiral tubes are presented in Figs. 

1 and 2, respectively, for a jet pressure ratio of 1.92 and a 
nozzle to cavity separation hid- 0.85. Figure 1 is seen to be 
characterized by the appearance of Mach reflection due to the 
inclination of the channel walls, whereas Fig. 2 is observed to 
exhibit a curved shock normal to the wall surfaces with no 
reflected shocks trailing behind. Here, the shock appears to 
eventually collapse into a relatively small region. 

Piezoelectric pressure measurements were carried out near 
the focus of the logarithmic spiral channel. From the recorded 
pressure histories, the pressure rise across the incident shock is 
found to increase, while that ahead of the shock decreases as 
the shock proceeds toward the vertex of the spiral which yields 
significant shock amplification. The pressure rise across the 
shock becomes twice that observed in the parallel tube walls. 
Thus, it appears that the logarithmic spiral cavity exhibits 
much better shock amplification than the square cavity. 

The thermal performance of the three insulated tubes is 
shown in Fjg. 3, which presents the maximum measured end- 
wall temperature vs jet pressure ratio. It is evident from the 
figure that the logarithmic spiral profile consistently achieves 
a higher temperature than the tapered tube over the entire 
range of jet pressures employed for the present tests. The 
increment in thermal performance is roughly constant. This 
implies that the shock amplification in the logarithmic spiral 
tube is significantly greater than in the tapered and rec- 
tangular tubes. It is evident from Fig. 3 that the peak thermal 
performance of the logarithmic spiral tubes employed in the 
present tests is rather low (320° C) compared to that obtained 
by other investigators, 4 ' but it is possible that higher tem- 
peratures may be achieved using a longer tube length. 


300 
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200 


100 
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Jet Pressure Ratio - P 0 j/P a 

Fig. 3 Knd-wall temperatures. 



Tests were also conducted to determine those locations of 
the tube mouth relative to the nozzle that produced maximum 
end-wall temperatures. Results for the logarithmic spiral 
profile at the jet pressure ratio of 3.77 show that the 
maximum temperature varies periodically with the nozzle- 
tube separation h due to the periodic structure of the exciting 
jet. The highest temperatures were recorded at hid =1.7. This 
means that the peak corresponds to a location where the tube 
mouth is situated downstream of the first pressure cell. This 
result differs from the findings of Przirembel and Fletcher. 8 
However this difference might be explained by the fact that in 
their case, the open-end walls were rounded in contrast to the 
30-deg chamfered entrances of the present tests. 

Conclusion 

Thermocouple measurements have demonstrated that a 
Hart 1 mn-Sprenger tube with a logarithmic spiral profile 
achi , a significantly greater base temperature than a 10- 
deg ; ered tube or a rectangular tube of comparable length. 
This improved performance is interpreted to result from a 
higher degree of shock amplification which produces a greater 
entropy rise in the shock-heated indigenous gas that remains 
trapped inside the tube. Therefore, it appears that an 
axisymmetric logarithmic spiral tube could be used to improve 
the performance further. 
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the following expression for the modes of a vibrating string: 


where T is tension and /t is mass per unit length. 

If a string is driven at the frequency of one of its natural modes, 
resonance can occur. There are many ways in which to apply the driv- 
ing force, three of which are illustrated in Fig. 4.6. The magnetic drive 
shown in Fig. 4.6(b) works only for a string of steel or other magnetic 
material. The moving violin bow applies a rather complicated driving 
force (to be described in Chapter 10), which has components at several 
different frequencies. The string may also be driven elec- 
tromagnetically, even when the string is made of a nonmagnetic metal, 
by placing a permanent magnet near the string and passing an alter- 
nating current of the desired frequency through the string. 



4.4 PARTIALS, ! ,RMONICS, AND OVERTONES 




FIG. 4.6 

Three ways to drive a string at one of 
its resonances: (a) a tuning fork; (b) an 
electromagnet; (c) a violin bow. 


it is appropriate at this point to clarify nomenclature in order to avoid 
confusion. We will use the term harmonics to refer to modes of vibra- 
tion of a system that are whole-number multiples of the fundamental 
mode, and also to the sounds that they generate. (It is customary to 
stretch the definition a bit so that it includes modes that are nearly 
whole-number multiples of the fundamental: 2.005 times the fun- 
damental rather than 2, for example.) Thus we say that the modes of 
an ideal vibrating string are harmonics of the fundamental, but the 
modes of a real string are usually so close to being whole-number 


multiples that we also speak of them as harmonics. Note that the term 
“first harmonic” refers to the fundamental. 

Many vibrators do not have modes that are whole-number 
multiples of the fundamental frequency, however, and the term over- 
tones is used to denote theii higher modes of vibration. Harmonics are 
therefore described as overtones whose frequencies are whole-number 
multiples of the fundamental frequency. Minor confusion arises, 
however, from the fact that the term harmonics includes the fun- 
damental, but the term overtones does not. Thus the second harmonic 
is the first overtone, the third harmonic is the second overtone, and so 
forth. 

There is another term in common use that refers to modes of 
vibration of a system or the components of a sound: partials. Partials 
include all the modes or components, the fundamental plus all the 
overtones, whether they are harmonics or not. The term upper partials 
excludes the fundamental and thus is a synonym of overtones, but use 
of the former will be avoided in this book. 

The actual motion of a vibrating string is a combination of the 
various modes of vibration. The way in which these modes or partials 
combine is given by the spectrum of the vibration. A vibration spec- 
trum is like a recipe that specifies the relative amplitudes of the par- 
tials. Similarly, the spectrum of a sound specifies the amplitudes of its 
partials, as we will discuss in Chapter 7. 


4.5 OPEN AND CLOSED PIPES 

The reflection of sound pulses at open and closed ends of pipes was 
described in Section 3.5. At an open end, a pulse of positive pressure 
reflects back as a negative pulse; at a closed end, it reflects as a 



positive pulse. These two end conditions can be used to arrive at the 
resonances for open and dosed pipes. 

The motion of the vibrating air in a pipe is a little harder to 
visualize than the transverse vibrations of a string, since the motion of 
the air is longitudinal. The motion of the air is greatest at an open end, 
but the pressure variation is maximum at a closed end. A pressure- 
sensitive microphone inserted into the tube will pick up the most 
sound at the points where the pressure variations above and below at- 
mospheric pressure are maximum. Thus in Figs. 4.7 and 4.8, both the 
air motion and the pressure variations are shown. 

In an actual pipe, the pressure variations do not drop to zero right 
at the open end of the pipe, but rather a small distance beyond. Thus 
the pipe appears to have an acoustic length that is slightly greater than 
its physical length. For a cylindrical pipe of radius r, the additional 


length, called the end correction , is 0.61r. Twice this amount should 
be added to the length of a pipe with two open ends to obtain its 

acoustic length. 
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FIG. 4.7 

Modes of vibration or resonances of an 
open pipe. At the open ends the pressure 
is equal to atmospheric pressure. The 
resulting modes include both odd- 
numbered and even-numbered harmonics. 
Minimum motion occurs at the nodes 
denoted by /V. 


FIG. 4.8 

Modes of vibration or resonances of a 
closed pipe. At the closed end. the air 
motion is minimum but the pressure is 
maximum. The resulting modes include 
odd-numbered harmonics only. 
Minimum motion occurs at the nodes 
denoted by /V. 


(a) 




FIG. 4.9 

Resonance of a closed tube demonstrated by (a) holding a tuning fork near 
one end; (b) blowing on an organ pipe. 


Resonance of a tube can be demonstrated by placing a tuning fork 
near one end of it, as shown in Fig. 4.9(a). A piston at the closed end 
makes it possible to change the resonance frequency. Blowing through 
a closed (“stopped”) organ pipe may excite several of its resonances. 
Gentle blowing excites the lowest mode, but blowing much harder 
causes the pipe to vibrate in its first overtone, which for a closed pipe 
ls the third harmonic (/, = 3/.), a musical twelfth above the lowest 
mode. 

In Chapter 15, which deals with speech production, we will be in- 
terested in resonances of the human vocal tract that allow us to enun- 
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Pianos and harpsichords have large soundboards with mans 1 
resonances of their own closely spaced throughout the playing range. 
Vibrations are transmitted from the string to the soundboard through 
the bridge as in the violin or guitar. We will discuss pianos and harp- 
sichords in Chapter 14. 


4.9 SUMMARY 

Resonance occurs when a force applied to a vibrating system varies 
with a frequency at or near the natural frequency of the system. 
Linewidth, Q, and maximum response are ways to describe the sharp- 
ness of the resonance, which in turn depends on the amount of damp- 
ing in the system. The phase difference between the vibrating object 
and the driving force changes near the frequency resonance. 

Normal modes of vibration of strings, pipes, and similar systems 
may be described as standing waves. Standing waves, which consist of 
waves propagating back and forth between the boundaries, lead to 
resonances in vibrating systems. Acoustic impedance, which is the 
ratio of sound pressure to volume velocity, is a quantity useful in the 
analysis of acoustical systems. Since sound radiation depends on air 
displacement, radiation from a vibrating string is greatly enhanced by 


the sympathetic 
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acoustic impedance A measure of the difficulty of gen- 
erating flow (in a tube, for example); it is the ratio of the 
sound pressure to the volume velocity due to a sound wave. 

amplitude The height of a wave; the maximum displace- 
ment of a vibrating system from equilibrium. 

damping Energy loss in a system that slows it down or 
leads to a decrease in amplitude. 

electromagnetic force The force that results from the in- 
teraction of an alternating electric current with a magnetic 
field. 

overtone A mode of vibration (or component of a sound) 
with a frequency greater than the fundamental frequency. 

partial A mode of vibration (or component of a sound); in- 
cludes the fundamental plus the overtones. 

phase difference A measure of the relative positions ot two 
vibrating objects at a given time; also the relative positions, 
in a vibration cycle, of a vibrating object and a driving 
force. 

Q A parameter that denotes the sharpness of a resonance; 
Q =/ 0 /. A /, where / 0 is the resonance frequency and Sf is the 

linewidth. 


vibration of a soundboard. 
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fundamental The mode of vibration (or component of 
sound) with the lowest frequency. 

harmonic A mode of vibration (or a component of a 
sound) whose frequency is a whole-number multiple of the 
fundamental frequency. 

Helmholtz resonator A vibrator consisting of a volume of 
enclosed air with an open neck or port, 
linewidth The width -A/ of a resonance curve, usually 
measured at 71 percent of its maximum height; a measure 
of the sharpness of a resonance (a sharp resonance is 
characterized by a small linewidth). 

resonance When a vibrator is driven by a force that varies 
at a frequency at or near the natural frequency of the 
vibrator, a relatively large amplitude results. 

soundboard A sheet of wood or other material that 
radiates a substantial amount of sound when it is driven in 
sympathetic vibration by a vibrating string or in some other 
manner. 

spectrum A recipe for vibratory motion (or sound) that 
specifies the relative amplitudes of the part ials . 



Resonance in Closed -ended Pipes 
(written for this work by a retired mechanical engineer) 
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Considering only pressure changes-- that is, compressions, C, 
and rarefactions, R--in a closed pipe whose length, L, is 
adjusted to be one-fourth the wavelength of the entering sound 
(Fig. 1), let us follow successive regions of C and R as they 
move into the pipe and are reflected from the ends, remembering 
that change of pressure phase occurs at the open but not at the 
closed end. 

• ••••• 


l. 


Ri 


>C. 


<c. 


Fig. 1) 
resonance 
in 1/4- or 
3/4 -wave 
closed-end 
tubes 


R, 


R, 


5 . 




t. 


>RC< 


RR > 


<RR 


R, 


CCC> 


The compression wave C-^ enters the pipe, is reflected at the 
closed end, and reaches the open end just as R, is about to 
enter. As R^ enters, C. is reflected back in the opposite phase 
as an R, and unites with entering R^ . The now amplified R^ is 
reflected at the closed end, and on returning to the open end, 
is reflected as an amplified C, which is further amplified by 
uniting with entering C£ • 

Thus the pressure disturbance traveling back and forth within 
the pipe is greatly augmented, and an intensity is built up within 
the pipe which is very much greater than that in a corresponding 
length of the incoming wave itself. In this way, a sound of such 
low intensity as to be inaudible outside the pipe becomes loud 
and clear on reinforcement be the resonance built up in an air 
column properly tuned to respond to it. 

The same phenomenon occurs likewise if the closed pipe is 
three-quarters wavelengths of the incoming wave. Compression C-^ 
then unites with rarefaction R£ , and these again with Co and R^ ; 
also, C^ meets rarefaction R^ and C , . Similarly one gets reso- 
nance in a closed pipe if the length is any ODD number of quarter- 
wavelengths (see Fig. 1 again). a.(2n-l) 

Resonance, closed pipe, for L = 1, 3, 5, 7,.., 

or for a given pipe length, response will occur for waves of length 


*1, 3. 5, 7. 9 


= 4L , 4L, 4L, 4L, 4L , 

3 5 7 9 
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Computing Propagation of Sound in Engine Ducts 

An updated computer program is easier to use. 

Marshall Space Flight Center, Alabama 


The Frequency Domain Propagation 
Model (FREDOM) computer program 
accounts for acoustic loads applied to 
components of engines. In particular, 
the program models the propagation 
of noise through fluids in ducts 
between components and through 
passages within components. It can 
be used not only to analyze hardware 
problems, but also for design purpos- 
es. FREDOM is an updated version of 
the FREQPL program, which was first 
developed in 1959. 

These programs were devised 
specifically for use in analyzing the 
acoustic loads in rocket engines, 
where noise generated in pumps, 
combustion chambers, and duct 
bends can be so intense as to cause 
mechanical failures, and where reso- 
nances in ducts can amplify the 
effects of noise. The underlying physi- 
cal and mathematical concepts imple- 
mented in these programs should also 
be applicable to acoustic propagation 
in other enclosed spaces; they might 
be useful, for example, in analyzing 
process plumbing and ducts in indus- 
trial buildings with a view toward 
reducing noise in work areas. 

The older program, FREQPL, is 
basically a matrix solver with plotting 
capabilities. Before running FREQPL, 
the user had to do all the preparatory 
theoretical and computational work, 
which included (1) writing the often 
complicated wave equations, flow 
equations, and other equations as 
they apply to the system to be ana- 
lyzed, (2) putting the equations into 
matrix»vector form, and (3) entering 
the matrix elements into the computer 
by responding to program prompts 
one element at a time. This preparato- 
ry work required an engineer with 
a high level of sophistication in 
acoustics and detailed knowledge of 
the system, consumed much time, 
and was highly susceptible to theoret- 
ical and typing errors. Also, the matrix 
for FREQPL was set up to solve for a 
specific variable. Most of the time, 
results are needed at many locations 
(different variables) of the system. 
When using FREQPL, the preparatory 
work had to be repeated for each new 
variable. 

FREDOM incorporates the matrix 
solver of FREQPL, but it also auto- 
mates the preparatory theoretical 
work, thereby greatly reducing the 



Figure 1 . Various Boundary Conditions are treated automatically and routinely. FREDOM even pro- 
vides for small changes in cross-sectional area of a duct due to the distensibility of the duct material. 




Figure 2. The Frequency Response plotted here is a transfer function that represents the 
acoustic response at one point along a duct to an acoustic excitation at another point along the 
duct. The peaks in the gain plot indicate the resonances of the duct. Also note a customary 
phase change in the accompanying phase plot at the resonance peaks, 
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probability of errors, FREDOM is a 
user-friendly and menu-driven code. It 
guides the user step by step through 
the analysis. The user needs only a 
good knowledge of acoustics. 

The us r '’' r ' r ovides data that describe 
the due ’ oassage geometry, con- 
nections, and properties of the fluid in 
the system via four input files. FRE- 
DOM also accounts for boundary con- 
ditions of the system automatically. 
The user merely states the type(s) of 
boundary condition(s) and the code 
incorporates the coefficients of the 
corresponding equation(s) into the 
coefficient matrix. Figure 1 shows the 
variety of bound : conditions handled 

by FREDOM. Fen xample, if an accu- 
mulator is needed, the user need 
merely specify its location, the type of 
gas, and initial pressure and volume. 
For most boundary conditions, it suf- 
fices to specify locations only. 

FREDOM then constructs the matrix 
from the data provided by the user — 
that is, FREDOM automatically does 
the equivalent of what the user would 
do in writing the acoustic equations for 
the system, but does it faster and more 
accurately. The matrix in FREDOM dif- 
fers from that in FREQPL: the FRE- 
DOM matrix contains the coefficients 
of the system equations that have 
been manipulated to put all variables 
on the left side; that is, the right sides 
of all the equations are set equal to 
zero. This matrix, along with the infor- 
mation used to create it, is stored in a 
"system file”. Once the system file is 
obtained, any of the variables can be 
set on the right side of the eq* ,; ons 
for calculation of the applicable fre- 
quency response, without having to 
create a new matrix for each variable. 
Results can be plotted (see Figure 2) 
and/or tabulated. 

FREDOM works directly in the fre- 
quency domain. In this respect, FRE- 
DOM offers a great advantage over a 
time-domain code, the output data of 
which would have to be postprocessed 
to convert them into the frequency 
domain. Postprocessing could add 
several hours to the overall computa- 
tional process. In a large iterative analy- 
sis (for example, probabilistic analysis 
of a system), a frequency-domain 
approach takes several days less than 
does a time-domain approach. 

This work was done by Silvia Saylor 
of Rockwell International Corp. for 
Marshall Space Flight Center. For 
further information, write in 248 on 
the TSP Request Card. 

MFS- 29971 


How the right switch anatomy 
can help you avoid an autopsy 



Switch failure at the 
wrong time can have you sorting 
through the pieces later. That’s why Janco 
makes switches that work every time... even in the 
most extreme environmental conditions. 

Janco rotary switches are designed to resist it all: 
shock, vibration, thermal variance, moisture, salt, and other 
contaminates. So if you can’t afford the price of failure, 
contact Janco for more information on our extensive line of 
miniature, power, push-button, and solenoid operated rotary 
switches today. 

3111 Winona Avenue, Burbank, CA 91504 
Phone: (818) 846 -1800 • FAX: (818) 842-3396 


I W ANESOPCORP. 


For Mors Information Writo In No. 4S7 



NASA Tech Briefs, March 1995 
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Simplified Derivation of Water-Hammer Formula 

By Lewis F. Moody' Philadelphia , Pa. 


As an introduction to the extension of the principles of water hammer to various problems arising in 
hydraulic practice, it is thought that it may tie useful to review and to rederive the basic relations on 
which this branch of hydraulics rests. Moreover, the author must admit that he has found many of the con- 
ventional methods of derivation somewhat difficult to follow and not fully convincing, with their complete 
acceptance depending mainly on the knowledge that the results have been verified by experience, rather 
than on the logic of the derivation itself. For example, some authors of well-known textbooks^ base the 
proof of the Joukowsky law on the assumption that there is no loss of energy in impact or otherwise in the 
action which takes place during sudden closure of a conduit outlet. Such a premise is not convincing as 
it is certainly not axiomatic. The method which will be used here is based on the principle of impulse 
and reaction, or balance of forces, and rests directly on Newton's second law of motion, thus leading at 
ODce to exact relations not dependent on assumptions as to energy losses. It may be useful briefly to 
review the broader 1 field of varying flow before considering the elastic actions of water hammer specifical- 


i- Accelerated or Retarded Flow in Closed Conduits. Consider the case of flow which is uniform 
along the conduit, that is, constant from point to point, but varying with respect to time. The condi- 
tions of steady or unvarying flow are given by writing the Bernoulli formula between sections 1 and 2 for 
the pipe shown in the sketch, Fig. 1, using the customary notation, 

h pi + z i + h L = >> p2 + * 2 + 

2 g 2 g 

This is the condition for equilibrium; that is, the body of water between 1 and 2 is not being accelerated 
or retarded and no work is being stored in it or withdrawn from its store of inertia. But if b is 

suddenly reduced to h ' p2 below the value required by the equation, as by partially opening the valve shown, 

equilibrium is destroyed, and we have an unbalanced force acting on the mass of water tending to accelerate 
it. • This force is 6 


P - P 
r i 2 


* wA (h 


Pi 


h pi * 


wA lh ' p2 - h p2 ) 


and is equal to the mass times the acceleration M (dv/dt) or to (wAl/g) (dv/dt), A being the cross-sec 
tional area of the pipe. Hence (h' pl - h pl » - (h ' p2 - h p2 ) = (L/g) Idv/dt), and substituting the valu 


Lue of 


"pi 


or; h' 


p 2 


Pi 


from the first equation, we have 


h 


Pi 


P 2 


= h„, - h 


"pi 


_L dv 

P 2 g dt 


+ V 


2 g 


- (Z, + 


1 


- h, ) + JL dv 
L g dt 


Z 1 + - 


2 g 


head) at section 1 - lost 
section 2 


b L dv 
1 g dt 


P 2 


2 2 + 


head 


2g 


that is: (pressure head + elevation head + velocity 


acceleration head = 


(pressure head + elevation bead + velocity head) at 
This is the Bernoulli formula generalized to cover varying flow. Evidently, for retarded 
flow the sign of the acceleration head is reversed. If the conduit is uniform, v and v are identical- 
if of varying section, a mean value of v must be used in dv/dt, and the formula may be 2 readily extended 
to cover this more general case. The foregoing relation holds for slowly varying flow, when the "con- 
tinuity of flow relation Aj v } = A 2 v 2 applies throughout. But if the change of velocity is sudden, as 

when a valve is quickly closed in a pipe line, the elastic properties of the water and pipe walls affect 
the conditions and all of the water will not be uniformly retarded. Only the portions of the flow near 
the valve w i 1 T immediately feel the effect of the closure; these portions will be compressed and the 


he 


adjacent walls expanded by the pressure caused by the closure, and the increased pressure will require 
time to travel along the pipe before the distant portions of the flow feel the effect. Hence some of tu 
elements of the water column will be brought to rest, and other elements will continue at their original 
velocity until the region of compression has reached them; and we shall have different velocities at 
different portions of the pipe at the same instant. The region of compression will extend itself progress- 
iyely upstream at a uniform rate, since the pipe is assumed to be uniform, and this rate is the "velocity 
of the pressure wave.” J 

v 

2 - Sudden Closure-- Water Hammer. Consider the long uniform pipe shown in the same figure, flowing 
full with uniform velocity V. Then imagine the valve at its end to be instantaneously closed. ’ After a 

time t, the conditions will be as shown in the sketch. Fig. 2. A portion of the column of length x will 
have had its velocity reduced to zero, the column thus being divided into two regions, one having the 
original uniform velocity and the other being at rest. Before closure and at the instant of closure the 
momentum of the entire mass of water is MV. After the closure the entire mass is no longer moving with 
velocity V, but only a portion of it having a mass m. The increase F in the force on the end of the 
column caused by the closed gate is equal to the time rate of charge of momentum of the column (Newton's 
second law), or the rate at which momentum is destroyed; that is: F - d (mV) / dt. The rate of change of 
momentum will be the rate at which the amount of the mass moving with velocity V is reduced, or 


Professor of Hydraulic Engineering, Princeton University; and Consulting Engineer, Baldwin-Southwark Cor 
poratlon, I. P. Morris Division. Mem. A.S.M.E. 

example, Mansfield Merrlman'a "Treatise on Hydraulics," which for many years was a standard work on 
hydr aul lco, and to which the author was Indebted for his Initiation In the subject. 
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r * a <mvi /at » v (dm/dt) * product of velocity and rate at which mass is brought to rest 

Hence 

F = r VA w/g) (dx/dt) 


Now dx/dt is the rate at which the pressure wave travels along the pipe, or the velocity o f the pres- 
sure wave. This is the same as the velocity of sound in the water within the pipe. Call this velocity of 
propagation of the pressure wave a. 

Then F * aVAW/g = wAy, y being the rise of pressure head due to the closure. Hence the rise of 
pressure head due to complete instantaneous closure is y * aV/g. This is Joukowsky's law. 

Since the pressure wave is traveling up the pipe with uniform velocity, momentum is being destroyed 
at a uniform rate, and the pressure at the valve remains constant until the wave reaches the open end of 
the pipe. At any instant all the water between the valve and the wave front will be in a state of uniform 
compression under a rise in pressure head y and at zero velocity, and all the water from the wave front to 
the open end will be uncompressed and under its original pressure and moving with constant velocity V. 

When the wave front reaches the open end, we shall have the pipe completely filled with water at rest, 

under a uniform pressure head y feet above its initial value, and therefore y feet above the pressure in 

the reservoir. 

The water in the pipe will then expand through the open end, creating a velocity - V, or V in the 
opposite direction, and relieving the pressure rise. The pressure at the open end thus drops to its normal 
or original value, and this pressure drop travels backward from the open end to the valve. The time re- 
quired for the round trip will b 2L/a. When the return wave or wave of unloading has reached the valve, 
the whole column will be again der normal pressure, but with a velocity - V, discharging water into the 
reservoir. The original actiot jen begins again, but with reversed sign, sinie we now have sudden stoppage 

of flow with velocity - V inste of + V as at first. A wave of rarefaction or of subnormal pressure - y 

now srarts at the valve and makes the round trip, restoring the original velocity + V. At the end of the 
foregoing cycle everything is again in the same condition as when the valve was initially closed, so that 
the cycle begins again and repeats itself continuously until the original kinetic energy is dissipated in 
friction, energy rejected into reservoir, etc. If the valve should be periodically opened and closed in 
synchronism with the cycle, additional pressure waves could be superposed until the pressure is increased 
indefinitely. Fig. 3 shows three diagrams of pressures plotted against time for three different points 
along the length of the pipe line, the closure in this being instantaneous. 


If the valve is only partially closed (but instantaneously) so that the velocity in the region near 
the gate is not completely destroyed but reduced to v, the change in velocity will be V - v instead of V 
as in the previous case; and the rate of change in momentum will now be F * (V - v) Idm/dt), and the rise 
in pressure head will be y = a (V-v) / g. This pressure rise will as before be maintained at constant 
value at the valve until the wave front reaches the open end and returns to the valve. If the valve closes 
only enough to reduce the velocity an infinitesimal amount dv Jf the corresponding pressure rise will be 
dy 5 * (a dVj/g). If, before the wave front completes the round trip, the valve is closed another similar 
amount, a second increment of pressure rise dy 3 is superposed 00 the first. The picture of the situation 
is indicated in Fig. 4. Here A - A is the wave front due to the first increment of closure, and B - B is 
the second wave front due to the second increment of closure. Hence, after any gate movement, the total 
rise in pressure head will be 

r adv , a. (V - V,) 

J g g 


This will be true so long as the first pressure wave has not had time to make the circuit and to add 
(algebraically I a pressure drop to the foregoing values. If therefore the valve is not instantaneously 
closed, but '.is partially closed within a time of al/a. seconds, the total rise of pressure head will be 
y*a(V-Vj)/g, as before; and the amount of total pressure rise is therefore independent of the time 
or manner of closing so long as the closure takes place within t¥e~time 2L/aT If the closure occupies a 
longer time than al/a, the pressure rise will in general be not so great; and its computation becomes 
more complex. (See, for example, N. R. Gibson, "Pressures in Penstocks," Trans. A.S.C.E., Vol. 83, 1920; 

R. S. Quick, "Comparisons and Limitations of Various Water-Hammer Theories," Mech. Eng., Mid-May issue, 
1927, p. 524.) For non-instantaneous closure in a time less than 2L/a the situation may be visualized by 
referring to Fig. 5. This sketch indicates the conditions when sufficient time has elapsed after the be- 
ginning of valve movement for the initial pressure wave to have traveled upstream as far as A - A. The 
effect of the completion of the valve closure has traveled as far as B - B. These two wave fronts are 
traveling upstream at the same rate, the space between them thus remaining of constant length and repre- 
senting a transition region in which the velocity varies from V at A - A to zero at B - B, and in which 
the pressure head is increasing from normal at A - A to normal plus y at B - B. All of the water upstream 
from A - A still has its original velocity V and original pressure, and all the water between the valve 
and B - B has been bronght to rest and is compressed and nnder the same pressure head, y feet above normal. 
As the region A - A, B - B moves uniformly upstream, that part of the column having the velocity V is being 
continuously and uniformly replaced by an equal portion having zero velocity, so that momentum is being 
extinguished at the same rate as in the case of instantaneous closure, and the pressure rise is the same, 
as just deduced. 


Velocity of the Pressure Wave. In the previous derivation we utilized the equation of forces, or^ 
principle" of impulse and reaction. We have not yet found the value of a, the velocity of propagation of 

the pressure wave in the pipe. To find this we shall now apply the equation for volume, or quantity of 

flow After the instantaneous closure of the valve, a region of compressed fluid is created in the end 
of the pipe, extending upstream a distance x at any time t seconds after closure. In this region the 
water is under excess pressure and the pipe walls are expanded as shown to an exaggerated degree in the 
sketch (Fit* 6). If A - A is a fixed cross-section of the pipe, which has at the instant considered been 

just reached by the wave front or boundary of the region of compression, then the volume of water which 

has passed through the section A - A during the time t, with the constant velocity V, is AVt cubic feet. 
Since during this time there has been no outlet, this volume must have been accommodated by compression of 
the water and by the stretching of the pipe walls, so that 


. A 1 d (Volume) x x it j d 2 _ 
1 Volume u 1 
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Here d (Volume) /Volume is the proportional decrease in volume, or the ratio of a small change in volume 
to the volume of the water; that is, the voluminal "strain" due to compression. Now the "bulk modulus" 
k of a fluid or solid is the ratio of unit compressive stress or pressure intensity, to the unit strain, 
or k = p/[d (Vol.) /Vol.]; so that in the equation d (Vol.)/Vol. may be replaced on p/k, k being the bulk 
modolus of water. 


Also the linear extension of the pipe walls depends on the tensile modulus E of the pipe material, 
and E = unit stress/unit strain = S/fd, - dl/d. But for a pipe 2 e S * pd, so that E = pd 2 /2eldj - d) and 
d, - d = pd 2 /2eE. In the equation, d 2 - d 2 may be written (dj - d) (dj + d), and since the difference 
between d, and d is small in comparison with d, we can put dj + d « 2d; hence the last term of the equa- 
tion becomes 

xJL (d 2 - d 2 ) = illl,-dl 2d = ill l! 2 ^ * x (SJ 1 \ = lM_ 

it ‘ 4 4 2eE \ 4 / eE efc 

(putting A for it d 2 /u. Our equation now takes the form: 

AVt = A xJL +_£*£??; or Vt * xp f-L * 

k e i^k eE J 

Since the elastic moduli k and E are customarily stated in pounds per square inch, and p may be taken 
in pounds per square inch, the foregoing equation will be dimensionally correct if all of its quantities 
are expressed in inch units of length. From the last equation we can obtain at once the velocity of the 
pressure wave 



t 



As just pointed out, V is in inches per second and a is consequently in inches per second, 
dividing both sides by 12, the same expression is correct if both a and V are expressed in feet 
But, expressing p in terms of pressure head; p = wy/144, w being in pounds per cubic foot and y 
and from the Joukowsky formula previously derived p = !w/i 44 l (aV/g), so that 


bu ? by 
per second, 
in feet, 



or a 


'J 


4720 

TU 

Ee 


Here a is in feet per second, 
both be in feet or both in inches; 


w is in pounds per cubic feet, g is in feet per second 2 ; e and d may 
and k and E are in pounds per square inch, k usually equals 300,000. 


If, finally, we now apply the principle of energy to the problem, and compare the.original kinetic 
energy of the moving water column and the work done in compressing the water and stretching the pipe walls, 
it may easily be shown that this work is equal to the original energy; so that during the action the kinetic 
energy of flow is merely transformed into elastic energy of strain, and none is "lost" or converted into 
heat. This, of course, assumes the pipe to be perfectly elastic, which is very nearly true for stresses 
less than the elastic limit. It also neglects the friction loss of the flow in the pipe, and neglects the 
variation of velocity across the pipe. Actually the original kinetic energy is finally dissipated in pipe 
friction, loss of velocity head due to non-uniformity of velocity over a section, and velocity head rejected 
into the reservoir; but the dissipation isslow, and the pressure waves may continue for a long period, 
being gradually damped. This damping affects but little, however, the pressure curve for a single cycle or 
several cycles, under usual conditions. 
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Water Hammer 

Assume a pipe with a valve at its outlet discharging from a 
reservoir. If the valve is suddenly closed, a dynamic pressure, 
in addition to the normal static pressure, is created within the 
pipe. .'This dynamic pressure is commonly called water hammer. 

•It is caused by the sudden transformation of kinetic energy to 
pressure energy. " If P is the pressure due to water hammer and 
.1/ is the mass of water in the pipe whose velocity is reduced the 
amount dY in the time interval (It, then 

r-vf mss; 

If dl = 0, the pressure becomes infinite. Instantaneous change 
in velocity is not, however, possible. In the following nomen- 
clature, all units are expressed in feet and seconds except as 
noted: 

b = thickness of pipe walls 
d = inside diameter of pipe 
a = cross-sectional area of pipe, sq ft 
e = modulus of elasticity of pipe walls, lb per sq ft 

e' = modulus of elasticity of water, lb per sq ft 

E = modulus of elasticity of pipe walls, lb per sq in. 

E' = modulus of elasticity of water, lb per sq in. 

g = acceleration of gravity 

AH = head due to water hammer (in excess of static head) 

H — total head producing discharge through valve 
H o = original head in pipe 
L — length of pipeline 
T = time of closing valve 

Fo = mean velocity of water in pipe before closure of valve 
p* = velocity of pressure grave along pipe 

Let li, h, h, . . . , In represent successive infinitesimally 
short sections of pipe beginning at the outlet. The instant 
the valve is closed, the water in section h is brought to rest, its 
kinetic energy is transformed into pressure energy, the water 
is somewhat compressed, and the pipe wall with which it is in 
contact expands slightly. Because of the enlarged cross- 
sectional area of h and the compressed condition of the water 
within it, a greater mass of water is contained within this 
section than before the closure. A small volume of water has 
flowed into 1, after the valve was closed. An instant later, a 
similar procedure takes place in 1* and then in U, causing a 
wave of increased pressure to travel up the pipe. The instant 
this wave reaches the reservoir, the entire pipe is expanded 
and the water within it is compressed. There is now no 
moving mass of water to maintain this high pressure, and the 
pipe begins to contract and the water to expand, with a con- 
sequent return to normal static pressure. This process starts 
at the reservoir and travels as a wave to the lower end. During 
this second period some of the water stored within the pipe 
flows back into the reservoir, but on account of the inertia of 
this moving mass, an amount flows back greater than the excess 
amount stored at the end of the first period, so that the instant 
this second wave reaches the valve, the pressure at that point 
drops, not only to the normal static pressure, but also below it. 
A third period now follows during which a wave of pressure 
less than normal static pressure passes up the pipe to the 1 
reservoir. When it reaches the reservoir the entire pipe is under' 
less than static pressure, buf since all the water is again at 
rest, the pressure in 1„ immediately returns to the normal static 
pressure due to the head of water in the reservoir. This starts 
a fourth period marked by a wave of normal static pressure 
moving down the pipe. When the valve is reached, the pres- 
sure there is normal, and for an instant the conditions through- 
out the pipe are similar to what they were when the valve was 
first closed. The velocity of the water (and the resultant 
water hammer) is now, however, somewhat less than it was at 


the time of closure because of friction and the imperfect 
elasticity of the pipe and the water. 

Instantly another cycle begins similar to the one above 
described, and then another, and so on, each set of waves succes- 
sively diminishing, until the waves finally die out. 

Equation (11-52) shows that, for instantaneous closure of 


valve, the pressure created would be infinite if the water were 
incompressible and the pipe inelastic. It is impossible to close 
a valve instantaneously, but if the valve is completely closed 
before the first pressure wave has time to return to the valve 
as a wave of low pressure, or in other words, if T is less than 
2L/v a , the pressure will continue to increase up to the time of 
complete closure, and the resulting pressure is just the same 
as if the valve had been instantaneously closed. If T is greater 
than 2L/v a , the earlier pressure waves return as waves of low 
pressure and tend to reduce the rise of pressure resulting from 
the final stages of valve closure. 

Hence, if T is equal to or less than 2L/iv, h wifi be the same 
as for instantaneous closure, but if T is greater :han 2L/v w , h 
will be diminished as T increases. 

The kinetic energy contained in the moving .inns of water 
within the pipe is 

T . T , a/ r o 2 , vy- 

HE = — ~ — = waL —y- ( 1 1-o.i) 

Z -{/ 


This energy is used in the work of compressing the water 
and stretching the pipe walls. If the increase in pressure 
head is AH ft, the compression of the water column absorbs 
(w &H)*aL/2e' ft-lb of energy because the average pressure 
force of the water is w AHa/2 and the distance through which 
this force moves is w AHL/e' . 

In a similar manner the work done in stretching the pipe 
walls is (w AH )~LaD/’2eb. Then, relating the energy of the 
water to the work done in stopping it, 

, TV (iv Aliy-aL (w MIY-LaD ni - n 

WaL 1g = 2? h 2ib 

or 

«*■»> 

As this pressure wave moves along the pipe with a velocity v w , 
the force may be equated to the rate of change in momentum 
as follows: 



w AH a = — v^aV » 

(11-50) 


9 


or 

AH = - Fo 

(11-57) 


9 


and 




v " ~ ^w/g{ 1/e' + D/eb) 

(11-58) 


This derivation is approximate in that the stress in the longi- 
tudinal direction is neglected. Yv hen this is included Eq. 
(11-57) becomes 1 


^ w/g( 1/e' - 


w/gt, 1/e' + Dc/eb) 


(11-59) 


where c is related to Poisson’s ratio and the manner in which 
the pipe is held in place. Values of c are somewhat less than 
unity; 2 however, it is on the safe side to use Eq. (11-57), in 
which c is assumed to be unity. The analysis presented here 
is credited to Professor N. Jankovsky of Moscow. 3 The 
analysis may also be carried out by writing the basic dynamic 
and continuity equations in a manner similar to the derivation 
presented for open channels in the first part of this section. 4,5 


1 John. Parmakian, “Water Hammer Analysis,” Prentice-Hall, Inc., 
Englewood Cliffs, N.J., 1955. 

J Ibid. 

3 0. Simin, Water Hammer, with Special Reference to Researches of 
Professor N. Jankovsky, Proc. Am. Water Works Assoc., vol. 24, 1904, 

4 Ibid. 

•George R. Rich, “Hydraulic Transients,” McGraw-Hill Book Company, 
Inc., New York, 1951. 

* N. R. Gibson, Pressures in Penstocks Caused by the Gradual Closing 
of Turbine Gates, Trans. ASCR, vol. 83, 1919. 



pressure w nen i < m/v*. ror ims cu.su an « 

computed directly from Eq. (11-55) or from Eq. (11-57) after 

v w is determined from Eq. (11-58). 

If for example, tv is 4,680 ft per sec and F 0 is 11.75 ft per sec, 
closure of a valve in a time T 1 2L/v u will create an additional 
pressure head of the gate of 

iff = V ^1 = 1,770 ft 
9 

Pressure When T ? 2 L/v w . For slow gate closure the 
computations may be carried out in steps by a method sug- 
gested by Gibson' or by means of a graphical solution. Both 
methods will be described and illustrated by means of a numeri- 
cal example. In the step solution Eq. (11-57) is written as 
follows: 

A//=-AF (11-60) 


where A H is the increase in pressure head resulting from one 
step in the valve closure, which reduces the velocity in the pipe 
by the amount aF. It is convenient to write Eq. (11-60) in 
terms of the ir d condition as follows: 


AH _ tvFo AV 
H o qH o 1 » 


(11-61) 


Relations between head and pipe velocity for various gate 
openings must next be developed. It is assumed that this 
relationship can be expressed as follows: 


Q = aV = C d A 0 V‘2 0 H (11-62) 


where a is the area of the pipe, V is the average velocity in the 
pipe, CdA, is the product of the coefficients of discharge and the 
area of the gate, and H is the head at the gate. Then, for full 

gate opening, 

Q = aF 0 = (Cd-4;. : V2yffo (11-63) 

and, in general, 

V CdA , .177 

Fo (Cd.4 8 )o >#0 

or, for simplicity, 

V f]i 

Fo T '//; 

where t is CdAJ{C d A „) 0 . Values of r vary from 1.0 for a fully 
open gate to zero when the gate is closed. Solutions of slow 
closure problems consist of the step-by-step simultaneous solu- 
tion of Eqs. (11-61) and (11-6.5), as will be illustrated by the 
following example. 

Example' 11-2. Determine the rise in pressure in a penstock 
leading to a power plant if the turbine gates are closed in 2.1 
sec. H o = 165 ft, L = S20 ft, Fo = 11.75 ft per sec, and 
Vw = 4,680 ft per sec. The solution of the problem will be 
carried forward in six time increments of 2.1/6 = 0.35 sec, 
arbitrarily chosen equal to the travel time of the pressure wave 
from the gate to reservoir and back to the gate 


(11-64) 

(11-65) 


2 L 

— = 0.3o sec 
v w 

A sketch of the penstock is shown in Fig. 11-5. The gate is 
closed in such a manner that six uniform decrements of area 


occur during each of the six time intervals. Then r has the 
values shown in column 2 of the accompanying table. In 
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solving the problem it helps to visualize that the head at the 
gate is increasing linearly during the time of uniform closure. 
Then, at t = 0.167 sec, the front of the pressure wave has just 
reached the reservoir and the pressure in the pipe is as shown by 



line 1-2 in Fig. 11-5. At t = 0.35 sec, the pressure is the sum 
of the positive pressures indicated by the line 3-4 and the nega- 
tive pressures indicated by the line 5-6. Thus, at t = 0.35 sec, 
the pressure wave has just returned to the gate G but has zero 
value at that point. Therefore, at time t — 0.35 sec, there is 
only an increase in pressure at G computed as follows, noting 
first that Vi = F 0 — AFi, Hi — H 0 + Affi, and r = 0.833, 


the subscripts indicating values at the end of the first time 
interval. Equation (11-65) becomes 


Fi _ Alj 
Fo T„ 


V/l: = °- s33 V 1 


AH_i 

//o 


( 1 1-6'w ) 


Also, Eq. (11-61) becomes 


Alii _ v w V o A Fi _ 4, (W0 X 11.7,5 AF, AF, 

H 0 gH o Fo 32)2 X 165" F„ ° 

( 1 1 -0 1 a ) 

These two equations are then solved by successive approxima- 
tions. As a first trial assume AF,/F 0 = 0.10. Then, from 
Eq. (ll-61a), 

~ = 10.36 X 0.10 = 1.04 
and Eq. (ll-05«) becomes 


1 - 


AF, 

Fo 


= 0.833 


A Hi 

~Ho 


1 - 0.10 = 0.833 V X 1 + 1-04 
0.90 X 1.19 


Because this is not an equality, further trials must be carried 
out. The final solution for this time interval gave values of 

AF,/F 0 = 0.033 

and AH/Ho — 0.34. Then F,/F 0 = 0.967 and 


Fi = 0.967 X 11.75 = 11.35 ft per sec 

Similarly, Hi/H 0 = 1.34 and Hi = 1.34 X 165 = 221 ft. 
These values are recorded in the table. 

During the second time interval the solution of Eq. (11 -G la) 
is carried out as before, except for the subscripts, because this 
equation deals only with the pressure increment resulting from 
another reduction in velocity. 

AH, AF> 

~ = 10.36 ffv (1 1-6 la) 

no b 0 

However, Eq. (ll-65a) must be revised because there is super- 
imposed upon the positive pressure wave the negative wave 
returning from the reservoir. At T = 0.70 sec this negative 




wave will cause a reduction in velocity equal to the increment 
from the first interval denoted here as Al'i (AlT/T 0 = 0.033) 
and a reduction in pressure A/7 0 (A//i//7 0 = 0.34). These 
values must be included in the discharge relationship at the 
gate as given by Eq. (ll-Goa) as follows: 


V ; 1 1 A V l 

l 0 l' 0 l 0 


AT; 

To 



(ll-65a) 


Inserting numerical values, 

0.9G7 - 0.033 - ^ = 0.G67 V 1.34 - 0.34 + ~ 

Vo ’ H o 

or 

0.934 - ^ = 0.667 \/l + ^ (ll-65a) 

V o ' -W o 

Solution yields AT 2 /T 0 = 0.067 and AH 2 /H 0 = 0.693. Then 
T 2 /Ta = 0.867 and H.JH o = 1.693. Corresponding values of 
T ; and H 2 are shown in the table. Computations were con- 
tinued in this manner, and the results are tabulated. 

Much less tedious than the step computations described 
above is a graphical method' of solving slow-closure problems. 
The equations upon which this method is based are obtained 
by solving the differential equations for dynamic equilibrium 
and continuity. 2 For any two points in a pipeline such as 
R and G in Fig. 11-5, the following equations can be written: 


hati ~ hn tl 

hm, — hat. 


tv T 0 
gH o 


(fed 


ttsij) 


tvTo 
gH o 


(t’Hi, 


ttCij) 


( 11 - 66 ) 

(11-67) 


In these equations h = H/H 0 and v = T/T 0 and ti, U_, and U 
are successive times separated by uniform time intervals. 
Furthermore, hat v t'ci, and hm,, v [ih are coordinates of points on 
a graph such as that shown in Fig. 11-6, and Eq. (11-66) repre- 
sents a straight line having a slope v w V 0 /gHt> connecting these 

two points. The points are designated in Fig. 11-6 by the 
subscripts only. A similar line is represented by Eq. (11-67), 
except that it has a negative slope. These two equations are 
derived from the dynamic relationship resulting from the water 
hammer and thus correspond to Eq. (11-61). The discharge 
relationship at the gate is taken care of in the graphical solution 
by plotting the parabolas represented by Eqs. (11-64) and 
(11-65) for various values of r as shown in Fig. 11-6. The use 
of these equations in making a graphical solution will be illus- 
trated by applying them to the same problem solved in the 
previous example by the step method. 



1 Parmakian, op. ext. 

2 R. W. Angus, Simple Graphical Solution for Pressure Rise in Pipes and 
urnp Discharge Lines, J. Eng. Inst. Canada, February, 1935, p, 72. 


Example 11-3. The graphical solution is conveniently made 
by writing the equations for time intervals equal to half the 
round-trip travel time of the wave, or L / tv = 0.175 see. Then 
the first four required forms of Eqs. (11-66) and (11-67) are: 


ha . o — liR.ns = 10.3 G( t'rr o — vr.us) (a) 

hR.ns — ha. 3a = — 10.36(rR. it'j — t’G. 35 ) (h) 

lie. 35 — Afi.525 = 10.36(l’(7.35 — VR. 5 25 ) (c) 

l!fi.5 2 5 — ha. ui = — 10.36(r’ft.5 2 5 — reus?) (d) 


Once the pattern of the equations is envisioned, they can be 
represented diagrammatically 1 as shown in Fig. 11-7, in which 

the four equations shown above are represented by lines a to d. 
It need only be remembered that equations showing a time, 
sequence from G to R have a positive slope, and those proceeding 
in the opposite direction have a negative slope. 

At zero time V = T 0 and II = //„; therefore h G . 0 and va.o 
both have the value 1.0 and the point G. 0 is plotted at that 
location in Fig. 11-6. Furthermore, at the reservoir, the origi- 
nal conditions exist at time t = 0 and t = 0.175 sec, so that 
points Tf.O and A. 175 are also plotted at v — 1, h = 1. Because 

both ends of the line a, as 
represented by Eq. (a), oc- 
cur at the same point, Eq. 
(a) is also represented by 
the same point as shown in 
the figure. The gate-dis- 
charge relationship also is 
satisfied by this point be- 
cause, with t = 1.0, Eq. 
(11-65) passes through 
h= 1, u = 1 . However, at 
time t = 0.35 sec, the con- 
ditions at the gate are rep- 
resented by Eq. ( b ) as well 
as by Eq. (11-65), with r 
equal to 0.S33. This value 
of r may be read from the 

tion of equaiTons,' Exar p irn-3.‘““‘ table shown , in Example 

1T2. Equation (11-65), 
with r = 0.833, and Eq. ( b ) are plotted, and their intersection 
locates point G.35. Value* of h = H.n/H 0 and v = T. 35 /To 
may be read from this point and compared with the tabulated 
values computed by the step method in Example 11-2. 

Equation (c) is plotted through G.35 with a positive slope, 
and the location of 77.525 is found at h = 1.0 because the head 
at the reservoir must remain constant at the original value H a . 
Point G.70 is then located by extending Eq. (d) with a negative 
slope from point 77.525 to its intersection with Eq. (11-65) 
with r = 0.667. The solution is then continued in the same 
manner until the valve is completely closed, after which, 
because friction has been neglected, the pressure continues to 
fluctuate from a maximum positive to a maximum negative 
value as shown by the loop at the left end of the figure. It mav 

be noted that the graphical solution gives results identical with 
the step solution shown in the previous example. 

The graphical solution is much faster than the step solution, 
which usually involves a considerable number of trials for each 
step. An additional advantage of the graphical solution is 
that it can be readily set up to find the pressure at other loca- 
tions in the pipe. This is illustrated by means of the following 
example. 

Example 11-4. For the same penstock used in Examples 
11-2 and 11-3, find the variations in the pressure at mid-length 
(point c in Fig. 11-5) as well as at the gate during valve closure. 
In this problem Eqs. (11-66) and (11-67) are written between 
G and C and also between C and R. The first few equations 
are given below, and the simpler scheme of representing the 
equations described in the previous example is shown in Fig. 



1.925 


2.275 ■< 


2.625 


2.45 
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1 This scheme is patterned after one shown by Parmakian, op. rit. 
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t 

0.0875 -< 


G 

t 

0.0 
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0.350 

0.525 

0.700 

0.875 



0.2625 



0.4375-s- 


0.6 I 25- 


0.7825 -e- 



0.175 

0.350 

0.525 
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0.875 


Fig. 11-8. Diagrammatic representation of equations, Example 11—t 


11-8. For this example the equations are written for time 
intervals of L/‘2v w = 0.0875 see. 


ho.Q - 

“ he. 0875 ~ 

10.3G(i-’g.o — tfc.osiO 

(a) 

Af. 0S7S 

— Ac. 175 == 

— 10.36(t'C.0875 — t'(7. 17s) 

(&) 

At’. 0375 

— A/2.175 = 

10.36Oc.0S75 — VR. 17.-») 

(c) 

Afl. 175 - 

” Ac. 2625 *** 

— 10.36(V/M75 — VC, 262s) 

W 

h(J, 175 * 

“ Ac. 2625 = 

10.360(7.175 — flC.2625) 

(c) 


he. 2625 — ha.ih = — 10.36(rc.2625 — va.si) (/) 

Several steps in the graphical solution are shown in Fig. 
11-9. The reasoning is the same as that of Example 11-3, 


except that points at the center of the pipe are located by the 
intersection of two of the straight lines. For example, C.4375 
is at the intersection of lines h and i. The location of line h is 
established by its starting point at R . 35 and its negative slope, 
and the location of i is determined by its starting point at 
G.35 and its positive slope. It may be noted that points such 
as G.35, G.70, A. 52.5, and R . 70 have the same coordinates as 
determined in Example 11-3 and shown in Fig. 11-6. 



Although the graphical solution is much quicker and easier 
than the analytical solution previously described, it also 
becomes cumbersome when very complicated conditions are 
encountered. Under such conditions water-hammer problems 
can probably be solved most expeditiously by means of the 
digital computer. 1 ' 2 A convenient method of applying the 
computer to water-hammer problems is the method of char- 
acteristics. 2 - 3 The procedure is carried forward in steps in a 
manner similar to that described in the first portion of this 
section, in which the method characteristics were applied to 
flow with a free surface. 

1 F. D. Ezekiel and IT. M. Puynter, “Computer Representations of Engi- 
neering Systems Involving Fluid Transients, ” Trans. AS ME, vol. 79, 
p. 1840, 1957. 

5 V. L. Streeter and Chintu Lai, Water-hammer Analysis Including Fluid 
Friction, Proc. ASCE, J . Hydraulics Div., May, 1902, p. 79. 

3 C. A. M. Gray, Analysis of Water Hammer by Characteristics, Trans. 
ASCE, vol. 119, p. 1176, 1954. 


69 


T/je. o f yv\zoh<r<.^ / r ca / E-Mtj/Mee-rZ f PrpczeJ Voi, / £ 5" 

/f’S'', />. V3-5~^ 


Water Hammer in a Pumping Main and its Prevention 

By A. M. Binnie., M.A. 5 M.LMech.E *, and D. G. Thackrah, M.A., Ph.D.f 

An earlier investigation has been continued on the protection of a rising water main against the 
consequences of a sudden interruption in the supply. In place of the air bottle previously tried, 
an automatic air-inlet valve was placed at the lower end of the pipe-line. Experiments on a laboratory 
scale showed that the cushion of air thus introduced brought the returning water column to rest 
without sound or shock. Theoretical calculations of maximum pressure and of air drawn in were 
verified. The lift of the pipe-line being small compared with its length, it was found that, over the 
range of velocity usually employed, the maximum pressure set up by the returning column diminished 
as the initial velocity of the interrupted supply was increased. 

With no protective device in use, a series of violent impacts took place owing to the column 
striking and rebounding from the closed valve at the lower end of the pipe-line. An electronic gauge, 
having a very high natural frequency, was used to measure the magnitude and duration of the shock 
pressures. Its indications were generally in accord with an analysis based on the usual theory of 
pressure waves, but the maximum pressures actually recorded somewhat exceeded the theoretical 
estimates, owing to the existence of additional pressure fluctuations caused by reflections from bends, 


sockets, and other discontinuities in the pipe-line, 
slighdy reduced the shock pressures. 

INTRODUCTION 

In a long rising main there is a risk that damage may be caused 
by a sudden stoppage of the pump. The stoppage may be acci- 
dental, due to a failure of the electric supply, or it may be 
designed, as when the pump motor is arranged to be switched 
directly on and off by means of a float in the suction well or the 
reservoir. In either case the pressure in the main may be so 
greatly reduced that collapse may ensue. A more serious danger 
is that the column in the main may break; the upper portion, 
after first coming to rest, will then return and All up the cavity 
with a heavy and possibly disastrous impact. Six methods of 
avoi din g these troubles are available — 

(1) Automatic air- inlet valves may be fitted at the positions 
where the column is expected to break. The valves admit air 
to the main as long as the pressure in the pipe falls below 
atmospheric, and this air serves as a cushion to break the 
force of the returning water. This method is of particular value 
in low-head installations in which the water contains solid 
matter. It was shown by Binnie (1945)+ how the positions 
where the column is most likely to break may be calculated. 
The column will break at the instant of stoppage if at all. 

(2) Non-return valves may be placed in the pipe. The 
merits of this method were explained by Angus (1945). 

v (3) An air-vessel may be attached to the main as investigated 
by Binnie (1945) and by Blair (1945). 

(4) Massive flywheels may be added to the pumps in order 
that the stoppage should be very gradual, but this plan is 
undesirable if the pump motor is started by direct connexion 
to the electric supply through simple control gear. 

(5) A relief valve may be fitted to the main and designed to 
lift when the internal pressure becomes excessive. 

(6) The pump may be by-passed by a short length of pipe. 

All these devices are capable of influencing the surges which 
result from a shut-down, but only methods (3), (4), (5), and (6) 
can be effective in reducing the pressure rise due to a sudden 
start. 

The first part of this paper is concerned with method (1); 

The MS. of this paper was first received at the Institution on 15th 
November 1949. For the Minutes of the meeting in London on 12th 
January 1951, at which this paper was presented, see Proc. I.Mech.E., 
1951, vol. 164, p. 103. 

* Fellow of Trinity College and University Lecturer in Engineering, 
Cambridge. 

f Mechanical Engineer, I.C.I., Ltd., General Chemicals Division, 
lately research scholar. Trinity College, Cambridge. 

$ An alphabetical list of references is given in Appendix II, p. 52. 
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its purpose is to show how an estimate may be made of the 
volume of air sucked in when the pump stops, the time taken for 
the column to come to rest, and the maximum pressure developed 
when the returning column first comes to rest under the influence 
of the air cushion. Of these quantities the first is required in 
order that a sufficiently large release valve may be fitted nearby 
to permit the escape of the air from the main. The second is 
wanted for comparison with the probable lapse of time before 
the pump is restarted by means of the float-operated switch, 
while the third is of value for comparison with the normal 
pumping pressure. The necessary theory has been derived and 
compared with the results of experiments with water on a 
laboratory scale. 

The consequences of omitting a protective device are next 
examined. Under these conditions water hammer of great 
intensity is set up, in contrast to the relatively gentle surges 
mentioned in the preceding paragraph. A fairly complete theory 
of shock waves in pipes has long been known, and numerous 
experiments have been performed on them, but many of these 
tests leave much to be desired because the mechanical pressure 
gauges employed were incapable of following accurately the very 
high rates of pressure change which occur. As O’Brien and 
Hickox (1937) remarked, measurements of water hammer when 
compared with theory are in good agreement, which might even 
be improved if more reliable instruments were available. The 
development of electronic apparatus has now made satisfactory 
observations of this kind possible, and tests were carried out on 
the same laboratory pipe-line with the aid of a capacity pressure 
gauge having a natural frequency of over 30,000 cycles per sec. 

Finally, with the same apparatus, trials were made of the 
performance of a relief valve in reducing the shock pressure 
(method (5)). Both theory and experiment revealed that when 
exposed to violent pressure rises this device is relatively 
ineffective. 

THE USE OF AN AIR INLET VALVE 

Theoretical Calculations. The profile of the main having 
uniform cross-section A sq. ft. and static lift H fee t is shown 
diagrammatically in Fig. 1. An air inlet valve is fitted at B, 
where the column is expected to break. After rime t seconds 
from the instant of stoppage, the length x feet of main from 
B to C is filled with air of volume V cu. ft. at absolute pressure p 
feet of water. The length L feet of the main from B to the outlet 
D is taken to be large compared with the maximum value of x.' 
It is not to be expected that the shape of the air bubble will be 
precisely as indicated in Fig. 1. In a large pipe it will tend to be 



elongated, and it will ascend the main faster than the water if 
the slope is sufficient. However, in mains of small gradient the 
upward movement will be slight, while if the main is relatively 
steeply inclined the water will come to rest quickly before the 
bubble can travel far. Thus in an approximate calculation it seems 
sufficiently accurate to assume that the bubble remains close to 
B throughout the motion, and that x is given by V/A. How long 
the pump takes to come to rest, and its frictional resistance (when 



stopping) to the passage of water through it, are matters on which 
information is usually lacking. It will, therefore, be assumed 
that the water in the main below B is instantly brought to rest 
at the moment of stoppage, and remains stationary during the 
motion of the portion CD, being prevented from falling back 
by the non-return valve which is almost invariably fitted dose 
to the pump. This assumption will lead to an over-estimate of 
the volume of air drawn in. 

On these suppositions the equation of motion for the broken 
column is 

H +P-p + kv*=-^ d £ ....(!) 


in which g ft. per sec. per sec. is the acceleration due to gravity, 
v ft. per sec. is the mean velocity of the column, regarded as 
positive in the direction BD, and P feet of water is the atmo- 
spheric pressure, assumed to be the same at B and at the reser- 
voir ; for simplicity the friction loss of head in the main is taken 
to be kv 2 feet of water, where k is a constant. 

The calculations can be pushed to a quick conclusion if the 
valve at B is so large and light that during the process of air 
admission, p may be regarded as a constant, say p 0 , the magnitude 
of which is controlled by a weight or spring; again this hypothesis 
will result in an over-estimate of V. Equation (1) is then easily 
integrated and shows that the total volume V x of air drawn in, 
measured at pressure p Q , is 


Vi = Ax x 


AL 
2 gk 



1 + 


fob 2 \ 
H+P-p 0 ) ‘ 


( 2 ) 


where v Q is the initial velocity, and x x is the value of x where the 
column comes to rest. To find the corresponding time t x it is 
necessary to replace vdvjdx in equation (1) by dvjdt. The 
integration is again simple and yields 

~ j{*(*+P-*o)}-» tan-* • (3) 

After momentarily coming to rest in its extreme position the 
column falls back under the action of gravity. It will be supposed 
that the valve at B instantly closes, so that an air bubble of 
volume V ls calculated from equation (2), is trapped at pressure 
p 0 in the neighbourhood of B. This bubble acts as a spring, and 
the column wall oscillate to and fro until brought finally to rest 
by friction. The maximum air pressure p 2 developed in the 
course of these oscillations is required, and this will occur when 
the column reaches its lowest position for the first time. On the 
assumption that the pressure-volume changes of the air are 
isothermal, p 2 can be determined by an adaptation of the method 
employed by Binnie (1943) for calculating the pressure changes 
in closed surge tanks. Equation (1) holds good, but dvjdt must 
be substituted for vdvjdx , and the sign of k reversed because 


the column is moving in the opposite direction. Thus, by means 
of the equation of continuity 


Av = 


dV 

dt 


(4) 


equation (1) transforms to 


k tdV\ 2 

H+P ~A^Ji) 


L dW 

+ Ag dt 2 ' • 


(5) 


At the start of the return motion the product of air pressure and 
volume is p$Vy, hence the left-hand side of equation (5) is equal 
to PoVijV. As shown in Appendix I the resulting equation can 
be treated to yield V 2 , the air volume when the column is in its 
lowest position. The expression for V 2 is 


EK-rVJ-EH-rVz) = -{exp(-rH 2 )-exp(-rF 1 )} (6) 


where 


r = 


2kg 

AL 


and q 


P*Vi 

H+P ‘ • 


(7) 


This equation for V 2 is readily solved by trial, the calculations 
being greatly facilitated by the tables of the exponential function 
exp(x) and of the exponential integral Ei(x) published by the 
Federal Works Agency (1939 and 1940). Finally, p 2 is found from 
the relation p 2 = p Q VijV 2 . If required, the corresponding time 
can be determined by the method explained in the paper already 
referred to. 


Experimental Apparatus, A diagrammatic elevation of the 
apparatus used to check the theory is shown in Fig. 2. The 



Fig. 2. Arrangement of Air Valve and Air Measurement 
Apparatus 


uniform pipe-line JL, which had a nominal internal diameter of 
2 inches, was provided with water from an elevated tank or, for 
velocities exceeding 6 ft. per sec., from a centrifugal pump. The 
supply passed through the orifice plate G and the control valve 
H to the disk valve J, which could be suddenly shut by means 
of a lever in order to simulate the abrupt stoppage of a pump. 
The air bottle K protected the supply pipe at the instant of 
closure. The length of the pipe-line JL was 164J- feet. The first 
46 feet was straight and almost horizontal; this was followed by 
a 180 deg. bend, a horizontal length of 85 feet, a 90 deg. bend 
and a further horizontal length of 21 feet, after which the pipe- 
line rose vertically to the reservoir M. 1 This tank, which was 
37 X 18 inches in plan and 24 inches deep, was fitted with a 
long overflow weir, so that the level in it remained nearly constant 
whatever the velocity in the pipe-line might be. The difference 
in height between the sill of the weir and the centre of the valve 
J was°ll-7 feet. The orifice-plate G was calibrated in position 
by leading the overflow from the weir to a measuring tank. The 
sluice valve W was inserted to isolate the reservoir when neces- 
sary, and a drain valve was attached to the horizontal part of the 
pipe-line. 

To admit air when the pressure in the pipe-line fell below a 
predetermined value, the mushroom valve N, which was 2 - 5 
inches in diameter and weighed only 0'42 lb., was fitted close 
to J. A stop limited its travel to 1 inch ; in this extreme position, 
the area of the opening was two and a half times the cross- 
section of the pipe-line. In order that the entry of air through 
the tube R should be obstructed to the least possible extent, the 
valve guide consisted merely of two sets of four light horizontal 



rods attached radially to the valve stem, their outer ends sliding 
on the inner surface of the tube. The valve was controlled by 
the light spring P, which was so long that the force it exerted 
was virtually independent of the valve opening. Two alternative 
springs were used, having tensions up to 0-8 lb. and 16 lb. 
respectively, and they could be adjusted by means of the rod 
Q, supported on a fixed horizontal bar. Thus, the valve was 
so designed as to conform as closely as possible to the assumptions 
which were found necessary in the derivation of the theory. 

To remove and measure the air admitted to the pipe-line, the 
valve S was opened and thereupon the water in the reservoir 
vigorously expelled the air into the graduated glass bell jar T, 
which was fixed in the shallow water tank U. The process was 
assisted by a slight fall in the pipe-line from J onwards, which 
also ensured that no air travelled up the pipe-line to escape at 
the reservoir. Before an observation was made, most of the air 
in the jar was sucked out by a water ejector connected to the 
tube V. The scale was calibrated by inverting the jar and pouring 
in known volumes of water. 

The pressure variations in the pipe-line, being slow, were 
satisfactorily measured by a calibrated steam-engine indicator 
(not shown in Fig. 2), which was attached close to the valve J. 
Its drum was rotated uniformly by an electric motor. From time 
to time the frictional loss of head was ascertained from readings 
on an open stand-pipe carried up from the lower end of the pipe- 
line to a position near the reservoir. To this pipe a mercury- 
water gauge was connected when the velocity, and hence the 
accompanying loss, was large. 

Experimental Results. A typical indicator diagram is shown 
in Fig. 3. When the disk valve was closed, the air valve instantly 
opened and the pressure fell abruptly to the suction value. This 
part of the diagram was obscured by instrumental vibrations 
but, being quickly damped out, they were of no importance. 
Thereafter the suction pressure was sensibly constant, as the 


theory requires; when the air valve had shut, the pressure 
performed damped oscillations with the static head as mean. 

Experimental and theoretical values of the volumes (at atmo- 
spheric pressure) of admitted air are given in Fig. 4 for spring 
tensions of 0-8, 8, and 16 lb., which correspond to pressures, 
required to open the valve, of 0-2, 3-5, and 7-3 feet of water. It 
was found that the latter pressures never differed from the 
suction pressures by more than 0-7 foot of water, an amount 
which itself is not far from the limit of accuracy of the diagrams. 
Equation (2) yields the theoretical volume at the suction pres- 
sure, and for the purpose of determining it at atmospheric 
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Fig. 3. Typical Indicator Diagram 
Air valve loaded to lift at pressure 7-3 feet of water. t>o = 6 ft. per sec. 

pressure the passage of air into the pipe-line has been taken as 
isothermal. The agreement between the observed and calculated 
volumes is good; even at the highest velocity in the pipe-line 
the discrepancies are slight. 

The maximum pressures developed at the first return of the 
water column are displayed in Fig. 5, to which the observed 
pumping head has been added for comparison. Again the results 
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Pumping head. 
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are satisfactory, and the important point emerges, both from But the acceleration is given by equation (1), hence 
experiment and theory, that the rebound pressure decreases as L' l L'\ 

v 0 is raised. If employed on a long pipe-line with a small static p' = j-(H+ P)— H'+ (l — TjPo • • (12) 

lift, the air valve is a device the efficacy of which actually _ ' 

improves when the conditions of working are made more Thus, p’ is independent of t and of k , and for a pipe-line of 
onerous. At high rates of flow the rebound pressure fell below constant gradient it varies linearly from p 0 at the lower end to 
the pumping pressure, and if this pipe-line were always operated P at the upper. In a similar way it can be proved that equation 
at 9 f t. per sec. there would be no point in reducing the spring ^(12) is true also for the return motion. 

tensionuelow 16 lb. Indeed, it would be undesirable to do so. At the instant when the cavity is filled up and the lower end 
for no more amfean is necessary should be admitted to the of the column suddenly brought to rest, a shock wave is pro- 

pipe-line because of the difficulty of removing it. duced which travels up the pipe-line, leaving each part in succes- 

The compression of the admitted air by the returning column sion at rest and abruptly increasing its pressure by a constant 
being taken to follow the relation pV n = constant, n was found amount, say p 2 feet of water. The wave reaches the open end 
from the observations of pressure and volume to lie in the range at the reservoir after a lapse of time Lja seconds, a ft. per 
1-1 1—1 -13 j hence the original assumption of isothermal com- sec. being the steady velocity of the wave, and there it is 
pression was not seriously in error. A few of the measured reflected as a wave of tension p 2 , relieving each section in turn 
suction times were compared with those calculated from of pressure and giving it velocity v 2 . After a further interval of 
couation (3) ; the latter were consistently smaller by amounts time Lja the wave arrives at the lower end of the pipe, and at 

varying from 10 to 23 per cent. this moment the entire column possesses velocity v 2 . At the 

The pipe-line was made up of new steam tubing, and the disk valve, which nay be regarded as forming a closed end to 
friction loss of head h increased perceptibly as the experiments the pipe-line, th ension wave is reflected unaltered in sign, 
described in this paper proceeded. After allowance had been and, if (H+P) e eds p 2 , it is able to travel up the pipe-line 

made for the bends and fittings, the coefficient A in the formula in this form. B_ in the present circumstances (H+P) is 

h — \Lv 2 [(2gD) rose, over the range 5-9 ft. per sec., from 0-021 smaller than p 2 , and the water column, being incapable of 

after 29 days to 0-024 after 911 days. withstanding tensile stress, bounces clear of the disk, again 

leaving a cavity behind it. Under the ideal conditions which are 


THE UNPROTECTED PIPE-LINE 

After the air valve and the air measuring apparatus had been 
removed, it was found that, when the disk valve was suddenly 
closed, a quiet interval of the order of a second ensued, followed 
by a succession of loud shocks. At the instant of closure a vacuous 
space was evidently formed near the va've while the water 
column in the pipe-line continued onwards. The water returned 
under gravity, refilling the space and striking the valve with great 
violence ; it then rebounded, again leaving a space, and the cycle 
of events was repeated until the water was brought to rest owing 
to the loss of its energy by friction and the transfer of vibrations 
to the surroundings. A detailed account of these effects will now 
be given. 


assumed in this account, the process may be summarized by 
stating that the column, after an interval 2 Lja, rebounds with 
unimpaired velocity. The whole sequence of events described 
above in this section is then repeated with the reduced velocity 
» 2 - 

It is shown in text-books that the relation between p 2 and v 2 is 



where the apparent bulk modulus K' of the water in the elastic 
pipe-line is 


K' = 


K 

DK (5 -4a) 
l + SE 4 


(14) 


Theoretical Calculations. Equation (1) defines the motion of 
the water column up the pipe-line until it comes to rest for the 
first time. An experiment to be described showed that during this 
process the absolute pressure p in the cavity near the disk valve 
possessed a constant value p 0 , therefore Ef and can be deter- 
mined from equations (2) and (3). Usually p 0 is so small that it 
can be neglected in comparison with (H+P). For the return 
motion, equation (1) holds good when the sign of k is reversed 
and p is taken to be the same constant p 0 ; and the velocity v 2 , 
with which the column strikes the valve, is required. The 
equation analogous to (2) is 

and with the aid of equation (2) the result is 

*2 = -Vo(l + H+ y p~-f) • • (9) 

the first minus sign arising from the convention that v is 
measured positively in the direction of pumping. The corre- 
sponding time is found by integration to be given by the 
logarithmic expression 

- 1+ ^2 

~{k(H+P~p 0 )}-i log, - • • (10) 


-{k(H+P-p 0 )}~i log. 


{H+P-pJi 

Next the variation in the absolute pressure, say, p' feet of 
water, must be determined at an intermediate section X (Fig. 1), 
which is H' feet above B and distant L' feet from it, measured 
along the pipe-line. If L' is large compared with the length of 
the cavitation bubble, then, during the first upward movement 
of the column, the equation of motion of the portion BX is 

L'v dv 

H'+p'-p 0 +kv*= -w . . (11) 


if the pipe-line is subjected to longitudinal as well as to hoop 
stress, or is 



if longitudinal strain is prevented. Here p lb. per cu. ft. and 
K lb. per sq. ft. denote the density and the true bulk modulus of 
water, E lb. per sq. ft. and a Young’s modulus and Poisson’s 
ratio for the material of the pipe-line, and D inches and 5 inches 
the diameter and wall thickness of the pipe-line. It is also 
well-known that the velocity of the shock wave is given by 



which, when combined with equation (13), yields 



(17) 


Equation (13) can be obtained from considerations of energy 
or of momentum, and the tacit supposition is made that no 
energy is lost within the system, or finds its way out in the form 
of vibrations. Among possible sources of internal loss are the 
imperfect elasticity of the water and the pipe-line wall, and also 
the friction associated with the longitudinal movement, due 
to compression of each cross-section of the water. In a long pipe- 
line this last effect can be considerable ; for example, the reduc- 
tion in length of a column of water filling a rigid pipe 100 miles 
long when stressed to 100 lb. per sq. in. is about 165 feet. We 
must also expect the wave to suffer partial reflection at any points 
such as bends, branches, and open valves, where the uniformity 
of the pipe-line is interrupted. 

It is now possible to follow the pressure changes due to the 
shock wave at the intermediate section X, and it is convenient 
for purposes of explanation to consider also another section Y 
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(lying between X and the reservoir), the quantities associated 
vdth this being distinguished by double dashes. The pressure 
at X remains at the value -given by equation (12) until the wave 
arrives after a lapse of time L'/a from the moment of impact, 
and it is then abruptly increased by an amount p 2 . After an 
additional passage of time (JJ-L^ja the wave reaches Y. 
Immediately ahead of it the pressure is still given by equation 
(12) with L ’ and H’ substituted for L" and H', and the pressure 
difference across the wave remains p z ; between X and Y the 
column is at rest. These effects do not completely reach X until 
a further time interval (Z/-Z/)/a, and at that instant the 
pressure p' there is given by their sum 

^-{H+P)-H'+[l- I ^jp^+p2+{H'-H') . (18) 

When Y is taken to be at the open end of the pipe-line, the 
expression (18) shows that 

p' = H+P-H'+pz .... (19) 
and at the moment this pressure is attained, the return wave 
passes X, leaving behind it only the static pressure ( H+P—H '). 
The pressure change at X, which occurs between the transits of 
the upward and downward waves, is linear if the gradient of the 
pipe-line is uniform. When the wave reaches the valve, it is 
reflected as a tension wave into a region where the pressure is 
(H+P), and the column breaks in circumstances akin to the 
separation of the two parts of the pressure bar used by 
Hopkinson (1914) in his researches on the pressures caused by 
the detonation of high explosives and the impact of bullets. In the 
ideal case in which the wave has a vertical front, the nose of the 
wave cannot penetrate into the column before fracture occurs. 
Thus, the sudden reduction in pressure which takes place at 
the foot of the rebounding column is from ( H+P ) to p a , and 
this is propagated as a tension wave of magnitude (H+P—po)- 
The experimental pipe-line was almost horizontal for most of 
its length, and in this portion the wave will not reduce the 
pressure below zero, but in the part close to the reservoir, where 
the pipe-line rose steeply, the column may be expected to 
disintegrate. At a point near the valve and continuously in 
contact with the column the pressure changes are relatively 
simple ; they are indicated in Fig. 6. 


velocity of propagation a'. On the assumption that conditions 
over each cross-section are uniform, text-book writers, for 
example O’Brien and Hickox (1937), showed that at 1 the 
transmitted wave is 


2 M „ 


M+ 1 • • ' • 

. . . (20) 

and the reflected wave is 


Af- 1 „ 


M+l' ' ’ * 

. . . (21) 

in which 


A a' 


M = . . . 

A a 

. . . (22) 


Where the portion 1-2 is a fitting, such as a socket, giving this 
length greater rigidity than the rest of the pipe-line, M may be 
expected to exceed unity, and in these circumstances the sign 
of the reflected wave is the same as that of the transmitted and 
the original waves. If, for convenience, the transmitted wave is 
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Fig. 7. Change of Cross-section and Rigidity of the Pipe-line 

denoted by (l + a)P, it follows from expressions (20) and (21) 
that the reflected wave is aP. The transmitted wave continues 
its course until it is partially reflected at 2, where the above 
equations hold good after M has been replaced by its reciprocal. 
Thus, the wave which passes onward along the pipe-line is 

I WVT (1 + “ )p = (1 - a ’ )P ■ • ' ' (23) 

A wave 

!$2=-J(l+«)P--«Cl + «)P - ■ ■ (24) 

is reflected back to i, where the consequences of the discontinuity 
are that a wave 



Fig. 6. Theoretical Pressure Changes in the Column Close 
to the Valve 

The final theoretical discussion, which is required, concerns 
the partial reflection which the shock wave undergoes at a 
section where the uniformity of the pipe-line is interrupted. 
Consider a wave of pressure P travelling with velocity a in the 
pipe-line, the cross-sectional area of which is A , and encounter- 
ing a length from 1 to 2 (Fig. 7), where the area is A' and the 


-i7MTI aC1+a)P== “ a(1 " a2)P • • (25) 

returns towards the origin of the original disturbance, and at 
the same time a wave 

~r/HTi a(1+a)P = a2(1+a)P • • (26) 

is reflected towards 2. In the absence of losses the passage of 
waves to and fro in 1-2 could continue indefinitely, but we see 
that, if a. is so small that its square can be neglected in com- 
parison with unity, the entire process may be epitomized in. the 
statement that a single pressure wave P is transmitted, while two 
reflected waves aP return towards the origin, the first being of 
pressure and the second of tension. 

Experimental Apparatus and Results. The pipe-line had been 
assembled with as few joints as possible, and any that showed 
weakness were welded up to withstand the high pressures now 
developed. Near its lower end the pipe-line passed close to two 
roof stanchions, and it was fixed to these with connexions of 
great rigidity in order to reduce the escape of energy in the form 
of vibration. Under the intense water hammer fine disk valve 
at first gave much trouble, and many trials, including its replace- 
ment by a heavy plug cock, were required before its performance 
became reliable. Ultimately the original valve was found satis- 
factory when its bronze disk had been exchanged for one of 
stainless steel, which passed water-tightness tests at 650 lb. per 
sq. in. without distortion. 

(7) Tests on the Cavitation Bubble. The pressure changes 
which occurred in the pipe-line near the disk valve were 
relatively small until the first impact occurred, therefore tire 
engine indicator, previously mentioned, was employed in the 
same manner to measure them. The experiments required care 
because it was obviously of vital importance that the indicator 
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cock should be shut before the shock pressure developed. 
However, all these hazardous tests, twenty-one in number, were 
performed without accident. After the initial instrumental vibra- 
tions had died down the cavitation pressure was found to be 
sensibly constant, and with values of v 0 = 3, 4, . . .9 ft. per 
sec. it varied inconsistently from 1-2 to 1-9 feet of water absolute. 
Eut these figures are not very accurate for, although the gauge 
pressure could be estimated to 1 per cent, the equivalent 
absolute pressure, being very low, was correct only to about 
25 per cent. Nevertheless, the mean result, 1-6 feet, may be 
taken as being substantially true, and this appreciably exceeded 
the vapour pressure of 0-5 foot of water at the temperature of 
the tests. The difference must be attributed to the liberation of 
dissolved air from the water. 

(2) Measurements of Shock Pressure. To determine the shock 
pressures developed in the pipe-line, an electronic gauge of the 
capacity type was employed, its indications being displayed 
on a cathode-ray tube. The pressure head or “pick-up”, which 
was described by Wolfe, Martin, and Judge (1946), was screwed 
into a boss welded to the pipe-line near the disk valve ; it is 
shown sectioned in Fig. 8. Essentially the head was a parallel- 



Fig. 8. Pressure Head 


plate air condenser, the upper plate being rigid and the lower 
plate consisting of a diaphragm exposed to the pressure in 
the pipe-line; thus the application of pressure deflected the 
diaphragm and increased the capacity of the condenser. A mica 
disk was introduced between the plates, and this, having a high 
dielectric constant, raised the sensitivity by reducing the virtual 
air gap and at the same time excluded the danger of the plates 
touching. The diaphragm A of thickness 0-0546 inch and effec- 
tive area 0-1901 sq. in. rested on an annular seating, and the 
space below it was connected to the pipe-line by several small 
perforations, rather than by one large hole in which resonance 
might more easily be maintained. They were packed with 
vaseline which, while allowing transmission of pressure, removed 
the possibility- of water penetrating into the interior of the head. 
On top of the diaphragm the spacing ring B was inserted, and 
above this a sub-assembly, consisting of the inner electrode C 
(or fixed plate of the condenser) to which the mica disk D was 


stuck, the centralizing ring E, and the three Corundite washers 
F; all were fixed in position by the retaining sleeve G. The 
setting of the air gap H was a delicate operation. Before fitting, 
the sub-assembly and the spacing ring were held together and 
lapped until their faces were level and truly perpendicular to 
their common axis. The sub-assembly was further lapped until 
0-0015 inch had been removed from the face of the inner 
electrode, and in the gap thus formed a mica disk 0-00075 inch 
thick was used. To ensure stability of calibration it was found 
important that the tightening down should be done evenly with 
a box spanner and should be. absolutely solid so that the 
diaphragm returned after the application of pressure always to 
the same position. This process necessarily involved distortion, 
and to obtain a linear calibration a trial and error procedure 
was required. The gauge was calibrated in position with_the aid 
of a hydraulic accumulator connected to the pipe-line through 
the drain. If plotted with capacity vertically and pressure hori- 
zontally, the calibration curve at a first trial could be convex 
upwards, showing that at high pressures th mica was being 
compressed by the diaphragm and that m e metal must be 
removed from the inner electrode. Alter ively, the curve 
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Fig. 9. Typical Pressure Gauge Calibration 

was convex downwards, indicating too large an air gap, and 
this could be put right by lapping the spacing ring. Thus, after 
a few attempts the typical straight calibration shown in Fig. 9 
was obtained. It was found necessary to recalibrate at intervals 
of a week, because the sensitivity slowly increased, probably 
owing to the components bedding down and reducing the air 
gap. 

With the aid of the analysis due to Lamb (1921) it was 
estimated that the natural frequency of the diaphragm, if damped 
round its periphery, would be 88,000 cycles per sec. in vacuo 
and 75,000 cycles per sec. when water was present on one side. 
The nature of the edge fixing is in doubt, but even if the disk 
were simply supported, the formula given by Prescott (1924) 
shows that the gravest frequency in vacuo would fall only 
by one half. It therefore appears probable that the gravest 
frequency of the diaphragm as used was not below about 
37,500 cycles per sec. Now the highest harmonic component 
of the pressure, which is to be measured, should have a frequency 
considerably less than that of the diaphragm if a response 
independent of frequency is to be obtained; Wolfe, Martin, and 
Judge (1946) and Wolfe (1949) stated that a ratio 0-3 between 
the two’ is suitable for design purposes. It appears, therefore, 
that the pressure head should accurately follow pressures 
fluctuating at 12,500 cycles per sec. 

A block diagram of the electrical apparatus is shown in 
Fig. 10. It was necessary to change the capacity variations in the 
pressure head into voltage variations suitable after amplification 
for display by the cathode-ray tube, and this was achieved by 
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means of an oscillator and detector unit employing A.C. 
polarization and amplitude modulation. A circuit diagram of the 
unit appears in Fig. 11 and a list of components in Table 1. 
The valve V l5 chosen because of its freedom from microphony, 
was connected as an electron coupled oscillator, its screen grid, 
control grid, and cathode being joined to a Colpitts oscillator 
circuit with a frequency of 1 megacycle per sec. The anode, in 
conjunction with its load Rj and coupling condenser C 4 , formed 
a buffer amplifier for the oscillator. The output from the buffer 
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Fig. 10. Block Diagram of Electronic Apparatus 


amplifier was modulated by a tuned circuit, comprising the 
inductance L 2 (thirty-six turns of 22 S.W.G. enamelled copper 
wire on a 1-inch former), the pressure head in series with a 
padder condenser C 16 , the capacity of the coaxal cable linking 
the pressure head to the oscillator and detector unit, and a 
variable capacity. Coarse variation was achieved by Si and the 
condensers C 5 to C 14 , fine variation by C 15 . After modulation the 
wave was fed to the grid of V 2 , selected for its high mutual 
conductance, 3-6 milhamp per volt; and constancy of bias 
potential was achieved by connecting the cathode to the po- 
tentiometer R4R5 across the high-tension supply. The filament 
required 4 volts, obtained by dropping the general 6-3-volt 
filament supply through the series resistor Rio- The anode of V 2 
was supplied from an unearthed 120-volt dry battery, the anode 


load Rg being connected between its negative end and earth, and 
the voltmeter and push-button switch S 2 were used to check the 
battery and give warning of its failure. The output terminal was 
joined to a tapping on the battery; thus the output potential, 
indicated by means of a centre zero micro-ammeter, was backed 
off so that zero output voltage corresponded to zero pressure on 
the pressure head. The shunt R^ served to protect the meter 
while the unit was warming up, and was disconnected when 
the output voltage had fallen to a safe amount. 


Table 1. Oscillator and Detector Unit Components 
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The D.C. amplifier had a gain of about 120. It included two 
valves connected in push-pull and matched in order to minimize 
drift. Nevertheless the drift, as measured over times of the order 
of minutes, was found to be considerable. It was of no importance 
during the experiments on shock pressure which were necessarily 
brief, and care was taken to apply each calibration pressure only 
fer a short time and to obtain the zero immediately before and 
after this operation, which was carried out as quickly as possible. 
During visual examination of the trace on the cathode-ray tube 
the time base was used, triggered by the lever of the disk valve. 
For the production of records the X plates of the tube were tied 
by a switch to the earthy focusing electrodes so that the spot 
could move in the Y direction only, and the screen was photo- 
graphed with a moving film camera. In this work the afterglow, 
which w r as useful for visual observation, proved somewhat dis- 
advantageous since it led to the production of “shadows” on the 
photographs. Measurements were confined to the central 2 inches 



Fig. 11. Oscillator and Detector Unit 
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| IS-INCH LENGTH REMOVED 

b With relief valve B set to lift at 265 feet of water (abs.). 

Fig. 12. Typical Record showing the First Three Impacts 
uo = 9 ft. per sec. 


of the 5-inch screen in order to avoid errors due to distortion. 
Movement of the disk valve was indicated by means of two 
contacts on its lever, which momentarily increased the spot 
brightness when the valve started to shut and again when 
closure was complete. To measure the film speed, which was 
usually about 4 in. per sec., a small 50 cycles per sec. voltage, 
checked by the frequency meter, was applied to the amplifier. 

The initial part of a typical diagram obtained with this 
apparatus is given in Fig. 12a, which shows the two bright spots 
indicating valve movement, the slight fall from the pumping to 
the cavitation pressure, and the succession, at decreasing time 
intervals, of impacts of diminishing intensity. Fig. 13a, which 
shows the first impact in detail, is an ordinary photographic 
enlargement, but for examination and measurement the films 
were placed in a Hilger projector having a linear magnification 
of 10. The start of the pressure rise was slightly rounded off, 
owing probably to the nose of the column not occupying the 
entire cross-section of the pipe • at the top the trace turned to the 
horizontal direction without overshooting, therefore the gauge 
was considered satisfactory, and it was not thought necessary to 
measure its natural period. The summit of the diagram was 


disturbed by waves, and nothing could be discerned of the 
changes shown in Fig. 6 which theory suggests ; in this pipe-line, 
however, in which the shock pressures enormously exceeded the 
pumping pressures, these changes were in any event small. 

Ideally, the velocity of the shock wave should be obtained by 
measuring the time interval from the start of the pressure rise 
to the start of the pressure fall. But here the latter point was 
obscured, and consequently the ends of the rise and the fall have 
been used. In this way the mean velocity of travel was found to 
be 4,450 ft. per sec., which is the average of twenty-eight 
observations of the first impact. The theoretical values, 4,600 and 
4,620 ft. per sec., obtained from equations (14) and (15) were 
so close together that the way in which the pipe-line was held 
was a matter of no significance. The fact that they were higher 
than the measured value showed that the wall of the pipe-line 
was slightly less rigid than expected. 

At least four records were obtained at each of the values of 
v 0 = 3, 4, . . . 9 ft. per sec. At all velocities the irregularities 
in the traces of the impacts followed the same basic pattern, 
and they were evidently due to the peculiarities of the pipe-line. 
In an attempt to account for the existence of these waves, the 



6 With relief valve B set to lift at 265 feet of water (abs.). 
Fig. 13. Enlargement of First Pressure Rise 
v 0 = 9 ft. per sec. 


a Without relief valve. 
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Fig. 14. Tracing of First Pressure Rise without Relief Valve 
t?o = 9 fc. per sec. 


positions of all the various fittings have been added to scale in 
Fig. 14, which shows a tracing made on the Hilger projector. 
Their order is : — 

(a) At 37-8 feet from the disk valve., a small set of 4 inches 
occupying a length of 11 inches. 

(i) At 46-0 feet, the 180 deg. bend made up of three sockets 
and two 90 deg. bends. 

(c) At 53-6 feet, the drain consisting of a tee and 18 inches 
of ^-inch piping closed by a valve ; this being a branch, the 
principal reflected wave should be one of tension. 

(d) At 73-5 feet, a socket. 

(e) At 93-2 feet, the sluice valve. 

(/) At 110-9 feet, a socket. 

( g ) At 130-7 feet, the horizontal 90 deg. bend. 

(h) At 151-6 feet, the vertical 90 deg. bend. 

It will be seen that the agreement between these positions and 
those of the waves is on the whole good, thus the waves can be 
explained as being due to reflection from the interruptions in 
the uniformity of the pipe-line. Their period is surprisingly long, 
perhaps because the fittings were of no great length and hence 
the assumption on which expressions (20) and (21) are based 
was to some extent violated. Moreover, each reflected wave, 
except that due to (a), had to pass on its return through at least 
one fitting, and no account has been taken of the complicated 
higher-order effects caused by the reflection of reflected waves. 
The maximum pressure, seen on the right in Fig. 13a, was due 
to the superposition of the waves set up by (/) and (g). 

The means of the observed shock pressure at the first five 
impacts are shown in Fig. 15a. For this purpose the height of 
the trace at the shoulder on the left was measured because there 
the diagram was undisturbed by reflected waves. For the first 
impacts the results agreed well among themselves, the largest 
departure from the mean being 2 per cent, but this figure 
increased to 11 per cent for the fifth impacts. Loud intermediate 
shocks were heard during the later impacts, suggesting that the 
column was disintegrating as theory predicts, and the decay of 
repeatability may have been caused by the critical nature of 
these effects. For the first and second impacts the theoretical 
shock pressures, obtained from equation (17) and the observed 
value of a, are also given in Fig. 15a. For the first the agreement 
is satisfactory, but for the second the theoretical values were 
too high, probably owing to the loss of energy at the first impact 
which was not allowed for in the calculation. The shock pressure 
rose with increase of v 0 , but w*as not linearly proportional to it 
because of the effects of friction on the moving column. 

The time interval between complete closure of the valve and 
the arrival of the shock wave was found to vary in a repeatable 
manner from 0-65 second at v 0 — 3 ft. per sec. to 1-64 seconds 
at = 9 ft. per sec. These times were about 7 per cent lower 


than the theoretical values ; they were, of course, large compared 
with the time 0-0739 second required for a shock wave to travel 
up and down the pipe-line, and with the time of closing the valve 
which ranged in a random way from 0-050 to 0-068 second. 

Tests were carried out with gradually decreasing values of v 0 
to determine the minimum velocity at which water hammer took 
place. A faint click could still be heard at a velocity of 0-38 ft. 
per sec., but at 0-36 ft. per sec. no sound was audible. Sir 
Geoffrey Taylor has since pointed out that, as the time for the 



D 0 -FT. PER SEC 


Fig. 15. Shock Pressures 
a Vfithout relief valve. 

O O Observed. • • Calculated. 

Pumping head. 

b At first impact with relief valve B. 

(1) Valve lifting at 265 feet of water (abs.). 

(2) Valve lifting at 496 feet of water (abs.). 

O O Observed. • • Calculated. 

Pumping head. 

valve to close was less than that taken for a wave to pass to the 
reservoir and back, a theoretical estimate of the minimum 
velocity can be obtained by means of equation (17). This 
velocity was very low, therefore the friction head was negligible, 
and the limiting pressure change to be considered is from the 
static absolute head 45-6 feet to the cavitation pressure 1-6 feet. 
It then follows from equation (17) that the required velocity is 
0-32 ft. per sec., which is in good agreement with the observa- 
tions. 

THE REDUCTION OF SHOCK. PRESSURE DUE TO A 

SPRING-CONTROLLED RELIEF VALVE 
Theoretical Calculations. The required modifications in the 
theory can be derived on the assumption that the shock pressure 
at the disk and the relief valves, which were placed near together, 
is of constant intensity (p 2 ) while it lasts. The wave as it passes 
up the pipe-line alters the velocity of flow from v 2 not to zero 
but to some value, say u 2 , therefore tue appropriate form of 
equation (13) is 



• • ( 27 ) 
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A calibration curve for the relief valve can be experimentally 
obtained, connecting the pressure difference across the valve 
with the discharge Q under steady condiuons. This curve will 
approximately hold good under shock conditions also, provided 
that the valve possesses sufficient damping and a natural period 
which is high compared with the durauon of the shock pressure. 
Now 

Q = -u z A (28) 

where A is the cross-secdon of the pipe-line, hence there is 
another relation between the unknowns p 2 and u 2 . The magni- 
tude of both can then be determined; a convenient method is to 
plot on the calibradon curve the relation (27), which is a straight 
line when v 2 is specified. Since the shock pressure lasts for a 
time 2 Lja, the total volume of water emitted by the valve is 
2 LQja. The shock wave while returning from the open end 
continues to change the velocity of each secuon by (,v z — w 2 ), 
therefore the column will rebound with velocity (v 2 — 2 uf). 

Experimental Apparatus and Results. The electronic gauge 
already described was again employed, and a relief valve was 
screwed into a boss welded on the pipe-line near the disk valve. 
Two commercial safety valves ^ inch in diameter were tried in 
turn. Valve A, which had been previously used for other 
purposes, possessed a flat seating, and it was set to lift at 100 lb. 
per sq. in. gauge (231 feet of water absolute). Valve B with a 
conical seating was new; it was provided with two altemadve 
springs which lifted at 100 and 200 lb. per sq. in. The valves 
were calibrated in posiuon under steady pressure with the aid of 
the hydraulic accumulator connected to the drain. A graduated 
class cylinder was used to measure the water ejected from the 
orifice in the body of the valve, the weir acdon of the orifice 
being made as definite as possible by blocking its lower part 
with a thin plate fixed with its top edge horizontal. 

A typical result appears in Fig. 12 6, which was obtained with 
valve B lifting at 265 feet of water absolute, the other condidons 
of the test being idendcal with Fig. 12a. An enlargement of the 
first pressure rise is shown in Fig. 136. It will be nodced that 
the series of impacts was not greatly modified by the relief valve. 
The sharpness of each pressure rise was, however, diminished, 
but on the right the superposition of two reflected waves again 
occurred, causing the maximum pressure in the pipe-line. It 
was not easy to decide which part of the diagrams should be 
measured for comparison with the pressure rises calculated from 
equations (27) and (28). The point at the right-hand end of the 
relatively horizontal portion, immediately preceding the well- 
marked reflected waves, was chosen, and the results for the first 
impact with valve B lifting at 265 and 496 feet of water absolute 
are displayed in Fig. 156; those obtained with valve A were 
almost identical. The curve obtained at the higher lifting pressure 
does not extend below i> 0 = 5 ft. per sec. because at lower 
velocities the valve failed to move. The theoretical pressures are 
seen to be somewhat in excess of the experimental, but the 
discrepancy is nowhere serious. Comparison with Fig. 15a 
reveals that the relief valve, when it lifted, caused only a small 
reduction in shock pressure. This device must therefore be 
regarded as almost useless when inserted in a pipe-line to mitigate 
the effect of water hammer. 

As the shock pressures followed in quick succession, it was 
impossible to isolate the discharges of water caused by the 
separate impacts. The total volumes actually ejected were found 
to be considerably less than the theoretical quantities. The 
difference cannot be attributed to the inertia of the relief-valves, 
for both were light, and von Freudenreich (1948) showed that the 
presence of a liquid in contact with a relief valve greatly increases 
its natural frequency. The discrepancy may be due partly to the 
lowness and partly to the variation of the observed shock pres- 
sure, which in the theory was assumed to be constant ; when the 
discharge was calculated from actual pressure diagrams, much 
better agreement was obtained. 
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APPENDIX I 


DERIVATION OF EQUATION (6) 


Equation (5) may be rearranged in the form 
d 2 V 
dt 2 " 


jdV\ 2 s 

s 

\dt! V~ 



where 


Ag 

i = -jrPo v i 


(29) 

(30) 


and q and r are given by equation (7). With the substitution 
dVjdt = yd, equation (29) becomes 



of which the integral is 

y . exp( — rV) — 2sjexp( — rV)-y-+— exp(— rV) + C . (32) 

where C is a constant. The limits of V are V j and V 2 , at both of 
which the left-hand side of equation (32) is zero. Hence, by 
subtraction, equation (6) is obtained. 
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Since Equation (33) contains only the F wave term, this equation 
defines the velocity and head relations in the pipe line in terms 
of F waves. Similarly, Equation (34) defines the velocity and 
head relations in terms of / waves. 

1 See Reference 3. 
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Vibrating Systems 


Nature provides many examples of vibrating systems: trees swaying in 
the wind, atoms in a molecule of water, the motion of the tides, elec- 
tric current in a flash of lightning, and so on. Although the motion is 
vastly different in each of these systems, they have several things in 
common: For one thing, the motion repeats in each regular time inter- 
val, which we call the period of the vibration; second, some type of 
force co antly acts to restore the system toward its point of 
equilibria 


2.1 SIMPLE HARMONIC MOTION 

Consider the very simple vibrating system in Fig. 2.1 consisting of a 
mass m attached to the end of a spring. We assume that the amount of 
stretch in the spring is proportional to the stretching force (which is 
true of most springs if they are not stretched too far), so that in order 
to stretch it a length / a force Kl is required. The symbol K is the spring 
constant or “stiffness” of the spring. 

Since the spring is vertical, the force of gravity on the mass 
stretches the spring by an amount that remains constant. In the 
equilibrium position shown in Fig. 2.1(b), the downward force of 
gravity on the mass (its weight) is just balanced by the upward force 
exerted by the spring: therefore the system is in equilibrium. The 
description of the motion is simplified if we specify the displacement 


t 

Force exerted by spring 


Force of gravity 


1 


1 


i 


of y of the mass from the equilibrium position and the net force Ftha’ 
acts on the mass. The relationship between Fandy is easily shown to be 


- AT. 


( 2 . 1 ) 


' t- v The minus sign in the above equation reminds us that when the 
mass is below ns equilibrium position (y is negative), the net force will 
be upward (F is positive) as shown in Fig. 2.1(a). Thus the force F 
C0U d be calied a restoring force, which always acts in a direction to 
restore m to us equilibrium position. When the restorine force is pro- 
portional to the displacement, as it is in most vibrating systems we 
study, the motion is given a special name: simple harmonic motion. 

or a system in simple harmonic motion, the period is independent of 
the amplitude (size) of the vibration. The frequency /of vibration is 
t e number of oscillations per second, which is obviouslv the 
reciprocal of the period T of one vibration: 


4 Net force ^ 

A -F 

FIG. 2.1 

A simple vibrator consisting of mass 
and spring. In (b) the upward force 
exerted by the spring and the force of 
gravity balance each other, and the net 
force Fon the mass is zero. 


f-\/T. (2.2) 

It is customary to use a unit called the hertz (abbreviated Hz) to 
enote cycles per second. In the case of the vibrating mass-sprin° 
system, the frequency of vibration is given by the formula 

/= h fw (2-3) 


Note that to double the frequency of vibration, the mass m may be 

re uce to one-fourth its original size, or the spring constant Afmay be 
made four times larger. 



Numerical Example 

Suppose a certain spring stretches 0.10 m when loaded with 
2 kg. What is its spring constant? At what frequency will it 
vibrate when loadpd with 2 kg? 0.5 kg? 

Solution: At rest Kl = mg, so 


v mg 

K - ~T 


2(9.8) 

0.10 


= 196 N/m (newtons per meter). 


Its frequency of vibration when m = 2 kg is 



1 


(2)(3.14) 


'196 

2 




6.28 


= 1.6 Hz. 


When loaded with 0.5 kg. the frequency is 

J_ IK = 1 /l9 6 = ' Jm 

2t rJm (2)(3. 14) V 0.5 6.28 


3.2 Hz. 


A strobe photograph of a mass-spring system was shown in Fig. 
1.9. Also shown was a graph of position y as a function of time. A 
graph of speed v versus time can be made by taking the slope of that 
graph at every time t. Graphs of y and v as functions of time are 
shown in Fig. 2.2. Mathematicians refer to curves shaped like these as 
sinusoidal or sine curves. The maximum value of y is called the 
amplitude. 


2.2 ENERGY AND DAMPING 


The formula for the kinetic energy KE of a moving mass was given in 
Section 1.7: 

KE = \mv 2 , (2.4) 


where m is mass and v is speed. Similarly, the potential energy PE of a 
spring, stretched or compressed a distance y from its equilibrium 
length, was given as 

PE =\Kv 2 . (2.5) 

From the graphs in Fig. 2.2, it is clear that v 2 reaches its maximum 
value when v 2 is zero and vice versa. Thus, the total mechanical energy 
is constantly changing from kinetic to potential to kinetic. At times f, 
and t v potential energy is a maximum, and at /, and / 4 , kinetic energy 
is a maximum. 



FIG. 2.2 

Graphs of simple harmonic motion: 
displacement vs. time; (b) speed vs. 
time. Note that speed reaches its 
maximum when displacement is zero 
and vice versa. 




Any real vibrating system tends to lose mechanical energy as a 
result of friction and other loss mechanisms. Unless the energy is 
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>• 



FIG. 2.3 

Displacement of a damped vibrator 
whose amplitude decreases with time. 


renewed in some way, the amplitude of the vibrations will decrease 
with time, as shown in Fig. 2.3. In many vibrating systems, a certain 
fraction (usually small) of the energy is lost during each cycle of vibra- 
tion; the result is a curve that decreases in amplitude in the manner 
shown in Fig. 2.3. The dotted curve, which indicates the change in 
amplitude with time, is called the envelope, or decay curve. A 
vibrating system whose amplitude decreases in this way is said to be 
damped, and the rate of decrease is the damping constant. 


2.3 SIMPLE VIBRATING SYSTEMS 


FIG. 2.4 

A simple pendulum. 



FIG. 2.5 

A piston free to vibrate in a cylinder. 


Besides the mass-spring system alre? y described, the following are 
examples of systems that vibrate in iple harmonic motion: 


1. Pendulum (small angle). A simpu: pendulum, consisting of a mass 
m attached to a string of length / (see Fig. 2.4), vibrates in simple har- 
monic motion, provided that .v <? /. Assuming that the mass of the 
string is much less than m, the frequency of vibration is 



( 2 . 6 ) 


where g is the acceleration due to gravity. Note that the frequency- 
does not depend on the mass. 


2. A “spring” of air. A piston of mass m, free to move in a cylinder 
of area A and length /, vibrates in much the same manner as a mass 
attached to a spring (see Fig. 2.5). The spring constant of the air in the 
cylinder is determined by its compressibility and turns out to be 
K = ypA/l, so the frequency is 


/= 



(2.7) 


where p is the gas pressure, A is the area, m is the mass of the piston, 
and 7 is a constant that is 1.4 for air. 


3. A Helmholtz resonator. Another common type of air vibrator, 
illustrated in Fig. 2.6, is often called a Helmholtz resonator, after 
H. von Helmholtz (1821-1894), who used it to analyze musical 
sounds. The mass of air in the neck now serves as the “piston,” and 
the air in the larger volume V as the “spring.” The frequency of vibra- 
tion is 



( 2 . 8 ) 


where a is the area of the neck, / is its length, V is the volume of the 
resonator, and u is the speed of sound (r'^344 m/s). 


In Fig. 1 .7(a). the speed remains unchanged in time; in Fia. 1.7(b) 

it increases at a steady rate. The acceleration in Fig. 1.7(a) is therefore 
zero; in Fig. 1 .7(b) it is 

35 cm/s 

a - — A . =70 cm/s-. 

O.a s 

Note the units for acceleration; the unit cm appears to the first power, 
but the unit s appears squared. 

An object in free fall in the gravitational field of the earth ex- 




(a) 


■ *■ ' 

* ■ 

W * 

a 


(b) 

FIG. 1-8 

Stroboscopic photographs of two 
objects in free fail: (a) object dropped 
from rest; (b) object thrown upward. 
(Photographs by Christopher 
Chiaverina.) 


FIG. 1.9 

Vibratory motion in which y, v, and a 
al! change with time. 


periences a constant acceleration of 9.S meters/second : . Thus if it 
begins with no initial speed (up or down), at the end of the first second 
it will have a speed of 9.8 m/s; at the end of 2 seconds its speed will be 
19.6 m/s, and so on. 

Note that acceleration does not always increase speed. If an object 
is thrown upward, the acceleration due to gravity acts to slow it down 
(sometimes this is called “deceleration”). Figure 1.8 shows strobo- 
scopic photographs of two objects in free fall. In Fig. 1.8(a), the 
object is dropped from rest; in Fig. 1.8(b), it is thrown upward, slows 
down, and then begins its descent. 

Figure 1.9 shows a stroboscopic photograph of an object with an 
acceleration that changes its direction with time. The object, a mass 
attached to a spring, is executing a type of vibratory motion called 
“simple harmonic” motion, which will be discussed in Chapter 2. The 
camera has been turned during the exposure, so that multiple images 
of the object create a photographic record of position as a function of 
time. Note that position y, speed v, and acceleration a all change with 
time. 

a 




1.5 FORCE, MASS, AND NEWTON'S LAWS 

Force is a quantity with which we are all familiar; it can be described 
as a “push” or a “pull.” Practically all human activity involves 
forces: running, lifting, eating, writing, even standing. 


The Helmholtz resonator can be thought of as having a mass 
m and a spring constant K that are: 

m = Pal and K- , (2.9) 

where p is the density of air. 


Helmholtz resonators can have a variety of shapes and sizes. For 
example, blowing air across an empty pop bottle causes the air in its 
neck to vibrate at a fairly low frequency. Note that the smaller the 
neck diameter a , the lower the frequency of vibration, which may 
seem a little surprising at first glance. 



FIG 2.6 

A Helmholtz resonator. 


2.4 SYSTEMS WITH TWO OR THREE MASSES 

The vibrating systems considered in the preceding section have one 
thing in common; A single coordinate is sufficient to describe their 
motion. In other words, they have one degree of freedom. In this sec- 
tion, we will consider vibrators with two or more degrees of freedom. 





Such systems have more than one mode of vibration, and the different 
modes will generally have different frequencies. 

Consider the system consisting of two masses and three springs 
shown in Fig. 2.7. The system has two “normal” or independent 
modes, as shown in Figs. 2.7(a) and 2.7(b). In one mode, the masses 
move in the same direction; in the other, they move in opposite direc- 
tions. Assuming equal masses and springs with the same stiffness K, 
the frequencies of the two modes are 
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Note that mode (a) has the same frequency as the simple mass-spring 
system shown in Fig. 2. 1 , whereas mode (b) has a frequency that is 1.7 
times that of mode (a). 


— ► — ► ■* — - — ► 

(a) (b) 

Modes (a) and (b) are virtually independent of each other. That is, 
the system can vibrate in mode (a) with minimal excitation of mode 
•(b), and vice versa. If one sets the system into oscillation by giving the 
two masses a push or pull, the resulting motion will nearly always be a 


FIG. 2.7 

Modes of vibration of . two-mass 
vibrator. The mode s wn in (a), in 
which the masses me in the same 
direction, will have t : . _ lower 
frequency. 


combination of modes (a) and (b). There are many recipes for combin- 
ing these two modes in different proportions, and thus many ways in 
which the system can vibrate. 



FIG. 2.8 

Two-mass vibrators using (a) a linear 
air track; (b) an air table; (c) masses 
and springs hung from an overhead 
rod. 


A great deal about the physics of vibration can be learned from 
watching the motion of a two-mass vibrator. Many physics labora- 
tories have linear air tracks or air tables, on which objects move on a 
film of air with negligible friction. These are ideal for studying two- 
mass oscillators. Another convenient arrangement is to hang the 
masses on long cords from the ceiling or an overhead rod. The cords 
must be as long as possible to minimize the tendency of the masses to 
swing like pendulums. These three arrangements are shown in 
Fig. 2.8. 

On the linear air track shown in Fig. 2.8(a), the masses are con- 
strained to move in one direction only. In the systems shown in Figs. 
2.8(b) and 2.8(c), however, the masses can move at right angles to the 
springs as well. Vibrations in this direction are called transverse vibra- 
tions, whereas vibrations in the direction of the springs are called 
longitudinal vibrations. Some systems vibrate only in transverse 
modes, some only in longitudinal. The air column of a musical wind 
instrument, for example, vibrates longitudinally, whereas the mem- 
brane of a drum vibrates transversely. A violin string normally 
vibrates transversely, althoughJiongitudinal vibrations (which sound 
like squeeks or squeals) are occasionally excited by the bowing of un- 
skilled players. 
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In addition to their two modes of longitudinal vibration, the two- 
mass systems shown in Figs. 2.8(b) and 2.8(c) have two modes of 
transverse vibration, which are shown in Fig. 2.9. In the mode of 
lower frequency, the masses move in the same direction; in the mode 
of higher frequency, they move in opposite directions; this is similar to 
the longitudinal modes shown in Fig. 2.7. 


Adding a third mass to the systems of Fig. 2.7 adds additional 
modes of vibration. In the case of the linear vibrator in Fig. 
2.7(a), which vibrates only longitudinally, a third mode of 
longitudinal vibration appears. The three independent modes 




FIG. 2.9 

Modes of transverse vibration of a 
two-mass system, (a) In the mode of 
lower frequency, masses move in the 
same direction; (b) in the mode of 
higher frequency, masses move in 
opposite directions. 


of vibration are those shown in Fig. 2.9. The systems in Figs. 
2.7(b) and 2.7(c) can vibrate transversely as well; in addition 
to the three modes of longitudinal vibration in Fig. 2.10, thev 
will have the three independent modes of transverse vibration 
shown in Fig. 2.11. 


(a) 


(b) ' 


(c) 

FIG. 2.10 

Independent modes of longitudinal vibrations of a three-mass 
vibrator. 



FIG. 2.11 

Independent modes of transverse vibration of a three-mass oscillator. 


The independent modes shown in Figs. 2.7, 2.9, 2.10 and 2. 1 1 are 
often called the normal modes of the vibrating systems. Getting the 
system to vibrate in a single normal mode requires special care°It is 
perhaps best done by driving the system at the frequency of the desired 
mode (this phenomenon, called resonance, will be discussed in 
Chapter 4). Carefully displacing the masses by the proper amounts 
and releasing them will also cause the system to vibrate in a single 
mode. 
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This is a minimum where dtjdt = 0. Performing the differentiation and setting 
to zero gives 



continuous piston motion with increasing acceleration; the point of shock 
formation moves closer to the origin and in the limit occurs there. 

Let the undisturbed region be denoted by subscript 1 and the region 
behind the shock by subscript 2 so that the notation is consistent with that 
of the shock conditions. With [u ] = -f- X, (7.33) gives 


X 

Cl 


- 2 
y + 1 


( 


M ln 



(8.37) 


This gives the shock Mach number M ]n in terms of Xjc u and all quantities 
of interest {P 2 IPi, T 2 IT U etc.) can then be calculated. 

The rightward-traveling shock shown in Fig. 8.23 is called a C + 
shock because it travels at a speed intermediate between that of the up- 
stream and downstream C + waves, as shown in Chap. 7 (in the limit of 
vanishing strength, the shock path becomes of course a characteristic). 
The motion shown in Fig. 8.23 is divided into two uniform regions by 
the shock path. The flow is then piecewise homentropic. The motion is 
also self-similar, depending only on xjc 0 t, like impulsive withdrawal. 

Each of the motions in this section and in the previous one was 
discussed in terms of C + waves. The results for simple C~ waves (e.g., 
with the piston at the right-hand end of the tube in Fig. 8.17 and ,v still 
positive to the right — the piston is withdrawn for piston velocities X > 0) 
are entirely analogous; it is only necessary to change the sign of u (and X ) 
in each of the solutions. This reflects the fact that simple C' waves and 
simple C" waves differ only in direction of travel, i.e., the sense in which 
x is chosen positive. 
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On the Necessity of Unsteady Flow 
in Fluid Machines 

It is shown that work can be done on, or by, a fluid in a frictionless flow process only if 
the flow is unsteady. The implications of this conclusion in the analysis of fluid 
machinery , and provisions which must be made in analysis to account for unsteady 
flow, are discussed. 


T he purpose of this paper is to demonstrate that the 
flow in frictionless fluid machines must be unsteady if work is 
done on or by the fluid. This fact has great significance in any 
analysis of these devices, since unsteady terms must be included 
in the momentum, continuity, and energy expressions applied 
through the impeller, rotor, or cylinder of the machine. The un- 
steadiness governs the limits reached in certain idealizations of 
these devices. Also, in measurement, the recognition of the flow- 
unsteadiness is important in selection of instruments and inter- 
pretation of data; utilization of unsteady properties also might 
offer new measuring techniques. 

The analysis leading to equation (30) in this paper is not origi- 
nal. It may be found in German in “Technische Stromungslehre” 
by B. Eck. 1 However, as far as the author knows, the analysis 
is not available in English and its implications have never been 
discussed in the literature. 

In order to illustrate the importance of transient effects an 
apparent paradox will be described first. Then the resolution of 
this paradox will be given through a development of the relevant 
general theory of unsteady flow. 

Consider the flow of water through a steadily operating, hy- 
draulic turbine as shown in Fig. 1. 

This machine may be imagined to be a frictionless, reversible, 
adiabatic device without compromising the conclusions. Now, 
for the control volume shown a common result, from the steady- 
flow-energy equation of thermodynamics may be written 

]V X = — A(y/ijoU) (1) 

This equation states that the shaft work per pound of water, ii „ 
is equal to the decrease in head, g/goH, as the water flows through 
the turbine. The head of the fluid is defined as 


1 Julius Springer, Verlag, Berlin, Germany, 1957, p. 50. 

Contributed by the Fluid Mechanics and Compressor Subcommit- 
tees of the Hydraulic Division and presented at the Semi-Annual 
Meeting, Detroit, Mich., June 15-19, 195S, of Thf. American So- 
ciety of Mechanical Engineers. 

Note: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME Headquarters, April 24, 
195S. This paper was not preprinted. 



Fig. 1 Schematic of a water turbine 


g/0i =- + —-+ g/gcz ( 2 ) 

P 2(?o 

where 

H = head (in length units) 
p = static pressure 
p = density 

C = velocity (in Newtonian space 2 ) 
z = elevation above a datum plane 
g = acceleration of gravity 

go = proportionality constant in Newton’s law- 
inserting (2) into (1) gives 

I It = -A (- + + g/got) (3) 

\P 2g a / 

In order to combine the variables p, C, p, and z in (3), Bernoulli’s 
equation is utilized for this frictioniess, incompressible flow; it 
states that 

A ( — + + y/gez \ = Ap/po/f = 0 (-1) 

\ P 2 f/o / 

Upon substitution of (4) into (3) it is discovered that 

W x = 0 (5) 

In other words, the turbine can produce no shaft work. 

Experience indicates that this conclusion is wrong. \\ here is 
the error in analysis? 

Fundamentals 

It must be concluded from this conflict that either the steady- 
flow-energy equation (1) for the control volume or Bernoulli’s 
equation (4) is wrong or misapplied. Both equations will be re- 
examined starting from their restrictions of derivation 


Restrictions 

Steady flow energy equation 

Applied through this machine 

1 Incompressible 

yes 

2 Steady 

? 

3 Frictionless 

yes 

4 Adiabatic 

yes 

Bernoulli 


1 Incompressible 

yes 

2 Steady 

? 

3 Along a streamline 

? (what hap- 

pens as the 
blades chop 


up the . 
stream- 
lines?) 

4 Frictionless 

yes 


1 Any space where: F ila. 
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From this examination, the question of flow steadiness arises 
with both equations. 

In the sense of the flow-energy equation, the flow is steady, for 

/. A £k v 

no energy is accumulating inside the control volume ( i.e., — - — = 


0 J over measuring times of several revolutions. Therefore, re- 
striction 2 of the stoady-flow-energy-equation tabulation can be 
answered “yes,” clearing that equation and its result, that the 
shaft work is finite. 

But, in the sense of the Bernoulli equation, is the flow stead}'? 
The Bernoulli equation will be rederived in an attempt to answer 
this question. 

Consider a particle of fluid tf.-t in cross-section area, ds in 
length, and of unit depth as shown in Fig. 2. Newton’s law states 
that the sum of all the forces acting thereon in the s or streamline 
direction must be proportional to the particle’s mass times its ac- 
celeration in that direction or 


IF. = - Ma. (G) 

go 

The sum of the forces in the absence of shear is 

AF. = - dsJA - g/go ~ ( pdsdA ) (7) 

os os 

net pressure force gravity force 

The mass of the particle is 

M = pdsdA (8) 

The acceleration is equal to the particle’s time rate of change of 
velocity as we follow the particle, for Newton’s law applies only to 
a particle of fixed identity. A special derivative sign D called the 
“substantial derivative” is used to indicate that the particle is 
being followed or 


Since the particle velocity C is in the s direction, the acceleration 
in the .s direction a. is equal to the time derivative of C. Then 
combining (7), (8), and (9) 


Now. the particle’s velocity can be thought of in terms of a velocity 
field at each point of which the velocity is different at any in- 
stant and is changing with time. The particle is visualized mov- 
ing across this field ; over a length of time At the particle moves 
from point a to point b. The difference in the particle's velocity 
when at point b, as compared to its velocity when at point a, can 
be thought of, approximately, as the difference in the velocity of 
the field between b and a, at the instant of time when the particle 


is at point a, plus the change in field velocity of point b during 
the time it takes the particle to travel from a to h. In other words, 
indicating the particle velocities by C and the field velocities by 

C', 

c„- C a = C\_M - C' h = Q + (CV - CV) _ 0 (II) 

Dividing by the time At and replacing the second term on the 

right-hand side by — where As is the distance between a 

As 

and b 

ACft-o = A cy As ACV. 

At At At As 

Then letting At —*• 0 which causes a and b to be separated by an 
infinitesimal distance (for finite velocities) and, therefore, As — ► 0 

As 

and — — > ► C , there is in the limit an exact expression for the 
particle’s acceleration 

dC dC dC 

— = — + C — (13) 

dt dl ds 

But, since the first term on the right was the time rate of change 

bC 

of velocity at one point in space , b, the partial derivative — should 


be written instead of to indicate this fact. Likewise, since 

dl 

the second term on the right was written for a given instant in 

bC dC 1 

time (t = 0) this term should be written C rather than C — 

’ ds ds 

dC , 

to indicate this fact. Also the term on the let t is for the par- 

dt 

DC . 

tide which is given the special symbol — as mentioned before. 

Dt 

Since the special derivative D has been adopted, all other deriva- 
tive symbols can be stated to refer to the field and the primes on 
C can be dropped without much confusion. 


F or particle 
DC 
Dt 


For field 


at a point at an instant in time 


Combining (14) with (10) gives an equation relating changes in 
velocity C, density p, pressure p, and elevation z to position along 
tlie streamline s and the time t 


±-(* + C* 

go \ dl ds 


1 dp , dz 

H — ~ b g/go — = 0 

p ds ds 


Now, to obtain a closed-form relation between pressure, elevation, 
velocity, and density, these properties must be made dependent 
on only one variable. In general, the properties depend on three 
space co-ordinates and time. First, time will be eliminated, 
which makes dC/dl = 0, and all other properties independent of 
time. If, in addition, there is integration only along streamlines 
so that changes normal thereto do not appear, two space varia- 
bles can be eliminated reducing the properties to dependence on 
only one variable, s. In this condition, the partial derivatives 
can be replaced by total derivatives or 

1 „ dC 1 dp , dz 
— C + — + g/go ~ =■ 0 (16) 

go ds p ds ds 

Canceling ds, simplifying, and rearranging into familiar order 

— + + g/go dz — 0 (17) 

p 2ff 0 
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Euler’s Equation is obtained. If the density is constant, it is 
possible to integrate along the streamline immediately giving 

A ( v /p + 0725-0 + g/goz) =0 (18) 

which is Bernoulli’s Equation. 

From this derivation we learn that both Euler’s and Bernoulli’s 

equations are limited to applications where: (a) They are applied 

only along streamlines, ( b ) the fluid properties and velocity in New- 

( dp dp 
i.e., — = — = 
dt dt 

dC \ 

— = 0 1. The latter restriction defines steady flow as far as 
dt / 

Euler’s and Bernoulli’s equations are concerned. 

Now returning to the water turbine or any other fluid machine 
that does not depend upon friction forces to do its work, our ex- 
perience tells us that, if we stand and observe the rotor or piston 
at a point stationary relative to us in Newtonian space, the fluid 
properties and velocity will not be constant in time. As the blades 
of a turbine move past the stationary point, first the high-velocity 
low-pressure region on the blade’s suction surface passes, then the 
relatively lower velocity and higher pressure on the blade’s pres- 
sure surface. At a point on the casing above the tip of the rotor 
blades, the velocity and pressure violently fluctuate in time. So, 
too, do the velocity and pressure at a stationary point in a re- 
ciprocating compressor or internal-combustion engine, or in any 
other fluid machine that does not depend upon friction forces 
to do its work. 

Experience and the restrictions on Euler's and Bernoulli’s equa- 
tions lead to the conclusion that those equations cannot be applied 
through a fluid machine in which shaft work is done on or by the 
fluid. The generality of this statement will be proved, thus re- 
solving the paradox discussed in the introduction. 


where Q = heat. If in addition, the process is adiabatic 

dQ — 0 = du + pdv (22) 

Differentiating ho, yields from Equation (19) 

dho = du -| — dp + pdv H dC 2 -1 — — dz (23) 

P 2go go 

or combining with Equation (22) to eliminate du + pdv 

dh 0 = 1 dp + — dC 2 + — dz (24) 

P 2 go go 

Then, if it is desirable to look only at changes in the direction of 
the streamline, 

d/t 0 1 dp 1 dC 2 g dz 

— - = — -f + — — (25) 

ds p d.s 2<7 o ds go ds 

This equation may be combined with Newton’s L. . tyo) to give 
finally 

d/to 1 dC 

— = — for reversible, adiabatic flow. (2G) 

d.s go dt 

This equation states the surprising fact that, in a frietionless 
adiabatic flow (compressible or incompressible), the spatial rate 
at which the stagnation enthalpy (or head, in an incompressible 
flow) is changing in the direction of the instantaneous velocity at 
a point is proportional to the lime rale of change of fluid velocity of 
the field at that instant, at that point. In other words, if the flow is 
steady, dC/dt = 0; the stagnation enthalpy /i 0 is constant along 
streamlines, and the steady-flow energy equation 

~ IF* — ho, — ho, (27) 


Analysis of the Flow in a Fluid Machine 

Since steady flow, in the sense that the fluid properties and 
velocity at each point are invariable, is not a tenable assumption 
for analysis of fluid machines-, particularly since it would lead to 
the embarrassing conclusion that the shaft work is zero, the un- 
steadiness of the flow must be included in the analysis. Accord- 
ingly, relations between the work done on the fluid and the un- 
steadiness of the flow will be constructed. 

Stagnation enthalpy In, is defined in Newtonian space as 

C 2 p C 2 g 

ho = h + - — + y/goz ^ u -f — + + 2 (19) 

zjo P 2g 0 go 

where 

■u = internal energy of the matter in the absence of the effects 
of motion, gravity, capillarity, electricity, and mag- 
netism. 

There may also be written for any kind of process in a pure sub- 
stance 3 

Tds = du + pdv (20) 

where 

T = temperature, absolute 
s = entropy 

1 

v = - 

P 

and for a reversible (implies frictionless) process 1 

Tds = dQ (21) 


where the subscripts 1 and 2 represent two stations along the 
streamline, indicates again that the shaft work must be. zero. 

Equation (20) turns out to be rather difficult to use in practice 
since the direction s is constantly changing as the streamline lashes 
about in an unsteady flow. This difficulty can be avoided 
by dealing with the rate of change of a particle’s stagnation 
enthalpy as it is followed along. Such '‘following’' brings in again 
the substantial derivative D, which is expressed as a result of 
arguments similar to those leading to equation (IT- by 


Dho dh 0 ^ ^ dho 

Dt dt ds 


(2S) 


The second term on the right is given by Equation (26) while the 
first on the right is easily obtained from Equation (23) by con- 
sidering changes in time only at a point 

dho 1 dp C dC g dz 

— = - — + + — — ( 20 ) 

dt p dt go cV g., dt 

The dz/dt term will be ignored, since a point stationary in New- 
tonian space cannot change elevation z in the earth's gravita- 
tional field. Then combining equations (2S), (26), and (20), 

— 0 = - — — for reversible, adiabatic flow. (30) 

Dt p dt 

For an incompressible fluid, since 

h = g/goH (31) 

then 


D(g/goEl ) 1 dp 

Dt p dt 


3 See Keenan, “Thermodynamics,” p. So. 
1 Ibid., p. 83. 


( 32 ) 
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These expressions state that the rale at which a particle’s stagna- 
tion enthalpy (or head) changes with lime is proportional to the in- 
stantaneous rate of pressure fluctuation at the point where the 
particle finds itself at any given lime. 

Thus if work is to be done on the particle by pressure forces in 
order to change its stagnation enthalpy (or head), the pressure must 
fluctuate at a point and the flow is, as a consequence, unsteady , 5 
If shear, electrical, magnetic, or capillary forces are utilized to do 
work on the particle, the flow may be steady. 

Physical Interpretation 

If a hole were bored through a turbine casing directly above the 
rotor blade tips, and the pressure were measured there, something 
like Fig. 3 would be recorded. 

This pressure -ace results from the pressure field of the blades 
moving under pressure-tapping hole. Looking at the trace, it 
o 

can be seen tha - is negative. This fact, according to Equation 
dt 

(30), implies that the stagnation enthalpy of the fluid particles 
passing through the turbine is decreasing. This conclusion is 
known to be true for a turbine. 

In an axial compressor, the reader can construct a similar argu- 
ment and arrive at a conclusion exactly reversed from the turbine 
case above. 

To drive home the theoretical arguments given, a ma- 
chine, say a centrifugal blower, will he analyzed in just the 
reverse fashion to the preceding development. Imagine a blower 
doing work steadily on a flowing gas. The only moans that the 
rotor has at hand to impart those forces which do the work on the 
fluid is through a local pressure difference between the two sides 
of each impeller vane. For simplicity, assume that the pressure, 
high on tin; driving vane face and low on the trailing face, varies 
linearly in the tangential direction across the space between the 
vanes. Then, as the rotor turns by, stationary observers note 
that the pressure varies relative to them. In fact, if the pressure 
difference across the blade is Ap and the time bfrjweon vane pass- 
ing; is 



Fig. 3 Schematic of the pressure variation with lime at a casing pres- 
sure tap in a turbine 



MERIDIONAL VIEW PLAN VIEW 


Fig. 4 Centrifugal compressor impeller 

where A p is the pressure difference across a vane at ratlins r. The 
rate at which work is done on the particle by this torque is ri). 
This power applied to the particle equals the rate at which the 
particle's stagnation enthalpy (or head) is increasing in adiabatic 
flow according to the first law of thermodynamics or 

I)hu (Ap)rdsdnCl 
IV pdsdn 2~riZ 

which is divided by the particle mass, since ho is the stagnation 
enthalpy per unit mass. 

Then, on combining (30) and (34), it can be seen again that 

Dho __ 1 bp_ 

Dt p dt 



(33) 


where 

r = radius of the observation point from the center of rota- 
tion 

Z = number of blades 

V. = angular velocity of impeller in radians/sec 


This is exactly the expression derived before from more abstract 
(but more exact) considerations. Thus, if work is clone on or by 
the fluid, it must be exposed to boundaries which impose a net 
pressure force and which must move creating by their motion 
a pressure that fluctuates in time at a stationary point. 

Some Further Considerations 


then the rate of pressure variation dp/dt is 

dp A p (Ap)Z'fl 

dt At 2ir 


(34) 


Xow consider the mass particle of fluid ds, dn, shown in 

Fig. 4. The net torque about the axis acting upon this particle 
in the direction of impeller rotation is 


2t = (Ap)rdsdn 


1 The pressure must fluctuate at a point stationary in the co- 
ordinate system relative to which we measure work and stagnation 
enthalpy (or head). Conventionally, this would be relative to the 
machine casing. Relative to the machine rotor (I.e. to an observer 
riding on the rotor), if the relative flow is steady, then no work can be 
or is observed by the rotating observer. If the relative flow is un- 
steady (e.g., due to the disturbances of stationary upstream or down- 
stream blades, struts, and so forth), the rotating observer will see 
work done, although a different amount than the stationary observer. 
It is also interesting to note that the rotating observer would regard 
the stators as the work-producing elements of the machine. 


PT 'When idealizing turbomachincs, an infinite number of blades 
is often assumed in order to produce axisymmetric flow. It is im- 
portant to realize that, if the total amount and spatial distribution 
of the work done by the impeller is assumed to be unaffected by 
this idealization, then the rale of pressure variation at each point 
dp/dt must be unaltered according to Equations (30) and 1,3. I. 
In other words, the pressure variations must resemble Fig. o. 
The slope of the “saw teeth” must be independent of the number 
of blades. 6 As a consequence, although the flow becomes axisym- 
metrio for an infinite number of blades, it does not become steady. 
Euler’s and Bernoulli’s equations do not apply in the idealized 
impeller. 

The idea comes to mind, that the distribution of work addition 
or extraction within fluid machines could be determined by 
measuring dp/dt. Modern pressure probes make this measure- 
ment fairly easy. 

8 This follows also from simple reasoning, i.e., twice as many blades 
give one half the distance and one half the passing time between 
blades and one half the load or pressure difference across each blade. 
Therefore, bp/bt is independent of the number of blades. 
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p 




FEW BLADES MORE BLADES 


vmm. 


MANY BLADES 


Fig. 5 Schematic of the alteration of the pressure variation with time 
at a casing pressure tap of a fluid machine as the number of blades is 
increased. 


Summary and Conclusion 

It has been shown that: 

1 The flow in fluid machines doing useful work by pressure 
forces must be unsteady. 


2 Euler’s and Bernoulli’s equations are for steady flow and do 
not apply through the working region of the machine even if the 
events therein are cyclic. 

3 The steady-flow energy equation does apply over the 
working region if the processes are cyclic, for such are “steady” 
as far as that equation is concerned. 

4 The time rate of increase of a particle’s stagnation enthalpy 
(or head), on passing through the working region of the machine, 
is proportional to the local time rate of pressure fluctuation at a 
point stationary relative to the machine casing. 

5 Conclusion 4 can be explained readily in terms of the forces 
that moving blades or vanes exert on the fluid in order to do work. 

6 Idealizing the machine to an axisymmetric flow produced 
by an impeller with an infinite number of blades cannot make the 
flow steady. 

7 Measurement of the local flow unsteadiness might be a 
useful new technique to determine work distribution in impellers.. 
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WAVE AND THERMAL PHENOMENA 
IN H.S. -TUBES WITH AN AREA CONSTRICTION 
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The propagation of waves in a Hartmann-Sprenger tube 
having an area constriction is studied both theoretical- 
ly and experimentally. It is found that there are seve- 
ral possible modes of oscillations within the tube. The 
observed mode depends on the diameter - and length - 
ratios of the upstream and downstream parts of the tu- 
be. The mode also depends on the Mach number of the ex- 
citing jet. In the "first" mode, the shock wave genera- 
ted in the upstream cavity is transmitted trough the 
area change and strongly amplified. In higher modes, the 
incident wave in the upstream cavity meets, near the 
area change, a shock wave traveling in the opposite di- 
rection in the downstream cavity. For a given total 
length, the oscillation frequency of the 3rd and 5th 
inodes are approximately 3 and 5 times the frequency of 
the 1st mode, respectively. The thermal effects, essen- 
tially due to repetitive shock heating , strongly depend 
on the oscillation mode. In the 1st mode, the tempera- 
ture steadily increases along the tube, with a steep 
temperature gradient in the downstream cavity. In the 
3rd mode, a steep temperature gradient is observed 
near the area constriction whereas the temperature is 
high and almost constant in the downstream cavity. The 
influence of the shape of the joining section between 
the upstream and downstream cavities has been investi- 
gated. The highest shock intensities and thermal ef- 
fects have been observed with a sinusoidal converging 
section . 


1 - INTRODUCTION 

At the 13th International Symposium on Shock Waves and Tubes, it was shown* 
that a Hartmann-Sprenger tube (HS-tube) presenting a sudden change in cross 
section could produce high amplitude pressure waves. The configuration studied 
is shown in Fig. 1. A high speed jet flows out of a "needle" nozzle^ and drives 
flow oscillations within a tube with constriction. A simplified wave diagram 
for the flow within the tube was proposed. A theoretical amplification was de- 
rived and experimental values were in good agreement with the predicted ones. 
Subsequently, it was found that several modes of oscillations are possible-. 

The wave diagram as well as the thermal effects are quite different for each mo- 
de. This is important when one considers the H.S. -tube as an igniter for pyro- 
technic compositions^. 
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FACTORIES BLOW 
SMOKE RINGS 

Factory chimneys of the fu- 
ture will be shorter and fatter, if 
an American scientist has his 
way. Instead of dribbling 
smoke into the air in a continu- 
ous stream, they will puff 
smoke rings, like a contented 
cigar smoker. Dr. Timothy Fohl, 
a physicist from Boston, has 
been studying industrial smoke 
production and air pollution 
problems for some years. He 
observed that smoke rings trav- 
el much more quickly through 
the air than "unstructured 
smoke emissions" and sought to 
apply this principle to the rapid 
dispersal of industrial smoke. 
Dr. Fohl has calculated that a 
short, relatively plump factory 
chunnev equipped with a shut- 
ter and pump mechanism to 
putt out smoke rings will push 
the smoke to heights of 10, ODD 
feet before it disperses. The 
tallest conventional chimnev 
does not send smoke higher 
than 2500 feet, he savs. It is also 
reckoned that chimnev blown 
smoke rings will be less affected 
by wind than are ordinary 
smoke plumes . — Compressed 
Air, May 19n8. 
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Postfix 

BY ALEXANDER ROCA 

Killer Air Ray 


M ammoth smoke 

RINGS DOWN EN- 
E M V BOMB E R S ! 

M O N S IT. R MA- 
i HIM CREATES 
MAN-MADE TOR- 
NADOES! It was a fantastic idea, 
and it approached reality during 
World War II. The client? The 
U.S. government. 

In recent decades wingtip vor- 
tices generated by jumbo jets 
have been known to disrupt the 
flights of smaller aircraft follow- 
ing nearby, turning them upside 
down in seconds. Imagine hostile 
pilots in close formation and 
laden with armaments and explosives encountering such 
whirlwinds. 

That was the goal of a wartime plan hatched by the 
U.S. Army Chemical Warfare Service (USACWS) at the 
Dugway Proving Ground in Utah. The department did 
research not only on chemical agents but also on the means 
of delivering them. That’s where vortices came in, because 
they can improve the penetration of a gas charge through 
the air. A Colonel Esmond of the USACWS hired Thomas 
Shelton to work on the idea with Philip Leighton of 
Stanford University. Shelton was a genius inventor and 
aerodynamicist noted for his twin-boom, twin-engine 
1935 Crusader, an aircraft that cruised at 200 miles per 
hour on just 312 horsepower. Its efficiency stemmed from 
Shelton’s careful wind-tunnel-tested streamlining. At cruis- 
ing speed the teardrop-shaped cabin alone generated sub- 
stantial lift. And fuselage-to-wing filleting prevented the 
interference drag rhat plagued other, later twin-boomed 
aircraft. (Eor example, early Lockheed P-38 Lightnings 
suffered from tail buffeting until fairings were added to their 
wing roots.) 

Shelton quickly developed a clever vortex maker. A bell- 
shaped chamber, acted on by a piston or an explosion, 
would send a packet of gas out its narrow end as a high- 
velocity vortex. The vortex enhanced the speed, distance, 
accuracy, and penetration through crosswinds of any 
vapor being shot out of the bell, just as rifling improves 
bullet performance. It could be used with vapor-based 
herbicides and insecticides, as well as with chemical and 
biological weapons. 

Shelton envisioned another way to use Ins device: Vor- 
tices from a huge generator might disrupt aircraft in flight, 
and the equivalent of a giant smoke ring might break 


up high-altitude enemy bomber 
formations. To test his thesis, 
he built a scale model. Powered 
by a shotgun shell, it shot an 
18-inch smoke ring 150 yards. 
And it made an eerie howling 
sound — potent psychological 
weaponry. 

According to Shelton, the 
Army Air Corps brass at Wright 
Field, in Dayton, Ohio, seriously 
considered a version with high- 
explosive power and a business 
end yards wide. It would rocket a 
deadly accurate and invisible ring 
of air at near-supersonic speed 
thousands of feet skyward. There 
the vortex would wreak havoc, spinning target aircraft 
wildly out of control, crashing into other craft, and falling 
to earth. 

This top-secret stratagem died unfulfilled as offensive 
weapons garnered war-research resources. But after the 
war Shelton patented miniature versions of his vortex gen- 
erator as children’s toys. The Flash Gordon Air Ray Gun, 
manufactured by the Budson Company of Chula Vista, 
California, knocked over targets with a puff of air and 
cost only $2.49. It was named “Toy of the Year” in 1949 
by Popular Mechanics magazine. 

Later, with backing from an investor named Paul Zifrin, 
Shelton developed the Magic Ray Gun, for the Nu-Age 
Company. It shot miniature smoke rings from small pellets 
developed by the Lion Match Company. The pellets, in 
book-match form, neither flamed nor flared but safely 
smoked prodigiously. Shelton based his idea on a toy made 
from a cylindrical oatmeal carton by his father, who was 
himself an industrialist and inventor of no small repute. 
When filled with cigarette smoke and tapped on one end, 
it would eject a ring from a hole in the other end. In the 
1960s Wham-O Manufacturing revived the Budson gun as 
the Air Blaster, cleverly configured with targets built into 
the cardboard box it came in. 

Today “death ray” fantasies are back as prototypes 
for ground- and satellite-based laser and particle-beam 
weapons. Tom Shelton’s idea — which also projected 
disruptive, destructive energy — looked forward almost 
half a century. Perhaps it was merely too far ahead of 
its time. ★ 

Alexander Roca is the author of Crusader: The Story of the 
Shelton Flying Wing (Rare Birds Publishing, Berlin, Mass.). 



The end result of the Army’s vortex-weapon research. 
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Ring Bubbles 
of Dolphins 

A number of bottlenose dolphins 
in Hawaii can create shimmering, stable 
rings and helices of air as part of play 


by Ken Marten, Karim Sharif Suchi Psarakos and Don J. White 



B elow the towering cliffs of Ma- 
kapuu Beach on the island of 
Oahu, Hawaii, is a unique lab- 
oratory dedicated to the study of dol- 
phins. Project Delphis, run by the non- 
profit conservation organization Earth- 
trust. in cooperation with Sea Life Park 
Hawaii, conducts research ranging from 
investigating dolphin self-awareness to 
exploring the animals' intelligence using 
an underwater computer touch screen. 
The scientists in the lab do not use food 
as a reward, so all the behavior observed 
is of the dolphins’ own volition. 

One of the most fascinating activities 
we have seen in our research involves 
no high-tech human toys at all. Instead 
the dolphins fashion their own enter- 
tainment by swirling the water with 
their fins and blowing bubbles into the 
resulting vortices to produce rings and 
helices of air. Furthermore, the physics 
behind the air rings turns out to be quite 
interesting. Few people doubt that dol- 
phins are highly intelligent animals, but 
these observations demonstrate just how 
imaginative they can be. 

As air breathers, dolphins blow bub- 
bles whenever they exhale underwater. 
Dolphin behaviorists have noted that 
when dolphins are excited, surprised or 
curious, they will sometimes expel air 
from their blowholes, generating large, 
amorphous bubbles that rise quickly to 
the surface. The animals occasionally 
emit streams of small bubbles when they 
make sounds; the bubbles might add 


KAIKO'O (left), an adult bottlenose dol- 
phin at Sea Life Park Hawaii, generates a 
rising ring of air to play with. The labora- 
tory at the park (above, right) offers a 
window into the world of dolphins. 


another component, detectable by sight 
or sonar, to the vocal message. 

Dolphins can also create more exotic 
types of bubbles for less prosaic reasons. 
In recent years, researchers at several 
oceanariums around the world have re- 
ported that a variety of marine mammals 
can blow smooth, stable rings of air that 
linger in the water for several seconds. 
Because of the intricate techniques and 
practice required to form such rings, as 
well as the helices we have seen, these 
bubbles are clearly not a spontaneous 
response to alarm or a standard part of 
communication. Wolfgang Gewalt of the 
Duisburg Zoological Gardens in Germa- 
ny observed untrained Amazon river dol- 
phins (Inia geoffrensis) producing bub- 
bles in unusual and playful ways. The 
animals emitted air from their mouths 
to yield necklaces of bubbles, which they 
would pass through or bite. 

Elsewhere, at Marine World Africa 
USA in Vallejo, Calif., Diana Reiss and 
Jan Ostman-Lind noticed bottlenose dol- 
phins ( Tursiops trancatus ) at the aquar- 
ium playing with rings similar to the 
ones we have seen. Kenneth S. Norris 
of the University of California at Santa 
Cruz described beluga whales ( Delphi - 
naptenis leucas) at the Vancouver Aquar- 
ium that expelled bubbles from their 
blowholes and then sucked them into 
their mouths as part of playful behav- 
ior. Some people have even witnessed air 
rings in the wild. The behavioral biolo- 
gist Karen Pryor observed male Pacific 
spotted dolphins ( Stenella attenuata) 
blowing rings during social encounters; 
Denise Herzing of Florida Atlantic Uni- 
versity has viewed similar displays in 
free-ranging Atlantic spotted dolphins 
(Stenella frontalis). And the marine 
mammal photographer Flip Nicklin has 


seen beluga whales living in Lancaster 
Sound in Baffin Bay, Canada, that re- 
leased rings of air as they repeatedly 
clapped their jaws together, possibly in 
a display of aggression. 

Hawaiian Ring Culture 

D uring the past five years at Sea Life 
Park Hawaii, we have studied 17 
bottlenose dolphins; nine of them, rang- 
ing in age from 1.5 to 30 years old, gen- 
erated air rings. Based on our observa- 
tions at other oceanariums and consul- 
tations with colleagues at various sites, 
we believe ring blowing is more com- 
mon at Sea Life Park Hawaii than at 
other aquariums; the dolphins here ap- 
pear to have created a “ring culture” in 
which novice dolphins learn to make 
rings in the presence of experts that, in 
a sense, pass down the tradition. 

Ring making is a leisurely pastime, so 
the animals generate rings only when 
they want to— not on command or for a 
reward of food. Furthermore, ring mak- 
ing does not seem to be associated with 
functional behaviors such as eating or 
sexual activity. Because ring blowing is 
unpredictable, we have had some diffi- 
culty documenting it. But over time we 
have captured most of the dolphins at 
play in photographs and on videotape. 
(Additional photographs can be seen on 
the World Wide Web at http://earthtrust. 
org, and a Quicktime movie can be 
found at http://www.sciam.com/) Un- 
fortunately, neither of these media does 
the dolphins justice. But the archived 
images do reveal important information 
about the physics of the rings. 

From what we have seen, the dolphins 
employ three basic techniques to form 
the rings. In rhe simplest method (also 
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used by human divers), dolphins puff 
our bubbles from their blowholes; these 
bubbles become halos of air that expand 
in radius while decreasing in thickness 
as they rise to the surface [ see illustra- 
tion above]. One of the dolphins we 
watched, the adult male Kaiko'o, could 
emit two rings in succession and then 
fuse them into a single, large ring. 

The physics behind this type of ring is 
relatively straightforward: any spherical 
bubble bigger than about two centime- 
ters in diameter will quickly become a 



VORTEX 


HALO OF AIR becomes thinner and expands 
in radius as it is carried to the surface by the 
vortex flow ( arrows ) through the center of the 
ring. The vortex also serves to stabilize the ring, 
preventing it from breaking up into smaller 
bubbles. Here Kaiko’o produces a ring and fol- 
lows it up through the water ( video stills, bot- 
tom to top). 


... 


X 



ring because of the difference in water 
pressure above and below the bubble. 
Water pressure increases with depth, so 
the bottom of the bubble experiences a 
higher pressure than the top does. The 
pressure from below overcomes the sur- 
face tension of the sphere, punching a 
hole in the center to create a doughnut 
shape. 

As water rushes through this hole, a 
vortex forms around the bubble. Any 
vortex ring travels in the same direction 
as the flow through its center; in the 
case of these simple air rings, the vortex 
flow, in combination with the air’s nat- 
ural buoyancy, propels the bubbles to- 
ward the surface. Although the process 
of making air-filled rings is fairly sim- 
ple, dolphins cannot blow stable ones 
without some practice, suggesting that 
additional factors— such as the viscosity 
of the water and how the air is ejected 
from the blowhole — need to be taken 
into account. 

In a more elaborate approach, the dol- 
phins fabricate rings that travel horizon- 
tally and sometimes even downward in 
the water. For instance, a dolphin might 
swim forward rapidly on its side so that 
its normally horizontal flukes (that is, 
its tail fin) are vertical. By thrusting its 
flukes vigorously to one side, the animal 
generates an invisible, ring-shaped vor- 
tex that travels horizontally and slight- 
ly downward. After quickly turning 
around, the dolphin finds the vortex and 
injects a bubble into it from the blow- 
hole. (The dolphin often produces an au- 
dible series of clicks before the release of 
air, suggesting that sonar may be em- 
ployed to locate the invisible vortex. ) The 
pressure inside a vortex is lowest in the 
center, or “eve, v of the sw irl; when the 
dolphin exhales into the vortex, the air 
migrates to the region of lowest pres- 


sure and is drawn out along the core of 
the ring-shaped vortex [ see illustrations 
at top of pages 86 and 87]. 

The resulting ring can be up to 60 cen- 
timeters (two feet) across and just over 
a centimeter thick, traveling horizontal- 
ly in the water. Once again, the move- 
ment of the ring reflects the direction of 
the flow through its center; in the case 
of the vortex created by the dolphin tail 
fin, this flow is horizontal and some- 
times even downward. With a sideways 
toss of its rostrum, or jaws, the dolphin 
can pull a small ring off the larger one 
and then steer it through the water. 

Making a vortex with flukes and 
planting a bubble in it can be done in a 
variety of ways — the adult female Laka 
often positions herself vertically in the 
water with her head pointed down and 
tail up. Laka then flips her tail fin down 
to stir up a vortex. She fills the vortices 
not only with air from her blowhole but 
also with air from her mouth. Occasion- 
ally, Laka will capture air from above 
the water with a downward thrust of her 
flukes. In addition, we have watched 
Laka release from her mouth small bub- 
bles that pass along her body; when the 
bubbles reach her flukes, she flips them 
into a ring. And she can even augment 
a ring by injecting more air into it. 

Experimenting with Bubbles 

T he third type of air-filled vortex 
dramatically reveals the dolphins’ 
capacity for experimentation. On a few 
occasions, we watched the young fe- 
male dolphin Tinkerbell, Laka’s daugh- 
ter, construct long helices of air, using 
the most complicated technique we have 
seen. These more complex structures no 
doubt result from considerable refine- 
ment through trial and error. Tinkerbell 
has developed two very different meth- 
ods for making helices. In one approach, 
she releases a group of small bubbles 
while swimming in a curved path near 
the wall of the tank. She then turns 
quickly, and as the dorsal fin on her hack 
brushes past the bubbles, the vortex 
formed by the fin brings the bubbles to- 
gether and coils them into a helix three 


84 Scientific American August 1996 


Ring Bubbles a* Dolphins 



[cT 







MARTEN, PSARAKOS AND WHITF: 



to rive meters (10 to 15 feet) long [see 
illustrations at bottom of these two 
pages]. We have also seen Tinkerbell 
swim across the tank in a slightly 
curved path, leaving behind an invisible 
dorsal fin vortex. She then retraces the 
path and injects a stream of air into the 
vortex, producing a long helix that 
shoots out in front of her. 

Again, because the pressure in a vor- 
tex is lowest in the eye, once the bub- 
bles are inside the vortex, they move to- 
ward the center, merge and elongate 
into a helical tube. Usually a tube of air 
in water is unstable and breaks up into 
smaller bubbles. But all the dolphins’ 
rings and helices are shiny and smooth 
because the variation of pressure inside 
the vortex (low pressure at the center, 
building up to higher pressure at the 
edges) actually works to stabilize the tube 
by smoothing the ripples that would 
otherwise break up the large bubble. 

Although we did not see any dolphins 


other than Tinkerbell generating helices, 
the practice of making rings spread 
through the population of dolphins, as 
some of the individuals learned the 
technique in the presence of their ring- 
blowing companions. We had the op- 
portunity to watch one young dolphin’s 
rings evolve over a period of two months 
from unstable, sloppy bubbles that dis- 
sipated rapidly to stable, shimmering 
rings that lingered in the water for sev- 
eral seconds. Older dolphins also need- 
ed time to acquire the talent. One adult 
male, Keola, lived in the research tank 
for two years with dolphins that did not 
produce air rings, and during that time 
we did not see him generate any. But 
when his younger, ring-blowing sibling 
Kaiko’o moved into the same tank, Ke- 
ola watched for long stretches while 
Kaiko’o blew rings; within a couple of 
months, Keola began making his own 
rings, which slowly progressed in quality. 

We have noticed that other dolphins 


also monitor their ring-blowing tank 
mates intently, suggesting that the exhi- 
bition interests the animals or offers a 
learning opportunity for them. On sev- 
eral occasions we saw the two brothers 
Keola and Kaiko'o lying side by side on 
the bottom of the tank, repeatedly blow- 
ing large doughnut rings either simulta- 
neously or within a second of each other. 
We have also seen one female, swimming 
closely behind another female who was 
blowing rings, produce her own bursts 
of small bubbles as she watched. 

The dolphins have drawn humans into 
their play as well: one day during a pe- 
riod of intense ring making, Tinkerbell 
repeatedly blew a ring and then came to 
the lab wrindow where one of us (Psara- 
kos) was videotaping, as if to include her 
in the activity. Once, we blew' soap bub- 
bles inside the lab in front of the dol- 
phins’ wfindow', and within a few min- 
utes one of the dolphins joined in by 
blowing simple, rising doughnut rings 




CORKSCREW'S OF AIR are uncommon: the authors have ob- 
served only Tinkerbell ( video stills, above) create them. In one 
technique, Tinkerbell emits a burst of bubbles while producing 


a spiral vortex off the dorsal fin (a). When the bubbles en- 
counter the vortex {b), they are drawn into it, then merge and 
elongate into a long helix of air (c). 
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HORIZONTALLY MOVING RING forms 
when a dolphin flicks its tail while swimming 
on its side. The horizontal and slightly down- 
ward motion of the fin creates a vortex travel- 
ing in the same direction (a). The dolphin then 
turns around and injects air into the swirling 
flow ( b ); the air is drawn along the core of the 
vortex (c), forming a ring that moves in the di- 
rection of flow through its center. The adult 
female Laka is pictured (left) exhaling into the 
vortex and examining her creation. 



near the lab window. The real surprise 
came when the dolphin swam away 
from the window and made several 
fluke vorrex rings — so different from 
what we each blew at the window. 

Our study of the bottlenose dolphins 
at Sea Life Park Hawaii continues in 
the hope of better understanding their 
behavior. As the only nonprimates that 
have shown strong indications of self- 


awareness, these dolphins may teach us 
about the nature of intelligence through 
their experimentation and play. But as 
we consider their remarkable abilities, 
we are haunted by the knowledge that 
many cultures, including our own, re- 
gard dolphins as expendable. Dolphins 
continue to be targeted by tuna nets, to 
become ensnared in expansive drift nets 
and gill nets, to be canned as mock 


whale meat and to be shot for crab bait 
or to be hunted for sport. Earthtrust and 
its sponsors work to address these is- 
sues, but we believe only a basic change 
in human behavior will make a perma- 
nent difference. It is our fervent hope 
that by providing new views into the 
dolphin mind, we may yet convince peo- 
ple to stop the indiscriminate slaughter 
of these fascinating creatures. Q 
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central portion of the air column continues to pass 
it, thus maintaining the rotary motion. 
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MECHANICS AND APPLICATIONS OF SELF RESONATING CAVITATING JETS 

G. L. Chahine and V. E. Johnson, Jr. 

Tracor Hydronautics. Inc. 

Laurel. Maryland 


ABSTRACT 

Self- resonating jets take advantage of the natural 
tendency of an axisynmetric jet to organize into large 
structures. If a feedback mechanism at the nozzle exit 
isprovided, the energy content of the jet natural 
frequency can be dramatically amplified with a resonant 
chamber. The submerged jet shear layer then organizes 
in large ring vortices anitted from the nozzle at. a 
discrete frequency and could cavitate to font toroidal 
bubbles. These cavities grow, collapse, and interact 
with a nearby impact surface and generate very high 
pressure fluctuations. Associated with these 
fluctuations are the usually undesirable phenomena of 
noise and erosion. These Effects can, however, be 
taken advantage of in important applications such as 
deep hole drilling, cleaning, cutting, underwater sound 
generation, .. .etc. In deep hole drilling, jet 
structuring is shown to improve on conventional 
cavitating jets by increasing the cavitation inception 
number, aiding in chip removal, and improving drilling 
rates. Cavitating jets are known to be noisy. 
Self-resonance of the jet makes these fluctuations not 
only two to three orders of magnitude higher but also 
concentrates the energy into discrete frequencies. 

This makes self-resonating cavitating jets very good 
candidates for controlled high intensity discrete 
frequency underwater noise generators. 

In order to understand the mechanics of these 
self-resonating cavitating jets, several fundamental 
studies have been conducted and are still on-going. 
These cover the areas of nozzle design, flow field 
description, and jet nozzle, jet target interaction. 

In this paper, the principle of operation of these jets 
will be described. The importance of the various 
parameters will be discussed, and results of tests as 
related to the various aspects presented. Analytical 
and numerical studies related to the bubble ring 
dynamics and interactions are presented. Finally, jet 
simulation by the replacement of the shear layer with a 
succession of ring vortices is shown to allow the study 
of important effects such as: jet structuring; jet 


shear layer - nozzle lip interaction; and jet wall 
interaction. 

I NTRODUCTI ON 

Cavitation is mainly known for its harmful 
effects, namely, loss of performance, erosion, and 
noise. The usual procedure to prevent these 
deleterious effects is to avoid the phenomenon oy 
proper design and by limiting the operating 
conditions. However, attempts to induce and harness 
cavitation for useful purposes have been increasingly 
successful. Ultrasonic cavitation methods take 
advantage of the erosive power of cavitation for 
cleaning, emulsification, and mixing. In high-pressure 
jets, cavitation has for seme time row been purposely 
induced in order to increase their drilling, cutting, 
and cleaning capabilities, (Johnson et al, 1972). The 
noise associated with cavitation is being used as a 
means of sensing cavitation when it becomes destructive 
and, hence, could allow for alleviating its damaging 
effects. More recently, ws have developed a more 
direct utilization of cavitation roise: underwater 

sound generation (Chahine, Johnson, Lindermuth, and 
Frederick, 1985). 

Experimental observations of submerged jets show 
the tendency of the turbulent eddies in their shear 
layer to organize in large structures. Excitation of 
an air jet with periodic acoustic signals produced 
upstream of the nozzle by transducers or loud speakers 
shows a remarkable change of the jet structure into 
discrete ring vortices when the excitation frequency, 
f, matches the predominant natural frequencies of the 
non-excited jet (Crow and Champagne, 1971 ard Kibens, 
1979). This corresponds to a Strouhal number, Ey, 
close to 0.3 or one of its first integer multiples. V 
and d being the velocity and the diameter of the jet, 
the Strouhal number is defined as 

S d = fd/V, (1) 



129 


This "natural tendency" of a submerged jet to organize 
in large structures is of great interest in 
aerodynamics for air jet noise studies. Crow *d 
hampagne (1971) studied this phenomenon ext' -ely 
and showed experimentally that forced excitai . <n of the 
jet at the preferred frequency (S^O.3) enhances the 
structuring. The vorticity is then mainly concentrated 
in ring-shaped, large structures. 

The potential of such a phenomenon for submerged 
water jets was recognized (Johnson et al, 1984), and a 
unique technology was developed related to achieving 
useful submerged jets having very high amplitude, 
periodic, oscillatory discharge without using moving 
parts in the supply system. The passive excitation is 
obtained hydroacoustically and structures the shear 
layer of the jet into discrete, well defined ring 
vortices when the excitation frequency, f, matches the 
jet's preferred value. This can be obtained by (a) 
feeding the final jet-forming nozzle with various types 
of acoustic chambers (for example, Helmholtz chambers 
or organ-pipe tubes) tuned to resonate at the desired 
frequency; and (b) shaping the nozzle so as to feed 
back the pressure oscillations which occur at the 
exit. Such devices are forms of "whistles" which 
self-excite and thus are totally passive. These 
STRATOJETs have shown enhanced erosivity frcm increased 
cavitation activity. The large pressure oscillations 
associated with the intensification of cavitation, with 
resonance in the nozzle assembly, and with the 
production and disappearance of large vortical 
structures greatly improve the erosion and cleaning 
capabilities of the jet and make a very effective noise 
generator. 

The principal types of STHATOJCT configurations 
that we have studied and developed are shown in Figure 
1. The PULSER (Fig. 1-a) and PULSF.R-FFD (Fiq. 1-c) 
configurations consist of a Helmholtz oscillator whose 
volime is tuned to coincide with the jet preferred 
Strouhal number. Another type of STRATOJET which has 
been the subject of extensive study and (development at 
Tracor Hydronautics uses an organ-pipe acoustic chamber 
as illustrated in Figure 1-b. In this device the 
organ-pipe resonant frequency is selected to match the 
desired structuring frequency defined by the critical 
Strouhal number of the jet. This concept offers the 
simplest system design and has been used successfully 
for erosion studies and noise generation. Although 
PULSER and PULSER-FED nozzle systems have many 
important uses and continue to be studied at Tracor 
Hydronautics, we limit our discussion in this paper to 
the Organ-Pipe STRATOJET. 



• . "PULSER" b.“ORCAN PIPE" c. "PULSER-FED" 

Figure 1 - Self-Resonating Muzzle System Concepts 


PRINCIPLES OF OPERATION OF AN ORGAN-PIPE STRATOJET 

Acoustic resonance is achieved in the nozzle 
feed-tube assembly when a standing wave forms in the 
"organ-pipe" section (length: L, diameter: D in Figure 
1-b]. This section is created by the upstream 
contraction, (Df/D) and the nozzle contraction 
(D/d). Peak resonance will occur when the fundamental 
frequency of the organ-pipe is near the preferred jet 
structuring frequency. The exact resonance frequency 
is dependent on the contractions at each end of the 
organ pipe; for instance, if both (Dp/D) and (D/d) 
are large, then the first mode resonance in the pipe 
will occur when the sound wavelength in the fluid is 
approximately four times L. When D/d is close to one, 
resonance occurs when the wavelength is approximately 
two times L. Although, there are interesting 
applications where quarter wave organ pipes are 
applicable, most of our studies to date have involved 
half wave length systems. 

The major problem when designing an ORGAN-PIPE 
system is the generation of the initial excitation of 
the oscillations in the organ-pipe chamber. This is 
successfully achieved through a feedback mechanism at 
the exit section of the nozzle. Empirical acoustic 
analysis (Chahine, et al., 1984) have led to the 
following approximation, useful for estimating the 
length of the orqan pipe: 

L/d = K j/MSj (1+8) (2) 


where the "mode parameter" K n is given, for Dp/D>1, 
by 

( 3 ) 


(4) 


In these expressions, n is the mode number of the organ 
pipe, Sy the preferred Strouhal number, M the Mach 
number (ratio of jet velocity to sound speed) , and 8 is 
an erd correction. Figure 2 represents some experi- 
mental results for half-wave modes and the relation 




where 8 has been evaluated. 

The principles of operation and noise generation 
of an ORGAN-PIPE STRATOJET are schematically 
represented in Figure 3.TVo predominant sources of 
pressure fluctuations can be distinguished in addition 
to the classical nonexcited turbulent shear layer 
between the jet and the surrounding liquid. One of 
these sources is easy to define and corresponds to the 
volume fluctuations of the moving vortex bubble rings 
formed in the center of the large structures of the 
self-excited jet. The other source of pressure 
fluctuations is more complex and relates to the exit 
area of the jet where high amplitude oscillations of 
the main flow characteristics are interrelated with the 
shear layer-nozzle lip interaction. Basically, the two 
sources of pressure fluctuations are simultaneously 
sources of jet excitation and noise generation. The 
acoustic signals fran both areas are forcing functions 
to the resonating chamber in the nozzle assembly. 
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rate of penetration. The actual drilling rate is 
always equal to or less than the rate which can be 
achieved if chips are removed as scon as they are 
created. Chips are removed from the bottom of the hole 
by the drilling fluid. Most roller bits utilize the 
hydraulic horsepower available by ejecting the fluid at 
the maximum velocity possible from single orifices 
located in the sectors between the roller cones. 

Further increases in the jet velocity would improve the 
cleaning and thus the drilling rate; but, for the most 
part, the problems involved in working with higher 
pressure equipment have postponed development of this 
option. Using cavitating jets is a very effective way 
of significantly improving cutting efficiency without 
expenditure of additional power. However, with 
increasing depth, cavitation becomes harder to achieve. 

A. Effects of Jet Structuring on Cavitation Inception 
~jhe parameter defininq cavitation is the 
cavitation number, o. 


REFLECTOR 


I 

I 

I 

I 



Figure 3 - STRATOJET-Noise Generator. Schematic of 
Cperation 

These signals strongly interact; they are both fed back 
and amplified by the organ pipe of the STRATOJET. 

CUTTING AND CLEANING APPLICATIONS 

It is generally accepted in the deep-hole drilling 
industry that bit hydraulics has a major effect on the 


where p a is the cmbient or far field pressure, Py 
is the vapor pressure of liquid, p is the liquid 
density, and v is the characteristic velocity, the jet 
mean velocity. In the case of high-pressure submerged 
jets used in deep-hole drilling, p a >> p^, and for 
well-designed nozzles 1/2pV 2 may be approximated by the 
pressure drop, Ap, across the nozzle. Thus 


Ap 


( 7 ) 


The particular value at vhich cavitation is 
incipient is defined as 


i Apat inception 


( 8 ) 


Thus, if the operating conditions for a submerged jet 
are such that o/o^<l then cavitation will occur and 
as o/oi continues to decrease below unity the amount 
of cavitation will increase. When a cavitatiang jet is 
impinged against a surface, so that the cavities formed 
by the jet are collapsed on that surface, very high 
pressures (>>l/2pl/ 2 ) are generated, and the resulting 
cleaning or cutting action is substantially greater 
than when the jet is not cavitating. In deep-hole 
drilling, the ambient pressure is hydrostatic and P a 
is directly related to hole depth. 

Figure 4 indicates zones of cavitation for jets in 
drill bits under various operating conditions. Fbr a 
given fluid density, the ordinate may be related to 
hole depth. In Figure 4, the right side ordinate is 
shown for a drilling mud density of 10 lb/gal (1.2 
gro/cm 3 ). For typical drill bit nozzles, o is 
approximately 0.7. Therefore, cavitation will take 
place at an operating pressure drop of 2000 psi (13.8 
MPa) at depths less than about 2500 ft (760 m). If 
oj^ is multiplied by a factor of three, as is the case 
with the STRATOJET, the depth over which cavitation 
contributes to improved drilling rate becomes 
significant even at current rig mud -pump operating 
pressures. Figure 5 shows for various nozzle 
geometries the improvement in cq as a function of the 
relative pressure fluctuations, p'/AP, at mid-length of 
the ORGAN PIPE. 
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NOZILt I'lUSSUKI IlHOI'. /,(), pM 

Figure 4 - The Influence of on the Possibility of 
Cavitation Utilization for Deep-Hole 
Drilling 



Figure 5 - Gain in STRATOJET Cavitation Inception Index 
as Compared to Conventional Nozzles as a 
Function of Relative Pressure Fluctuations 
at Organ Pipe Mid-Length 

There are two reasons why structured jets should 
be of interest in the context of improving drill bit 
performance. Both of these reasons may be revealed by 
estimating the expected differences between an 
unstructured and a structured jet in local pressures in 
the flow field and at the surface on which the jets are 
impinging. 

Figure 6 shows the velocity distribution in a line 
vortex having a rotational core radius denoted as r c 
and a velocity at r c equal to V c . Such a Rankine 
vortex is a reasonable approximation of vortices which 
exist in real, viscous fluids. For such a single line 
vortex, the value of the pressure drop from the ambient 
pressure, p a , to the minimum pressure, p„i n , which 
exists at the center of the core (and reaches the vapor 
pressure, py, at inception) is 


1 2n 2 V 2 r 2 Oi '( vortex center) 

jpV 2 c 

where r is the circulation around the vortex. 

Since oi is desired to be as high as possible in 
order to cause increased cavitation and erosion, it is 
preferable for a given nozzle, liquid and speed to have 
r as large as possible. For unexcited jets, the shear 
zone has many small vortices ( r is small and of order 
of the product jet velocity by turbulence length-scale) 
which roll up into a statistical distribution of sizes; 
whereas, for an excited jet that becomes structured 
into discrete ring vortices, r is of the order of the 
product jet diameter by jet velocity. 

Hie preceding analysis is not exact because of the 
various simplifying assumptions made, (for example, the 
detailed pressure distribution in a ring vortex system 
is more complex) but the important result shown is 
that, qualitatively, oj- structured is much greater 
than oj-nonstructured. 



Figure 6 - Distribution of Tangential Velocity and 
Pressure Associated with a Vortex 

B. Effects of Jet Structuring on Hole Bottom Cleaning 
The normal procedure to resist "blow outs" in bit 
holes is to operate with an overbalance in pressure. 
That is, the mud density is selected so as to produce a 
hydrostatic pressure at the hole bottom greater than 
the formation pore pressure (Johnson, et al., 1984). 
Thus there is a flow of fluid ; nto the formation and a 
resulting "filter cake” buildup across the hole 
bottom. Consequently, the fluid passing over the hole 
bottom to remove chips must "lift" them against a very 
large pressure gradient - creating a hold-down force 




(15) 


which is many times larger than the force caused t 
gravity. Drilling rate decreases dramatically . e 
pressure difference across the hale bottom is 
increased. The formation of large structures in the 
STRATOJET (vortex rings) has the advantage of 
providing, even in the absence of cavitation, the 
needed lift force. The passage of each ring vortex 
across the hole bottom generates a substantial pressure 
fluctuation at this surface. These pulsations of 
negative pressure, estimated (for an inviscid fluid) to 
be an order of magnitude greater than the rozzle 
pressure drop, should readily overcome most normal 
hold-down pressures, and thus lift the chips previously 
created by the mechanical bit action, even at depths 
where cavitation is suppressed. This phenomenon has 
been modeled as follows (Geroux and Chahine, 1984a.) 

Neglecting the fine structure of the jet's shear 
layer and assuming that the major features of the flow 
field are associated with the large structures, the 
excited jet is modelled with large vortex rings emitted 
from the nozzle orifice at constant time intervals 
which correspond to the frequency of the pulsing jet. 
Each of these vortex rings is allowed to move and 
change shape under the influence of its self- induced 
velocity as well as the velocities induced by the other 
ring vortices present in the field. The presence of 
the wall is accounted for by applying the method of 
images since no viscous effects are considered. 

Knowing the location and characteristics of all rings 
at any time allows the computation of the generated 
pressure field. 

For a bubble ring of initial section radius, R 0 , 
stall as compared with the initial overall ring radius, 
A o , the self-induced velocity of the ring is 


v 


r 


self 4 tiR c 


- [bog 
A 


8A 


" 2 + 


■] + 0( f) 


( 10 ) 


W s ~ 1 is the inverse of fehe Weber number and is 
the circulation number. They measure respectively the 
relative importance of the surface tension arri the 
pressure drop due to circulation relative to the 
ambient pressure. These are written as: 

w s _1 = -r/(R 0 PJ , (11) 

fig = p( r/2*R 0 ) 2 /P„ ■ (12) 


The other terms in (10) are defined as follows 

e = Ro/A 0 «1 , A = A/A o . (13) 

A similar expression is known for vortex rings with a 
viscous liquid core. In that case e is replaced by 6, 
the ratio between the core size and the torus radius, 

A, the surface tension term is null, and the constant 
1/2 is replaced by 1/4 (Safman, 1970) 

V . =q-4- « [Log^ - + 0( 6) (14) 

VISCOUS 4 ttR o 4 J 
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ds 


S is any surface limited by the circular vortical line 
representing the vortex ring. P is a field point on 
this surface and thus PM is the vector connecting P and 
M. is the unit vector on the ring axis. The 
translation velocity of vortex ring due to the incoming 
flew generated by the other vortex rings is equal to 
the fluid velocity at its location. 

Figure 7 shows an example of vortex trajectories 
obtained numerically. T\vo distinct trajectories for 
odd and even number rings are seen. As a result, the 
even number rings get very close to the wall and 
therefore generate great tangential velocities and 
pressure drop on the wall along their trajectory. 
Comparison with high-speed movie observation of a 
STRATOJET ring shows a fair qualitative agreement. 



Figure 7 - Trajectories of Structured Jet Vortices, 

S = 0.25 

Figures 8 and 9 show the tangential velocities and 
pressures induced in five different locations on the 
solid wall versus time. The spikes correspond to the 
passage of a vortex ring. It is interesting to notice 
that both the amplitude and the frequencies of the 
locally sensed velocity and pressure spikes are 
dependent on the point of observation. The occurence 
of the pressure and velocity spikes are very well 
correlated as well as their amplitudes. Both seem to 
indicate that the predominant factor is the distance at 
a given time between the ring and the wall. This 
explains why the highest negative pressure peaks are at 
an optimun value of X. Both the high intensification 
of the tangential velocities and of the suction 
pressures on the wall due to the passage of organized 
structures indicate a significant enhancement of the 
cleaning capability of water jets through excitation. 
Ibis is possible through an increased shearing action 
and lifting forces on the solid chips at the wall. 


The interaction of different vortex rings is obtained 
usirg the velocity potential induced by each vortex 
ring at a field point M: 
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Figure 8 - Tangential Velocities Induced at Different 
Locations, X, on the Impacted Wall Versus 
Time 
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irrotational , one has to solve in the liquid the 
Laplace equation subjected to the boundary conditions 
at the torus wall, at infinity and eventually at a 
nearby target. The presence of the solid wall is taken 
into account by the use of the method of images: the 

wall is replaced by a symmetrical singularity 
distribution. 

In a matched asymptotic approach, the two boundary 
conditions at the bubble wall - namely the kinematic 
and dynamic conditions, apply only in the "inner 
problem," vhile the at-infinity condition and the solid 
wall boundary condition are to be verified only in the 
"outer problem." The two problems are linked with a 
matching condition. 

In order to simplify the analytical formulation, 
we follow the bubble during its motion, and we subject 
the frame of reference to move with the vortex. 
Therefore, the origin of coordinates in the inner 
problem is also the location of the vortex center. The 
first order equation for the normalized cross-section 
radius, R° = R°(t)/R 0 , of a bubble ring interacting 
with a solid wall initially at a distance dp = d(/A 0 
from the ring center is the following: 



In this expression, T is the characteristic time of the 
bubble dynamics. P„(t) being the time dependent 
imposed anbient pressure of characteristic amplitude, 
AP, the parameters in (16) are defined with: 

u 0 = d-ki 0 r 1/2 , pjt) = [pjt> - pjo) ]ap, 

P = P /AP , w e = RqAP/t, « = P( r/2nR 0 ) 2 / AP , ( 1 7 ) 

W = (p (0) - p )/ap = p + si 2 /2 - w _1 . 

^ » v ' go e 


Figure 9 - Pressures Induced at Different Locations, X, 
on the Impacted Wall, Versus Time 

c . E ffect of Jet Structuring on Cavitation Ero sion 

The aspect of the shear layer of a submerged 
cavitating jet is modified when excited from scraggly 
randan quasil inear vapor vortices to structured 
we 11 -organized ring vortices moving along the jet 
periphery. Figure 10 shows a comparison of cavitation 
appearance in three nozzle types. The self-resonating 
nozzle is seen, even at a higher cavitation number, to 
heavily cavitate. Figure 11 shows a high-speed film 
strip of the bubble ring motion. The behavior of these 
bubble rings as they collapse on a target is 
responsible for improved cavitation erosion when using 
the STRATCdET. We have modelled the dynamics of these 
bubbles when the ratio, e, between the characteristic 
size, Rq , of the cross-section of the ring and the 
overall ring radius, Ag, is small (Chahine and Genoux, 
1983; Genoux and Chahine, 1984b.) 

The existence of two length scales R 0 and A 0 
allowed the use of the method of matched asymptotic 
expansions. In the "outer region" (macroscale A 0 ) the 
torus appears as a distribution of singularities on a 
moving circle of radius A(t). In the "inner region" 
(microscale Ro) the torus is reduced to a cylinder of 
variable section R(t). At high orders of approximation 
deviation from the circular symmetry, motion of the 
vortex ring and modification of its overall radius, 
A(t), also appear. Assuming the fluid inviscid and the 


Pq „ is the initial gas pressure in the bubble, K the 
polytropic coefficient of compression, and P v is the 
vapor pressure. The characteristic time T is obtained 
by balancing the predominant terms between the left and 
right hand side of the dynamic equation (16), 

T = R 0 [pAP -1 bog(8A) ] 1/2 ; d 0 =0(1) , (18a) 

This result is the same as for a torus collapsing 
in an infinite medium since it assumes that atf<( a 0 ) is 
of order 1. In a more general case where d 0 is small 
compared to A 0 , the characteristic time T increases 
when d 0 decreases. Its expression can then be written, 
provided that d 0 «1, 

T Q = R 0 [pAP _1 Log (32/ed 0 ) ] 1/2 ; d 0 « 1. (18b) 

Equation (16) as well as higher order equations 
describing bubble ring motion and deformation can be 
solved numerically using a Runge-Kutta procedure. This 
was achieved for a large set of values of Pg 0 , e, ft, 
and d 0 . (Chahine, Genoux, Johnson, and Frederick, 
1984). The most important conclusions are the 
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SMITH NOZZLE 


Pa = 2 8 psi 

AP = 90 psi 

o “ 0.31 



CAVIJET NOZZLE NO LIP 


Pa = 27 psi 

A P = 102 psi 
0 “ 0.26 


CAVIJET NOZZLE THICK 
LIP 

Pa = 90 psi 

A P = 72 psi 

o » 0.56 
f ~ 700 Hz 


FIGURE 10- COMPARISON OF CAVITATION APPEARANCE IN THREE NOZZLE TYPES 
(ONLY THE LAST IS SELF-R ESON AT I NC ) 





FICURE 11 - 


HICH SPEED MOVIE SEQUENCES OF RINC MOTION AND 
EMISSION FROM A STRATOJET 
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following. The intensity of the collapse of a babble 
ring is of the same order as a spherical bubble but is 
spread on a larger area. Therefore, a higher erosive 
power is expected with bubble ring collapse However, 
the efficiency of this implosion depends 1 ... :y on the 
circulation parameter ft. Let us restate that SI is a 
measure of the amount of the pressure drop due to the 
vortical motion relative to the characteristic 
amplitude of dynamic pressure variations. Figure 12 
ilustrates the importance of Si. Increasing the 
vortical motion is seen to have the same effect on the 
torus behavior as increasing the amount of 
noncordensible gas in ' the bubble. 



Figure 12 - Vortex Strength Influence on Bubble Ring 
Collapse 

This result is logical once we consider Equation 
(16). It is seen that changing Pg or ft has a 
comparable effect, and for K = 1 (isothermal 
compression) the two parameters, P„ and ft, can be 
ccmbined into a single variable. Physically this is 
urderstardable since increasing ft implies decreasing 
the local ambient pressure at the torus boundary and 
thus reducing the driving collapse pressure. The same 
remark applied for Pn Q , since increasing the initial 
gas pressure, while keeping the ambient pressure 
constant, will also reduce the driving pressure for the 
implosion. This observation on the influence of ft has 
important consequences. For example, when attempting 
to increase the erosion caused by a jet, increasing ft 
has two opposite effects: (a) cavitation inception will 
ap>pear earlier, and (b) the collapse is damped. This 
indicates the existence of an optimum ft for given 
conditions and the need for an in-depth examination of 
the criterion required to determine this optimum 
value. The research of this optimum constitutes our 
main effort toward the design of an optimized STRATOJET 
device. 

In the presence of a solid wall, numerical 
solutions have shown the possibility of the formation 
of a reentering jet as observed for an initially 
spherical bubble. We have not as yet observed 
experimentally the existence of such a "conical" jet 
due to the very short time duration of the phenomenon 
and to camera frame rate limitations. Figure 13 shows 
the numerically generated bubble ring sections at 
different times during the collapse. Three cases are 


shown: a bubble ring in an infinite medium ana a 

bubble ring at distances of the same order as the 
overall ring radius and of eV2 times this radius. 

It is interesting to notice that during the collapse 
stage, there is a change in the direction of the 
reentering jet as the bubble approaches the wall. 


I** im trt.a.i 





c. d = O ( \J R o A o ) 


Figure 13 - Collapse of a Bubble Vortex Ring Near a 
Wall. Cross Section Shape at Successive 
Times 
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In the infinite medium case, the formation of this 
asymmetry in the bubble shape is only due to the vortex 
bubble self-induced motion, i.e., its translation 
notion relative to the ambient fluid due to both > 
vortex and the source distribution on the ring oemier 
line. In this case, the reentering jet is directed 
towards the ring axis in the direction of notion of the 
bubble ring vortex. In the presence of the wall, an 
opposite force is generated by the bubble ring image. 
Its importance increases with the_closeness to the 
wall. For instance, in the case d 0 = 0(e'/ 2 ), the 
force exerted by the wall overtakes the self-induced 
force, and the reentrant jet is directed away from the 
ring axis (i.e., the jet axis for the case of a _ 
submerged self- resonating cavitating jet). Por d 0 = 
0(1), the wall influence is weak, and the same behavior 
as in an infinite mediim is observed. Hie comparison 
of Figures 13b and 13c leads to the speculation that 
there is a distance, d 0 , for a given circulation for 
which the reentering jet of the collapsing ring is 
directed perpendicular to the wall. Since reentering 
jets constitute a basic erosive mechanism of cavitation 
bubbles, there is a practical application for this 
observation. 

Actual erosion tests of STRATOJETs versus 
conventional nozzles show a significant intensification 
of cavitation erosion with the jet oscillations. Table 
1 shows some typical results in Indiana Limestone for 
pressure drops across the nozzle of about 2500 psi. 
Except at the very low cavitation number of 0.067, the 
organ-pipe CAVIJET nozzle system cut this limestone to 
a depth 2 to 2.5 times greater than the standard 9nith 
system. Figure 14 shows similar results in aluminun at 
a nozzle pressure drop of 7000 psi. After a 5-minute 
jetting period, the STRATOJET does much better then 
both conventional CAVIJET and SMITH nozzles. 

NOISE GENERATION APPLICATIONS 



Figure 14 - Comparison of Mean Depth of Erosion for 
Four Nozzles in Alunintm 



where the jet acoustical pressure coefficient a is 
given by 

a = (n/2)S d (v'/V). (21) 

v'/V is the relative amplitude of the velocity 
oscillations, and Ap is the pressure drop across the 
nozzle. Similarly, the acoustic efficiency n of the 
jet, defined as the ratio of the radiated sound power 
to the hydraulic power, can be written: 

n = 8 q 2 M, (22) 

where M is the jet Mach number, M = V/c. 


Noise from submerged cavitating jets has been 
extensively studied. The early study of Jorgensen 
(1961) covered most of the ground to tie investigated 
and remains highly relevant to any study of cavitating 
jets. Jorgensen studied the intensity and frequency 
content of the noise generated by a submerged jet, 
investigating the scaling effect associated with nozzle 
size and the influence of the cavitation number. The 
acoustic energy upstream was found, for all cases 
investigated, to peak at a preferred normalized 
frequency or Strouhal number Sy. The value of 
was found to be strongly dependent on the cavitation 
number o. This "natural tendency" of a submerged jet 
to organize in large structures, and thus of its 
acoustic frequency content to peak near a discrete 
normalized frequency, is taken advantage of in the 
STRATOJET. The self-excitations enhance pressure 
oscillations in the Organ-Pipe feed tube and the 
formation of large structures in the jet shear layer. 
The nozzle-orifice area can be considered as a piston- 
like noise radiator, and the acoustical pressure p' 
generated at a distance r from the radiator can be 
approximated by 


p' = pcv'K d 2 /8r, (19) 

when the wavelength is large compared to the piston 
diameter. In Equation (19), v' is the velocity oscilla- 
tion of the radiator, K is the wave number of the 
oscillations, c is the speed of sound in the fluid, p 
is its density, and d the jet diameter. This can be 
written 


Equation (22) contains the basic information 
needed to optimize such a noise source. The highest 
coefficient of proportionality between p'/# 5 and d/r is 
needed as well as the highest possible jet velocity (or 
Mach number). These two factors are, in fact, 
dependent since maximizing the nozzle oscillations 
corresponds to a determination of the Mach number (see 
Equation (5)). Therefore, criteria for compromise 
depending on the sound generation requirements and on 
the physical constraints must be sought. 

Maximizing the nozzle oscillations also depends on 
maximizing the forcing function arising from the bubble 
ring collapse described in the preceding section. A 
measure of the cavitation amplification factor, A cav , 
is the ratio between the pressure generated by a bubble 
ring collapse and the pressure fluctuation, p' . One 
can show that A cav is bounded by (Chahine, Johnson, 
Lindenmuth, and Frederick, 1985) 


A > R 2 

cav 2 a 


(23) 


When the jet is excited and a is dramatically 
increased, A^y is reduced. This result reflects the 
fact that a part of the cavitation energy is then 
transferred into jet oscillation. Since the pressure 
driving the collapse is the ambient pressure around the 
jet, while the pressure scaling the jet noise is tp, 
the relative importance of the cavitation noise is 
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inversely proportional to the degree of cavitation. 
Fortunately, a lower bound for Ac av is proportional 
to R 2 . Maximizing R 2 is a logical requirement to 
maximize Piq^, since it means an increase.! the 
collapse intensity. Numerical values of P . order 10 
were obtained. This would indicate for the 
experimental tests reported below a significantly high 
cavitation amplification factor, Ac av > 125. 

In order to evaluate the STRATQJET noise 
generation characteristics, minimize reflections, avoid 
standing waves in the range of frequencies of interest, 
and be able to vary o for constant jet velocity, we 
used the HydronauticS-Ship Model Basin. The basin is 
128 m long, 7.6 m wide, and 4 m deep. In addition, the 
presence in this basin of a well 13 m deep and 3.35 m 
in diameter allowed us, by locating the nozzle at 
various depths, to control the ambient pressure at the 
nozzle exit. This enabled the variation of the 
cavitation number while keeping the jet velocity 
constant without the drawbacks of a small enclosed tank 
and wall reflections. 

The nozzle assembly tested for noise generation 
consisted of a plenum with an 8.25 an diameter feeding 
a straight stainless steel pipe of length L and 2.5 cm 
diameter. This pipe, which acted as an organ pipe, 
discharged into the ambient liquid through a nozzle 
having a hemispherical shape with a sharp edge. The 
nozzle lip was shaped and sized such as to provide a 
feedback of the preferred natural frequencies of the 
jet oscillations to the organ pipe which amplified 
them. That lip consisted of a cylindrical section of 
length e, followed by a conical bevel of length t. A 
pressure transducer, located at the midpoint of the 
organ pipe, allowed for comparisons between the 
pressure oscillations in the pipe and the noise 
measurements in the field surrounding the jet. An 
example of the correspondence between the signal of the 
transducer located in the pipe and the hydrophone 
signal can be seen in Figures 15 and 16. A good 
agreement was observed at the low frequencies, 
especially at the main organ-pipe frequency. However, 
the high frequencies sensed by the hydrophone which are 
probably due to cavitation were filtered by the 
organ-pipe chamber and hence were not strong in the 
pipe. 
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Figure 15 - Frequency Spectrun of Pressure Fluctuations 
in Pipe and of Acoustic Pressures at r = 
0.92 m, 8 = 18', Hydrophone 3.4 x 10 4 
Pascals/volt, Transducer 1.4 x 10 5 
Pascals/volts, d = 1 cm, AP 0.35 MPa 



Figure 16 - Frequency Spectrum of Pressure Fluctuations 
in Pipe and of Acoustic Pressures at r = 
0.92 m, 8 = 18’, Hydrophone: 1 volt = 1.4 

x 10 4 Pascals. Transducer: 1 volt = 1.4 x 

10 5 Pascals/volts, d = 1 cm, AP « 0.35 MPa 

For a practical noise generator it is of great 
importance to determine the decay function of the 
generated noise with distance frcm the nozzle. Figure 
17 shows, for two nozzles of diameters 1 cm and 0.5 cm, 
the results obtained at distances of up to 500 nozzle 
diameters. These results give the most important 
indication of the effectiveness of the STRATCJCT as a 
noise generator. Both the rms value of the acoustic 
pressure and the amplitude of the fundamental frequency 
of the noise detected by the hydrophone decayed with 
the distance r frcm the nozzle as 1/r, with a slope a 
{Equation 20) . 

The measured value of a for the conditions 
presented in Figure 17 was for 9<60°: 


a =0.18, for d=1an, 
a =0.11, for d=0.5 cm. 


(24) 
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Figure 17 - Attenuation of the RMS Value of the 

Pressure Fluctuations with Normalized 
Radial Distance frcm Nozzle Orifice 
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The difference between these two values was 
related to the acoustical performance of the two 
nozzles with relation to the cavitation nunber. '!' > 
acoustical performance of the self-resonating nozzle 
was much better than a nonresonating cavitating jet. 
Expressed in terms of the coefficient a, such 
non-resonating jets would have the value (Jorgensen, 
1961) 


0.008 < oK 0.02. (25) 

The acoustic efficiency, n, for these STRAIOJETs is: 
n a 0.268 M, for d=>1cm, 


(26) 


n=0.093 M, for d=0.5 cm. 

Jorgensen reports for a conventional cavitating nozzle: 


0.0005 M < n < 0.003 M. (27) 


These values for conventional cavitating jets, even if 
an order of magnitude higher than noncavitating jets, 
are two to three orders of magnitude lower than the 
efficiency of the STRATQJETs used in these 
experiments. In terms of decibels, the 1-cm- diameter 
nozzle, working at a pressure drop of 0.35 MPa, 
generated about 176 dB re: 1 kiPa at a distance of 1 m. 
This value can be increased by increasing nozzle 
diameter, pressure across the nozzle, and by improving 
the nozzle efficiency with better design. Limitation 
on this increase would, however, come frcm scaling 
effects and practical limitations on geometrical size 
and punp power. Recent tests have al90 shown that the 
generated signal could be easily modulated with 
relatively high frequency (more than half the main 
frequency for these tests) . 

NUMERICAL SIMULATION OF JET STRUCTURING 


The understanding of the basic mechanisms by which 
self-resonance and intense oscillations of the flow are 
initiated and sustained is essential for further 
development and use of these new nozzle systems. Our 
first efforts have used systematic experimentation to 
determine the satisfactory nozzle shape and piping size 
which achieve strong oscillations at predetermined flow 
conditions. Problems appeared later when scaling 
effects, different liquids, (e.g., drilling mud), and 
lower cavitation number conditions had to be dealt 
with. This has justified our concentrated effort on an 
analytical and experimental program to improve our 
knowledge and establish general guidelines for the 
optimization of the nozzle/piping assembly design for 
any prescribed conditions. 

Many investigators have theoretically examined jet 
structuring, and, based on experimental observations, 
severed attempts have been carried out to model a jet 
with a discrete distribution of vortex rings. Axially 
symmetric vortex rings with a visoous core have been 
investigated, and the numerical schemes were able 
satisfactorily to reproduce experimented measurements 
or observations. Ft>r instance, Acton (1980) 
represented the upstream portion of the jet with a 
double row of ring vortices which were allowed to 
interact and move freely downstream of the nozzle 
exit. The roll ip of the vortex sheet into large 
axisymmetric structures was numerically demonstrated. 
Artificial excitation of the jet by periodically 
modulating the strength of the emitted vortices 


exhibited the same features as observed experimentally, 
namely the existence of an optimal Strouhal number. 

In an ongoing program, we follow the same approach 
and introduce, in addition, a potential flow model for 
the flow inside the feed section to the nozzle 
orifice. The interaction between the free emitted 
vortices and the nozzle walls in the separated area of 
the nozzle (Figure 18) is allowed. The interaction in 
this region between the jet shear layer and the walls 
controls the feedback mechanism, and its study is 
fundamental for the understanding and the control of 
this mechanist. 

As sketched in Figure 18, the flow field under 
study is divided into three regions, a) The first 
region, noted U for ipstream, is constituted of a 
semi-infinite cylinder. The flow in that region is 
reproduced by a cylindrical vortex sheet of uniform 
intensity wyp. b) The second region is noted N and 
corresponds to the nozzle walls both in the wetted and 
separated areas. In that region the flow is reproduced 
by cylindrical distributions of vortices of intensity u 
and of sources of intensity s. c) The last region is 
the jet region and is noted S for the shear layer. The 
jet shear layer is replaced by a succession of discrete 
rings emitted from the nozzle, of spacing X, which are 
free to move and interact together and with the 
upstream singularities. The circulation of any ring k, 
r(k) can be modulated, oscillate with time, and be tied 
to the organ pipe oscillations. 



Figure 18 - Schematics of Nozzle Assemblies for Flow 
Visualization and Numerical Simulation 

The distribution of singularities in region U is 
constant with time and constitutes a "boundary 
condition" for the rest of the flow. However, the 
intensity of the singularities in region N constantly 
varies with time and is such that there is no 
cross-flow at the nozzle walls. This intensity is 
mainly controlled by the positions of all freely moving 
rings, and by the influence at the considered point of 
the rest of the singularity sheet. Similarly, the 
motion of each vortex ring is a combination of its 
self-induced velocity, the velocity induced by the 
other rings, and the velocity due to the singularity 
distributions in regions U and N. 

For the numerical simulation, the solid boundaries 
are discretized into axisymmetrical panels of small 
size on which the distribution of singularities is 
uniform. For a given nozzle shape, the flow induced by 
a distribution of sources and vortices on the panel "i" 
of the nozzle is conputed at any discrete point "j". 
This allows the definition of matrices of influence 
which are a function only of the nozzle geometry and 
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independent of the rinqs emitted at the exit. This 
simplifies considerably the numerical procedure and the 
computation time. These matrices and their inverses 
are computed once and for all for a given nozzle and 
then become inputs for the actual flow probl 
Outside the nozzle, the ring-vortices emitted interact 
together and with the upstream singularities in the 
same manner as described in preceding sections. 

The computations are done in the following 
sequence (Chahine, Genoux, and Liu, 1984). TVo 
matrices of influence, defined as follows, are 
computed. The matrices, and are such that the 
components AS^j and AVjtj are the normal velocities 
induced at the location of an element E[ by a uniform 
distribution of sources and vortices of unit intensity 
on the element Ej . We define, in addition, the 
vectors, AVy ,, and V TO , whose components are the 
normal velocities at each element E^ due respectively 
to a vortex sheet of unit intensity in the region U and 
to the discrete ring vortices in region S. AV, AS, and 
AVjp are not time-dependent and are computed once for 
W given nozzle geometry. However, has to be 

computed at each time step. 

With the above definitions, the condition of zero 
normal velocity at the nozzle wall can be written: 


0 




s + A V • 


OJ + 


0) AV 

Up —Up 


+ - TO' 


(28) 


where the vectors s and u have as ccmponents the actual 
unknown source and vortex intensities at each element 
E;. The system fas one degree of freedcm, and s or w 
could be imposed. If we preselect the value of w, then 
the source intensity distribution is obtained by the 
following expression: 

s = - ^T 1 - V AV up+ V, m ). (29) 


Since changes at each time step, s is recomputed 
accordingly. We have implemented such a computation 
code, and two test cases are presented in Figures 19 
and 20. Figure 19 shows the case of a Smith nozzle 
with a cylindrical whistler mouth added to it. The 
time evolution of the ring position, the strong 
fluctuations in the whistler, and the tendency to 
structuring can be seen. However, large computation 
times are needed to introduce enough ring elements to 
observe distinct structuring. Figure 20 shows a case 
where t>>1. Structuring is here very clear. In order 
to use the code as a reasonable tool for resonating 
nozzle design, we are presently working on means to 
reduce the computation time. This involves automated 
schemes to recognize clustering of vortices in a large 
structure and the replacement of this structure with a 
single high circulation vortex. 

CONCLUSIONS AND ACKNOWLEDGMENTS 

In this paper we have described seme aspects of 
the mechanics and principles of operation of 
self-resonating cavitating jets. We have tried to 
highlight the principal advantages of these jets, 
namely, early cavitation inception, improved cleaning 
capabilities and chip removal, enhanced cavitation 
erosion, and intense noise generation. We have tried, 
within the space limitations, to describe our 
contributions to the modeling of the various aspects of 
these jets. We ask the interested reader to consult 
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Figure 19 - Simulation f Jet Nozzle Interaction. Smith 
Nozzle with Whistler Mouth - Oscillating 
Vortex Strength 
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Figure 20 - Simulation of Jet Shear Layer Structuring. 
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earlier publications on the various details. Several 
aspects of the mechanics of these jets are still 
unclear and our ongoing efforts are concentrating on 
these aspects. FOr instance, the fundamental 
understanding of the feedback mechanism at the nozzle 
orifice, viscous effects, the optimization of ring 
vortex circulation, and the influence of the cavitation 
number on self-oscillations are important questions 
which need investigation. 
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ONE-LITTER AMBULANCE for frontline use is mission envisaged for Para-Copter rotcrcraft by inventor, George Schmidt. 
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SCHMIDT PULSE-JET ENQINE 


“ALMOST SILENT” PULSEJET ENGINE used by Schmidt delivers 105 lb. of thrust on special kerosene-base fuel. 


Noise Is Cut in New Pulsejet Engine 


A small pulsejet-powered helicopter, 
designed and built by George Schmidt 
as a combat evacuation craft, is nearing 
completion at the former Focke-Wulf 
engineer’s Brentwood, N. Y., factory. 

Military interest has been shown in 
the Para-Copter 2 (Aviation Week 
Sept. 7, p. 9), intended as a one-litter 
ambulance to fly wounded soldiers from 
front line areas to points of safety. 

► “Silent” Operation— Schmidt engi- 
neered and constructed the 105-lb. - 
thrust twin pulsejet engines for his 
new craft, developed a special kerosene- 
base fuel called JPX, and devised a 
silencing system that he claims will 
make the helicopter powerplants ap- 
proximately 80% quieter than other 
existing types. 

“Almost silent” operation of the . 
machine is projected as a top selling 
point for military operation and in the 


civilian market, if the engineer develops 
a civil configuration of his copter. 

► Adjustable CG— Schmidt based the 
Model 2 on an earlier, smaller strap-on 
type of helicopter that he built and 
flew successively. 

The Para-Copter is constructed of 
welded steel tubing, stands 8 ft. 7 in. 
high, and has a folding tricycle landing 
gear of 7-ft. base. 

An aluminum seat, located directly 
under the rotor hub, is adjustable fore 
and aft for center of gravity balance. 
The single litter is carried directly un- 
der the pilot seat. 

► Little Torque— Schmidt designed ro- 
tor controls that extended from the 
rotor hub but now is installing conven- 
tional floor configuration demanded by 
the latest military specifications. 

Practically no compensation for 
torque is required in the Para-Copter’s 


control surfaces, the designer says, be- 
cause the only torque existing in the 
craft is caused by bearing friction. A 
small fin may be added to correct the 
friction effect, however, if the need is 
indicated by future flight tests. 

► Powerful Pulsejet— Schmidt developed 
several innovations in design, operation 
and mounting of the helicopter pulse- 
jet engines. 

Each powerplant is 32 in. long, 5 
in. in major diameter and weighs 17.5 
lb. Although smaller than the new 
British pulsejet developed by Saunders- 
Roe for copters (story on next page), the 
Schmidt engine delivers more than 
twice the former’s 45 lb. thrust. 

The engine shell is constructed of 
welded 347 stainless steel. 

► Menthol Mixture— To preserve the 
life of four valves in each pulsejet, a 
mixture of menthol crystals dissolved 
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in alcohol and mixed with water is 
bled from a high pressure tank into a 
fuel line mounted just ahead of the 
throttle regulator. One pint of this 
coolant is required for each hour of 
engine operation. 

In operation, a small amount of 
coolant mixed with fuel is sprayed into 
the pulsejet just ahead of the venturi 
through a ,004-in. hole. The same fuel- 
coolant mixture is injected into an 
expansion chamber of the engine 
through a No. 60 hole. 

This combination, Sclunidt claims, 
keeps valves much cooler, preserves 
their life and gives a better fuel mix- 
ture for operation. 

► Muffled Pops— The Para-Rotor’s si- 
lencing system probably is the first to 
cut down noise of a pulsejet engine. 

Schmidt deadens the loud, pulsating 
pop in two operations: 

• Nozzle of .007 in. diameter sprays 


► Simple Overhaul— Major overhaul of 
the Schmidt pulsejets after 60 hr. of 
operation requires less than an hour 
Or each powerplant: 

• Nose venturis are removed by tak- 
ing out four screws. 

• Valve plate assembly is lifted out. 

• New assembly is installed. 

Unlike conventional rotor-tip jet 
units, the Schmidt engines are not at- 
tached solidly to the rotor blade but 
arc installed on a ball-bearing mount 
that permits the powerplant to swing 
freely within limits in an up-and-down 
motion. 

Small fins at the tail of the pulsejets 
and the center-of-gravity mounting 
keep the engines in a horizontal posi- 
tion regardless of the angle of attack 
changes of the rotor blades. 

A manufacturer has quoted Schmidt 
a price of $85 per pulsejet for the first 


lot o£ 50 powerplants. 

► Strap-On Copter— Schmidt flew his 
early strap-on helicopter to an altitude 
of nearly 12,000 ft. in July of last year. 
The craft was demolished on a second 
flight when the designer, unfamiliar 
with the flying characteristics of his 
new copter, allowed the rotor-craft to 
be blown sideways and tipped over bv 
a strong gust of wind. 

But Schmidt says he expects to re- 
build the small copter. 

► Achgelis Engineer— The German-born 
engineer worked on Focke-Wulf’s 
Achgelis helicopter before he came to 
the United States in 1939. He has been 
employed by several major U. S. air- 
craft firms during the past 14 years. 

Working with Schmidt on the Para- 
Coptcr arc T. McGeough, chief engi- 
neer of the independent company, and 
engineer M. Simbrow. 


water into the ejection chamber, where 
the 1,700F heat turns it into a steam 
mass that muffles noise. Each engine 
uses approximately one quart of steam 
an hour. 

• Thirty-two tubes of 7-in. long stain- 
less steel are mounted in the tail of 
the pulsejet ejection chamber, where 
they change noise that escapes the 
steam muffle into a high-pitched sound 
above the level of human reception. 

Schmidt claims Ins silencing sys-^T 
tem also gives the pulsejet 15-20WL 
greater thrust. 

► Folding Package— Rotor blades on the 
new copter have a straight chord of 
9 in., arc hinged to fold at the semi- 
rigid, sce-saw hub. 

When the rotors and landing gear 
are folded, the Para-Copter is a small 
package that can be stowed aboard a 
cargo transport and dropped bv para- 
chute into forward combat areas. 



The Clark Y airfoil of the copter’s 
rotor blades, made of eight layers of 
laminated sitka spruce and seven layers 
of nylon, may be shaped and mass- 
produced by unskilled workers on a 
special machine designed by Schmidt. 

A 7-ft. blade can be cut in six 
minutes on this machine, he reports. 
► Pressurized System— Fuel system of 
the new helicopter is a pressure type 
up to the rotor head, where it becomes 
centrifugal. 

Rotor blades currently have integrally 
built fuel lines, but they will be re- 
movable in the future to facilitate 
maintenance. 

An 18-gal. fuel tank, located behind 
the pilot seat, gives the Para-Copter 
sn endurance of 1 hr. 10 min. Air pres- 
sure of 35-40 lb. in the tank is used 
to start the pulsejets, but centrifugal 
forces keeps the fuel flowing after the 
engines begin burning. 

A bottle containing 92 cu. in. of 
compressed air is carried for starting 
purposes. 


British Pulsejet 1 

The pulsejet engine is the most 
promising means for tip propulsion of 
rotors, says the Helicopter Division of 
Saunders-Roe, Ltd., British aircraft firm. 

After a year of experimenting, the 
company has developed a 45-lb. -thrust 
unit that has logged many hours of 
test and is working on a 1 20-lb. -thrust 
engine. 

► First British Pulsejet— This is the first 
British pulsejet to propel or lift any air- 
craft, although similar U.S. engines 
have powered American Helicopter’s 
XH-26 and Marquardt’s M-14 for some 
time. 

The pulsejet, one of the family of 
aerothermodynamic engines that in- 
cludes the ramjet, was used by the Ger- 
mans during World War II to power 
V-l buzz-bombs. 

► Operation— A pulsejet engine consists 
of simple parts: A cylindrical combus- 


3 Power Copters 

tion chamber with a spring-loaded, nor- 
mally open intake valve, and a tailpipe. 

In operation, a mixture of fuel and 
air is ignited in the combustion 
chamber. The resulting explosion forces 
gases out the exhaust pipe at high speed 
and closes the intake valve. 

The discharge lowers combustion 
chamber pressure to a point where the 
spring-loaded valve is able to open 
again, and a fresh charge of air is sucked 
and rammed in. The cycle then repeats. 
► Acoustic Tuning— The Saunders-Roe 
unit is “tuned” by the acoustic prop- 
erties of the chamber and tailpipe to a 
frequency of 120 cycles per second. It 
can be varied to a maximum of approxi- 
mately 200 cps. or a minimum of 46 
cps. 

Engine weighs only about 15 lb. for 
its 45-lb. thrust. It is 47.5 In. long, 
and 5.5 in. in diameter. 
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Recommended noise 
exposure limits are a major 
factor in the decision to 
institute a noise control 
program. Air equipment 
produces vibrations 
ranging from the 
subaudible (below ' Hz) 
to the limits of man 
perception. Depenc, g on 
the environment, noise 
control usually should be 
considered for equipment 
vibrating at frequencies 
over 30 Hz and producing 
noise above 90 dBA. 



Beating OSHA to the punch by 



Air compressors and vacuum pumps are prime 
sources of noise in offices and plants. With OSHA 
about to issue a standard limit-ing allowable 
exposure to 90 dBA, users and manufacturers 
alike are scrambling to find ways to quiet their 
pneumatic components. Generally, the 
technology to reduce noise in air-powered 
equipment is readily available and easy to apply. 


ROBERT MOFFATT 

Gast Mfg. Corp. 
Benton Harbor, Mich. 
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PNEUMATIC components in 
equipment as diverse as pack- 
ing machines and computer 
peripherals can produce noise 
exceeding daily exposure limits 
mandated by OSHA. Besides 
being dangerous to hearing, 
such noise can interfere with 
communications, contribute to 
lowered productivity, and cause 
frequent errors and inattention. 
Also, the vibrations that pro- 
duce the noise can damage 
machinery and pipes. 

Fortunately, noise generated 
by air systems is comparatively 
easy to control and can be kept 




How Noise Dissipates 



Noise exposure often can be reduced to 
acceptable levels by simply separating 
the noise source from personnel. This 
graph shows free-field dBA reduction 
produced by'increasing the distance 
between a noise source and a listener. The 
chart's usefulness depends on knowing 
the dBA level at 3 ft. For example , a 
piston-driven vacuum pump rated at 73 
dBA at 3 ft produces only 62.6 dBA at 10 ft 
and 54.7 dBA at 25 ft. 


well below OSHA limits. Al- 
though the burden of noise con- 
trol normally falls on the user, 
information on noise charac- 
teristics and noise control is 
readily available from man- 
ufacturers. The decision on how 
to control noise generally de- 
pends on a consideration of 
costs, operating environment, 
worker comfort, and legal re- 
quirements. 

When Is Noise a 
Problem? 

The determination of 
whether noise is a problem de- 


pends on how people react to the 
sound emitted. Many factors 
about noise are relative and 
subjective, a fact that should be 
taken into account when plan- 
ning a noise-attenuation pro- 
gram. For instance, a compres- 
sor operating in a metal- 
stamping plant might go un- 
noticed, but the same compres- 
sor in a dental office could be 
quite annoying. 

To establish whether a prob- 
lem exists at the objective level, 
the sound must be analyzed 
with a sound level meter or an 
octave-band analyzer. Such 


equipment is useful in deter- 
mining potential hearing dam- 
age, annoyance, and speech in- 
terference. 

Annoyance level and loud- 
ness may not be synonymous, 
although a loud noise is usually 
more annoying than the same 
sound that is somehow muted. 
Given the same sound pressure 
level, higher-frequency sounds 
(above 2,000 Hz) are usually 
more annoying than lower- 
frequencv sounds. 

The vibrations produced by 
air equipment range from the 
subaudible to the limit of 
human perception. A low-speed 
pump, for instance, can produce 
low-frequency vibrations (less 
than 15 Hz) that are felt rather 
than heard. As frequency rises, 
the vibrations are accompanied 
by noise, with the audible range 
beginning at about 30 Hz. In 
this range, noise control usually 
becomes an important consid- 
. eration. 

Two Ways To Reduce 
Noise 

Generally, noise is reduced at 
two locations: the source or the 
listener’s ear (by changing the 
sound path). The first correction 
involves installing equipment 
such as dampeners and mufflers 
at the intake or discharge end of 
a compressor, vacuum pump, or 
air motor. Or the equipment can 
be placed in a specially built en- 
closure. 

Noise reduction at the lis- 
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tener’s ear can involve simply 
moving either people or equip- 
ment to another location. The 
relationship between noise 
level and distance is shown in 
the noise reduction graph. 

The sound path also can be 
changed by modifying the 
acoustical environment or plac- 
ing appropriate barriers be- 
tween the source and listener. 
Sound follows a line-of-sight 
path; obstructing the path re- 
duces noise. The intent of these 




Oil and contaminant drain 


Mu r based on the phase-shift principle lowers frequency peaks of 
en: ig air by splitting the air stream into two branches and shifting them 

IS i, jut of phase. In operation, some air passes through the small holes 
in the first baffle plate, and some rebounds and enters the flared ends of 
the tubes. The tubes are of different lengths, so some air flows to the 
muffler exhaust and some bounces off the second baffle plate. The 
second air stream then enters another tube and flows to the exhaust. 


methods is not necessarily to si- 
lence the equipment, but to re- 
duce noise to levels that are 
physically and psychologically 
tolerable and to reduce vibra- 
tions to levels that pose no dan- 
ger to the equipment itself. 

There are many types of noise 
attenuation devices for air 
equipment. Dampeners and 
mufflers covering the range of 
pulsations from 1 to 20,000 Hz 
are available. Some are based 
on the phase-shift principle, 
which is analogous to cancel- 
ling out or stabilizing three- 
phase alternating current. 


Also, absorptive and dissipat- 
ing mufflers, in which the com- 
pressed air passes over perfora- 
tions that lead to sound- 
absorbing materials, can be 
used. 


One of the most economical 
and most effective devices to at- 
tenuate noise over a wide spec- 
trum is the low-pass filter. This 
filter, which can be as simple as 
a glass or plastic jar , allows only 
lower-frequency sound waves to 
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How Low-Pass Filters Work 



In a low-pass acoustic filter, the ratio of the 
cross-sectional areas, A\ and A2, to each other and to 
lengths, Li and L2, determines the frequency range of 
sound that is attenuated and passed through. Such a filter 
can be as simple as a glass or plastic jar mounted on a 
pump outlet. 



orifices acts as a high-pass fil- 
ter. 

The effect of a low-pass filter 
is the same as that produced by 
passing the sound waves 
through alternating layers of 
different media, each transmit- 
ting only certain frequencies. 
The power transmission 
ratio — the ratio of average 
power emerging from the 
medium to that entering 
it — can be calculated tc deter- 
mine the influence of the layer. 
This ratio depends on medium 
thickness and specific acoustic 
resistance for a plane wave. 
(For instance, water has a 
specific acoustic resistance 
3,800 times that of air.) Various 
types of media thus could be 
used to build an acoustic filter 
that either blocks or selectively 
transmits waves of certain fre- 
quencies. 

Compressors and vacuum 
pumps also can be modified by 
inserting a plastic-foam liner 
directly into the units. This can 
reduce pump and compressor 
noise by up to 14 dBA, depend- 
ing on the amount of pressure or 
vacuum being generated. Mod- 
ification of a unit with a readily 


available kit is easy; the end 
plate is removed and the liner 
simply slipped between the 
shroud and body of the pump or 
compressor. 

Acougtically insulated cab- 
inets are also available to 
reduce air-equipment noise to 
office levels (about 50 dBA). 
One version of such an enclo- 
sure quiets — by as much as 17 
dBA — the equipment by hold- 
ing it on a vibration isolator in- 
side a foam-lined, baffled 
cabinet. 

Compressors and vacuum 
pumps engineered to run 
quietly are, of course, the best 
choices for applications that re- 
quire noise levels even lower 
than those mandated by any 
"official” standard— in dental 
and medical offices, for exam- 
ple. One such unit reduces noise 
output by as much as 15 dBA 
because its motor is suspended 
on brackets within a plastic- 
foam-insulated canister. The 
unit also is fitted with rubber 
mounting feet. (The standard 
unit, by contrast, is rigidly 
mounted directly to the motor.) 
The quiet model also can be fit- 
ted with a double-length car- 


tridge muffler to reduce noise 
even further. 

Noise Trade-offs 

The technology and data re- 
quired to reduce noise in air 
equipment are readily availa- 
ble. The major criterion for se- 
lecting a method is to buy only 
as much quietness as is re- 
quired. Indeed, it may be more 
difficult to define required noise 
attenuation — regulatory, com- 
fort, and aesthetic considera- 
tions could have significant 
impact — than it is to select a 
muffling technique. 

The least expensive cure is to 
simply increase the distance be- 
tween the noise source and the 
listener. However, this may not 
be possible if space is already at 
a premium. 

Applying the various sound- 
deadening techniques costs 
somewhat more than merely 
moving equipment around. 
Minor modifications such as 
quiet boxes or kits are the least 
expensive, add-ons such as 
larger mufflers are moderately 
expensive, and specially en- 
gineered units are the most ex- 
pensive. 
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(word processor parameters LM=8, RM=75, TM=2, BM=2) 
Taken from KeelyNet BBS (214) 324-3501 

Sponsored by Vangard Sciences 
PO BOX 1031 
Mesquite.- TX 75150 

There are ABSOLUTELY NO RESTRICTIONS 
on duplicating, publishing or distributing the 
files on KeelyNet except where noted! 
February 26, 1993 

NOISECNC. ASC 


This file is from the Dallas Morning News of 
Saturday, February 13, 1993 - Business section. 


The following file is most intriguing and applies to one of the root 
studies of Vangard Sciences and KeelyNet . Keely used compound 
vibrations and we believe at various phasing differentials. This 
phase adjustment is known in modern terms as interference or phase 
conjugation and can be either constructive (additive) or 
destructive. Use of such principles can lead to phenomena beyond 
simple cancellation of a signal, but we will stay within the context 
of the article as given. 


Technology Report 
Silence is gold for companies 
Experiments seek noise reduction 
by Rebecca Smith 

San Jose, Calif. - Imagine summer without the earsplitting roar of 
power lawn mowers, flying without hearing the whining din of jet 
engines, running the dishwasher without having to shout over the 
noise . 

Companies that make noisy products are looking at an ingenious 
breakthrough technology that enables them to dramatically reduce 
unwanted sounds . 

"Active noise reduction" systems represent the broadest assault on 
noise pollution since invention of the muffler. Products employing 
the technology will begin showing up in mass consumer products later 
this year. 

In what could be a boost for U.S. competitiveness, about 75 patents 
covering product - speci f ic applications of the technology are 
controlled by two U.S. companies. 

"It'll be part of the next wave of luxury home appliances and cars," 
says technology analyst Dan Hutcheson, president of VLSI Research in 
San Jose. "Right now, it's expensive, but the cost should come down 
fairly fast as volumes build." 

The idea of active noise reduction dates back to a 1936 patent for a 
"process of silencing sound oscillations" by German inventor Paul 
Lueg . But commercial execution required invention of the microchip. 
The concept is simple. Fight fire with fire. 
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Here's how it works : 

Start with an irritating noise, say, that of an air compressor. 
A microphone samples the noise and sends a signal to a microchip 
analog-to-digital converter, which translates the sound waves 
into computer language. It creates a digitized blueprint of the 
sound, then passes the blueprint to a digital signal processing 
chip that analyzes the sound's magnitude and frequency. 
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The digital sound processor (DSP) then computes WHAT pattern of 
sound waves would be EQUAL but OPPOSITE - a MIRROR IMAGE. (180 
degree shifted for perfect cancellation ... Vangard) 

That sound profile (signature) is then converted back from a digital 
to an analog representation. An audible sound wave is created and 
amplified through speakers or a headset. When the yin-and-yang 
sound patterns collide - the whole process takes a fraction of a 
second - they largely cancel each other out . The loud air 

compressor noise becomes a quiet hum. Advanced systems continuously 
sample the ambient noise level, correcting and re-correcting. They 
follow noises as they change, such as engines accelerating. 

"Most of these systems are still in the experimental stage," says 
Mike Collins, a manager of digital signal operations at Motorola 
Corp . in Austin. "In a couple of years, it has the potential to 
become a commonplace consumer item that will make quieter 

dishwashers, quieter airplanes and so forth." 

Most of the patents required to make a DSP-based noise reduction 
system work are owned by two small U.S. companies : privately held 

Active Noise and Vibration Technology in Phoenix and publicly held 
Nose Cancellation Technologies in Stamford, Conn. They are working 
with Motorola and Analog Devices of Norwood, Mass. to develop the 
best chips. Other companies are expected to become active players 
in coming years, at least at the chip-production level. 

Active Noise says it will have headsets on the market this year 
carrying a price tag of about $250. They'll be marketed primarily 
to workers in noisy industrial settings, where they will replace 
earplugs as safety equipment. 

Unlike earplugs, which block sounds indiscriminately, the headsets 
will be programmed to target low-frequency, repetitive sounds that 
bother most people. A jackhammer operator, at least theoretically, 
will be able to operate his or her machine and still converse with a 
co-worker without shouting. As brilliant as noise reduction 
technology sounds, not everybody is getting on board yet. 

"The more mainstream approach is to try to reduce noise at the 
source," says engineer Tom Mitchell, head of the Noise Technology 
Center at General Motors Corp ' s Power Train Division in Pontiac, 
Mich. "Anything you add to a vehicle to susceptible to failure, and 
there's some added cost associated with a noise reduction system." 
Mr. Mitchell says GM has devoted considerable resources to reducing 
unwanted vibrations and noise in automobiles. Fluid- filled engine 
mounts, developed eight years ago, are an example of a simple device 
that eliminates noise by reducing the transfer of vibrations from 
the engine to the frame of the car. 
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"I could bore you for hours with all the stuff we've developed," 
says Mr. Mitchell. "Like double-wall, insulated exhause pipes and 
sound absorbers in the headliner. We devote a lot of energy to 
reducing noise, because customers tell us it's one of the best 
mechanisms they have for assessing quality." Reducing noise can get 
kind of tricky. Nobody wants road noise, thus tire makers have 
sound standards to meet . 

"On the other hand, nobody wants absolute silence. They want some 
feedback that the engine is there. That's where an electric car can 
be eerie," Mr. Mitchell says. 

The key to the noise-cancellation feature of DSP systems rests in 
the software that's built into the chips. That's where most 
research and development energies are being expended. The chip has 
to be told which sounds to attack and which to leave alone . 

"You don't want to cancel the sound of a great Porsche engine," says 
Tom Hesse, founder and chief executive officer of Active Noise in 
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Phoenix. In some cases, noise suppression systems actually do more 
than eliminate an irritation or a safety hazard. They can improve 
performance. Systems are being developed to replace the entire 
exhaust system and muffler, thus eliminating the back pressure that 
the muffler creates. 

"As a result, you get up to a 5 percent improvement in gas mileage," 
says John McCoy, president of Noise Cancellation Technologies. 

Many companies - including Boeing Co., Mercedes Benz and Whirlpool 
Corp . - confirm they are looking at noise reduction systems but 
don't want to comment on their research progress because they don't 
want to tip off their competitors about the timing of the 
introduction of improved devices. "I know there's a lot of work 
going on here, but it's highly proprietary," says Boeing spokesman 
Tom Cole . 

Meanwhile, designers at Noise Cancellation and Active Noise are busy 
drumming up interest. So far, most of their income has come from 
licensing big companies to develope their own systems for specific 
products. When they come to market, there may be some interesting 
twist s . 

"We could use it to do more than reduce noise . It could create 
sounds," says Active Noise's Mr. Hesse. "We probably could make a 
Honda Civic sound like a Ferrari . " 

Better yet, how about one that can tune out, "Are we there yet?" 


If you have comments or other information relating to such topics 
as this paper covers, please upload to KeelyNet or send to the 
Vangard Sciences address as listed on the first page. 
Thank you for your consideration, interest and support . 

Jerry W. Decker Ron Barker Chuck Henderson 

Vangard Sciences/KeelyNet 


If we can be of service, you may contact 
Jerry at (214) 324-8741 or Ron at (214) 242-9346 
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